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Editors’ foreword 


This book fits closely into the network of topics 

of the Oxford Physics Series and at the same time 
preserves an integrity of its own. The earlier con- 
plementary core texts include Radtatton and quantum 
phystes, and Atomic nuclet and thetr particles with 
supporting texts on Electromagnetism and Special 
relattvity, and Atoms and thetr structure. 

The volume Basic quantum mechantes will provide 
the necessary introduction to many sections of this 
new text on The properties of nuclet, the level of 
which is pitched to take over where the corresponding 
sections of Atomte nuclet and thetr particles leaves 
off. 

Dr. Jones concentrates his attention on the 
available data and their interpretation, developing 
a natural sequence from the properties of nuclear 
forces, through nuclear models to nuclear decay and 
reactions including fission and fusion. The study 
of these topics is undertaken in second- or third- 
year courses in physics, applied physics, and nuclear 
engineering and challenges the student to appreciate 
both the power of the available experimental and 
theoretical methods and also the wide range of oppor- 
tunities for research in the nuclear many-body pro- 
blem. The sureness of approach found in this treat- 
ment illuminates the scene in a rewarding fashion 
which reflects the author's experience of teaching 
and developing the subject over a number of years. 


March 1977. E.J.B. 


Preface 


The level I have aimed at in this book is that of the 
nuclear-structure course given to Oxford students in 
the final term of their second year. In their third 
year some, but not all, pursue a more advanced course, 
according to their choice of options. The book there- 
fore reflects the requirements of 'mainstream' under- 
graduate physics courses. It presupposes a certain 
familiarity with particle accelerators and radiation 
detection. 

The content is chiefly devoted to the static and 
quasi-static properties of nuclei and nuclear excited 
states. But I firmly believe that all physics students 
should have some notion of how a nuclear reactor works, 
and so a chapter on nuclear reactions has been included 
in addition to a brief account of fission (and fusion). 

In a book of this length much has to be taken on 
trust; quantum-mechanical formulae have been used 
freely with no more than a passing reference (the par- 
ticular reference is to a text which I found useful as 
a student and still do as a teacher, but a more recent 
publication may appeal to the reader). Whatever the 
choice, I hope my book will encourage the student to 
acquire a working knowledge of basic quantum mechanics, 
which forms the framework of his (her) studies. 

I am indebted to members of the Department of 
Nuclear Physics, Oxford, who read and commented on the 
script, and especially to Dr. M.G. Bowler with whom I 
had many valuable discussions. I should also like to 
thank one of the editors of this series (E.J.B.), who, 
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I feel, put more time than his editorial duties deman- 
ded into a critical appraisal which I much appreciated. 
Finally, a number of the figures were taken from, 
or based on, published work. The figure captions in- 
dicate, where relevant, the author and the work. My 
thanks are also due to the publishers for permission 
to reproduce published work: They are: Fig. 1.2, 
Springer Verlag; Figs. 2.1, 3.8, 3.9, 4.3(a),(b), 
North Holland; Figs. 2.4, 6.3, American Institute of 
Physics; Fig. 3.2, Pergamon Press; Fig. 3.3, Annual 
Reviews; Fig. 3.4, Academy Munksgaard; Fig. 4.2, 
McGraw Hill; Fig. 5.2, John Wiley and Sons; Fig. 6.4, 
Brookhaven National Laboratory; Fig. 7.3, W-A. Ben- 


jamin. 
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1. General properties of nuclei 


INTRODUCTION 
From early measurements on Coulomb scattering it is 

known that the atomic nucleus is a small central par- 
ticle (radius <10721 


whole mass of the atom and which carries a positive 


cm) which possesses virtually the 


charge equal to the atomic number of the atom. From 
later work it is known that the nucleus is built up 

of protons and neutrons; the former supply the positive 
charge, whilst both contribute to the mass. From gene- 
ral quantum-mechanical considerations such a confined 
system will possess energy levels characterized by 
quantum numbers. Under normal circumstances, the nu- 
cleus will be in its ground state, and this chapter is 
concerned with nuclear properties in this state. To 
determine the properties of the nucleus in other (ex- 
cited) states it is necessary to set up such states 

by means of nuclear reactions. 

Thus a particular nucleus (in its ground state) 
must be described in terms of the quantum numbers 
"angular momentum' (also loosely called spin) and 
'parity' (to be explained later in the chapter). This 
is the minimum number of quantum numbers required by 
the symmetry properties of the system. There are 
others, but they are approximate or model-dependent, 
i.e. they depend on a simplified picture of the nucleus, 
and, unlike the corresponding atomic system, no one 
model stands out sufficiently to attach extra signifi- 
cance to its quantum numbers. 

In addition, a system of finite extent can be des- 
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cribed in terms of moments of the matter distribution. 
In the case of the nucleus these moments are usually 
determined from its interaction with the atomic elec- 
trons which is electromagnetic in origin, so the moments 
measured are electric and magnetic multipole moments. A 
system can have an electric dipole moment only if it has 
an up-down directional asymmetry in space, e.g. the 
molecule HCl, but not Cl,, can have a permanent electric 
dipole moment. For nuclei no such asymmetry is believed 
possible, so nuclei do not possess such a moment. Since 
the components of the nucleus, protons and neutrons, 
have intrinsic spins and magnetic dipole moments, the 
nucleus itself can also have a magnetic moment. The two 
important moments of the nucleus, apart from the elec- 
tric monopole, are therefore its magnetic dipole and 


electric quadrupole moments. 
Whilst nuclear physics deals with all nuclei, of 


especial interest are the so-called 'stable nuclei' 
These are the nuclei which exist in the world around 
us, and this is quite specific, though 'stable' itself 
is much less definite. A nucleus may be unstable 
against emission of a neutron or a proton, or a com- 
posite particle such as an a-particle, or against B- 
decay. Of these neutron emission is very fast 

(ce aor? 
but still very fast, whilst the latter two can have 


or less), proton emission generally less so 


lifetimes varying over a vast range depending on the 


energy available (tT v us to %v 102° 


years). By 'stable' 
a lifetime of the order of the age of the universe 

(v 1029 
stricted meaning 'stable against the @-decay process’. 
THE SIZE 


The original measurements made by Rutherford and 


years) is implied; sometimes it has the re- 
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his co-workers have now been augmented by more accurate 
work, which has not only measured an effective radius 
for the nucleus but has also obtained a measure of the 
rate of fall-off of density, or charge, at the peri- 
phery. The measurements fall into two classes: (1) 
those that determine the charge distribution and (2) 
those that determine the matter distribution. 


The charge dtstrtbutton 

These measurements are made with charged leptons 
as probes. Leptons form.a class of fundamental par-' 
ticles which interact only weakly with protons and 
neutrons, except for the long-range Coulomb term where 
relevant.’ Thus the interaction between an electron 
Or a muon and nuclear matter is completely dominated 
by the Coulomb interaction. The greatest source of 
information is the elastic scattering of electrons. 
Accurate measurements of optical and X-ray spectra 
and muonic X-ray spectra have also given information on 
the nuclear size. 


Electron elastte scattering 

Outside the nucleus of charge Ze the electron ex- 
periences a force « Z Ir 2, but inside the nucleus 
(assumed spherical) the force is « Z*(p)/r®, where 
Z'(r) is the charge contained in the sphere of radius 
yr. Thus small-angle scattering, which corresponds to 


TInteractions between fundamental particles can be 
put into the categories strong, electromagnetic, weak 
and gravitational, which are roughly in strength as 


10 : 1/137 : 10°’ : 10°” (see e.g. Perkins 1972). 
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a classical trajectory rather distant from the nucleus, 
will obey the Rutherford formula and therefore give no 
information about the distribution of charge, whilst 
large-angle scattering will show a significant devia- 
tion from the formula. The ratio to Rutherford scat- 
tering as a function of angle is determined by the form 
of Z'(r) and therefore by the charge distribution. 

From the uncertainty principle, if we wish to attach 
Significance to changes of density over a distance Ar, 
then the incident particle must be capable of trans- 
ferring momentum hq to the target nucleus such that 

qhr ~ 1; obviously the wavenumber of the incident beam 
must exceed q. For Ar % 10° +5cm, electrons of kinetic 
energy > 200 MeV are required. Some twenty years ago, 
measurements were made at 50 MeV and were capable of 
defining a suitably averaged charge radius only; more 
recently, the energy range has been extended to 600 - 
1000 MeV and the radial dependence of the nuclear 
charge distribution has been revealed in fine detail. 


X-ray and opttecal spectra 

Since the nucleus has finite extent its field of 
force on the atomic electrons will not be « 1/r? when 
the electron is within the nucleus. If therefore the 
atomic. energy levels can be calculated accurately for 
a point nucleus and compared with their measured values, 
the nuclear size can be deduced. But the calculation 
is an unsolved many-body problem. For X-ray transi- 
tions the screening effect of the outer electrons is 
not large and can be approximated to; but for the 
tightly-bound inner electrons there are uncertainties 
even in the solution of the two-body problem, whilst 
for the valence electrons, although these latter un- 
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certainties do not arise, the screening effects cannot 
be calculated to sufficient accuracy. In certain cases 
useful information can be derived. It is reasonable 

to assume that, for a comparison between spectral lines 
from isotopes of the same element, the differences will 
arise from the differences in nuclear mass and size 
only. For light elements the mass difference produces 
effects large enough to measure, but for heavier ele- 
ments the measured effects must be ascribed to dif- 
ferences in nuclear charge distribution. Thus the 
method of isotope shift gives information on the way 
the charge distribution changes on adding neutrons 


to a nucleus. 


Muonte X-rays 

Negative muons u , like electrons, can be captured 
into atoms and interact electromagnetically with nuclei. 
Since they have a mass %v 200 me (where m is the mass 
of the electron), their Bohr orbits will have 1/200th 
of the radius of the equivalent electron orbits. For 
a heavy element such as lead, the lowest Bohr orbit 
(since no orbits are occupied by other muons, the 
exclusion principle does not operate here) is such 
that the muon spends most of its time within the nucl- 
eus. The size effect is therefore large, reducing the 
lowest p > s transition from ~ 16 MeV for a point 
nucleus to ~% 6 MeV as measured. The muon is, of course, 
unstable and decays spontaneously to an electron and 
two neutrinos (v) with a mean life of % 2 us. In the 
lowest Bohr orbit in lead it decays in about half this 
time because of the competing mode p +p->n + Va 
but during its lifetime it traverses several metres of 


solid nuclear matter (density ~ 3 ~x 10° kg fs thus 


6 GENERAL PROPERTIES OF NUCLEI 


emphasizing the weakness of the weak interaction. 

From a theoretical aspect, the computational dif- 
ficulties are the same as those encountered in inter- 
preting electronic X-rays, but the uncertainties are 
reduced. Obviously the screening effects of the 
atomic electrons will only be of the same order as for 
the K-electrons, whereas the nuclear interaction ener- 
gy is increased by the factor 200. It also turns out 
that the uncertainties in the two-body problem are re- 
duced by the same factor. Since the bound muons have 
effectively a long wavelength (% size of the Bohr 
orbit) the measurements lead to an average nuclear 
radius and give no information on the charge distri- 
bution. But the average they give is very precise, 
and nowadays it is used as a datum in the analysis of 
electron scattering, acting as a constraint when look- 
ing for the finer details of the charge distribution. 


The matter dtstrtbutton 

If the nuclear probe interacts strongly with pro- 
tons and neutrons then the size determined from the 
interaction will be that of the nuclear matter as a 
whole. We have also the added complication that the 
strong interaction has a range which, though short, 
is not short compared to the nuclear size, and due 
allowance must be made for this interaction when de- 
ducing the mass distribution. 

The historical example of this kind of measure- 
ment is the elastic scattering of a-particles, which 
led to the first determinations of the nuclear size. 
Other charged particles have since heen used, e.g. 
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protons, carbon and oxygen nuclei (in the form of par- 
tially stripped ions), and other heavy ions. From a 
theoretical point of view, the simplest process is the 
scattering of neutrons, but experimentally this is 
rather difficult, since it involves three nuclear inter- 
actions - one to produce the neutron, another to scat- 
ter it, and a third to detect it. Other methods in- 
volving the electrostatic term are of historical im- 
portance: they are the lifetime for o-radioactivity 
and the determination of the electrostatic energy of 
the nucleus from the energy difference between mirror 
pairs. Recently measurements of K-mesic X-rays have 
proved to be a promising method of probing the nuclear 
surface.' 


Scattering of charged particles 

Experimentally, the scattering cross-section (see 
Appendix A) is measured as a function of angle for a 
fixed beam energy or as a function of energy at a fixed 
angle. In either case the ratio of measured cross- 
section to Rutherford cross-section is plotted as a 
function of the variable. The ratio will be unity over 
a range of angles from near zero up to a certain angle, 
beyond which it drops off, indicating absorption of 
the incident particle. Alternatively, for constant 
angle the ratio will be unity up to a certain energy, 
beyond which it drops off. However, for protons or 
even a-particles on light nuclei it is possible for 


the cross-section to rise above the Rutherford value, 


Tk-mesons, which can be positively or negatively 
charged or neutral, have a mass v 1000 m, and inter- 


act strongly with nucleons. 


8 GENERAL PROPERTIES OF NUCLEI 


because such particles stand a good chance of being 
emitted after absorption. For heavier ions, or protons 
and a-particles on heavy nuclei, emergence of the ab- 
sorbed particle is improbable. 

The region of fall-off is correlated with the 
classical trajectory passing through the periphery of. 
the nucleus, immediately giving an idea of the size 
of the nucleus. How to obtain a more precise deter- 
mination of the size was uncertain before the develop- 
ment of the optical model (see Chapter 3). We should 
notice that, in the interpretation of the results, in 
addition to making allowance for the range of nuclear 
forces it is also necessary to take into account the 


finite size of the incident particle. 


Neutron scattertng 

When fairly energetic neutrons (% 10 - 30 MeV) 
impinge on nuclei they are strongly absorbed from the 
initial beam, though they may subsequently be emitted 
from the resulting excited nuclei. The degree to which 
they are emitted depends upon the number of competing 
modes of decay (see Chapter 6). If there are many 
such modes, and there are for E ~ 10 - 30 MeV, elastic 
scattering will be reduced to a minimum. Elastic scat- 
tering cannot be zero since an extension of Babinet's 
theorem in optics to this nuclear wave problem in- 
dicates that absorption without scattering is impossi- 
ble. This scattering is predominantly through small 
angles (up to 1/kR as in optical diffraction). The ab- 
sorption cross-section can be conveniently measured by 
a beam-attenuation technique. A neutron detector 
placed in the path of a collimated neutron beam deter- 
mines the beam flux before and after interposition of 


GENERAL PROPERTIES OF NUCLEI 9 


the absorbing target. If the angle subtended by the 
detector at the target is greater than the angle of 
diffraction, then diffracted neutrons are detected, 

and the beam attenuation is a true measure of nuclear 
cine St Re A 
where the additional term X, the wavelength of the pro- 


absorption. <A simple theory gives o 


jectile, expresses the uncertainty with which we can 
locate a neutron. A more sophisticated analysis can 
be carried out in terms of the optical model. 


a-Radtoacttvity 

From consideration of the mass equation (See pp. 
12-20), all nuclei above A ~ 150 are unstable, in 
principle, and can break up into two smaller units. Be- 
cause of its high intrinsic binding, the a-particle 
is usually one of these units. The process involved 
will be discussed in more detail in Chapter 4, but for 
the present it can be summarized by stating that the 
probability of a-particle emission depends on two fac- 
tors: the intrinsic nuclear factor and the barrier 
penetrability. The latter is a rapidly varying func- 
tion of the nuclear radius, which can therefore be 
determined if the former is known - the measured life- 
time is known with great precision. Unfortunately, 
the intrinsic factor is a notional one and cannot be 
measured; theoretical estimates are uncertain, with 
the result that the nuclear radius deduced is quali- 
tative only. It is probably more meaningful to use 
the data, and a nuclear radius from other data, to 


deduce the intrinsic nuclear factor. 


Electrostatic energy differences of mtrror nuclet 


Mirror nuclei are nuclei in which the proton num- 
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ber and neutron number are interchanged, e.g. AG and 
DON havine OS 6: iv ST and B S57 WS Gu MOSt- oF tne 
measurements refer to mirror pairs having |w - Z| = 1. 
It is believed (Chapter 2, p.42) that the speci- 
fically nuclear interactions will be the same for such 
pairs, so their mass difference will arise from the 
electrostatic interaction and the n-p mass difference. 
“The former term is discussed in greater detail in the 
next section, for the moment-a uniform distribution 
of charge in the nucleus is assumed to give for the 


mass difference 


1 (22+1)e° 


ATE R 2 


| 3 

AM = M(A,Z+1) - M(A,Z) = my Tt Pe 
where the M and m are defined in the next section. 

The mass difference can be measured accurately 


from the g” -decay spectrum of (A, Z+1), since 


where Plan is the maximum kinetic energy available to 
the Bg"; or from the threshold energy for neutron pro- 


duction by protons on AyZ) 


M(A,Z+1) + m. + Q 


M(A,Z) + my + m 


AM =m.+m.-m. - Q 


Thus the @-value of the (p,n) process is the negative 
of the electrostatic energy difference. 

From the data, R = Rae 
but the value deduced for Ro» namely 1°46 fm is rather 
large compared with deductions from other processes 


is found to hold good, 
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(see the end of this section). This reflects the 
dubious supposition that the nuclear charge is uni- 
formly distributed within the nucleus. Since the data 
on mirror nuclei refer only to light nuclei (A < 40) 
surface effects will be expected to be large. 
K-meste X-rays 

From the discussion of muonic X-rays, it follows 
that the K -mesic Bohr orbits will have an even smal- 
ler radius. In addition, the K -meson will interact 
strongly with the nucleus, so nuclear absorption will 
compete with X-ray emission. Absorption will occur 
from highly excited Bohr orbits when the K-meson over- 
laps the tail of the matter distribution. These meas- 
urements therefore focus attention on the details of 
the matter distribution at the surface of the nucleus. 


Conelustons on stze 

Analyses of the earlier data, which were capable 
only of defining an effective nuclear radius, gave: 
the result that the radius is proportional to the cube 
root of the mass number A, i.e. that the nuclear dens- 
ity is a constant independent of the nuclear size. 
Later measurements indicate that, although this is 
true for the central region, near the periphery the 
density falls off smoothly and quickly to zero. The 
precise nature of the fall-off is not determined from 
the data, but it has proved convenient to represent 
it by the (Fermi) function 


(o)r = [1 + exp (24) 71, 


where R, the radius at half-density, is referred to as 
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the nuclear radius and ais a surface thickness para- 
meter which defines the region over which the density 
is falling off. This latter parameter appears to be 
independent of nuclear size, whilst R « alls 

The measurements referring to the charge distri- 
bution are of greater precision than those referring 
to the matter distribution, but both agree on the 
values of RF and a; there appears to be no significant 
difference between the distributions of neutrons and 
protons within the nucleus. The data may be summed 


1/3 where 


up by giving the charge radius as RF = RoA 
Kq 


not known accurately. The matter distribution will 


= 1-1 fm and a‘ 0°5 fm. The latter figure is 


have the same parameters, but both are less accurate. 
An objection may be raised, by the more dis- 
cerning student after reading on, that nuclei have 
been treated as spherical whereas some of the more 
massive nuclei have quite large quadrupole moments. 
The reason is that the nuclei have not been aligned 
in their targets during the experiments described, 
and so the results give a time- and space-averaged pic- 
ture of the nucleus which will have spherical sym- 
metry even when it is intrinsically non-spherical. The 
averaging process may well be expected to produce a 
wider region of fall-off.i.e. a greater value for a. 


There appears to be some indication of this. 


THE MASS AND STABILITY 

Historically, the masses of neutral atoms were 
determined relative to (atomic) 165 = 16 atomic mass 
units and later relative to (atomic) 120 = 12 U (uni- 
fied atomic mass unit), by chemists using chemical 


reactions. The measurements were later improved upon 
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using physical techniques involving the trajectories 
of ionized atoms through electric and magnetic fields. 
With this information available the nuclear physicist 
found it expedient to use atomic masses in his equa- 
tions for nuclear reactions, since conservation of 
charge ensures that the same number of electrons will 
appear on each side of an equation (though not for 

B* -decay, see p. 21). Balancing the electron number, 
however, does not allow completely for the difference 
between atomic and nuclear masses, since no account 
has been taken of the binding energy of the atomic 
electrons - which, though small, can add up to 

~ 0-5 MeV for uranium. If this appeared directly as 
an error in the determination of @-values in nuclear 
reactions then it would be serious. But most reactions 
are concerned only with small changes of A and Z and 
therefore only with small changes in this term, which, 
if neglected, could introduce errors of order 20 keV 
for a-decay in the region of uranium when using nuclear 
masses - though the correction is well known and ac- 
curate to a fraction of 1 keV. On the other hand, use 
of atomic masses largely eliminates such effects since 
the mass measurements are done on ions which are ion- 
ized in the valency shell and thus already include all 
the large contributions to electronic binding. For 
this reason the atomic scale is to be preferred. How- 
ever, the student may come across the usage of nuclear 
mass (Oxford University Examination Papers are a case 
in point), and so, although the atomic scale will be 
employed here, comparisons will be made between the 
two scales. The nomenclature employed will be that 

M will refer to atomic masses and m will refer to 
nuclear masses. Since the masses quoted for the pro- 
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ton and neutron are usually m_ and m rather than 


’ 
Ms (mM, es anyway), the mass seuatione will contain 
2 (m,, + m,) rather than 2M), 

From the equivalence of mass and energy, the mass 
of a nucleus consisting of Z protons and W neutrons 
(Z + N = A) will be less than the sum of the masses 
of the individual nucleons by the energy of binding 
(expressed in u (1:66 x 10°77 kg) - it is in fact more 
convenient to re-express the masses in MeV and since, 
in u, the mass values are close to the integers A, it 
is usual to give the mass excess in MeV, i.e. 
{M(A,Z) - A} u converted into MeV). Thus 


M(A,Z) = Z(m, tM) A (A-2)m, - B. 


It is necessary therefore to consider all con- 
tributions to the binding energy. This is done in an 
empirical way, but guided by physical facts and theor- 
ies, in the mass equation of von Weizsdcker. 

From the fact that nuclear density (neglecting 
surface effects) is a constant a nucleus can be com- 
pared with a liquid drop, for which it is known that 
the binding energy (again neglecting surface effects 
which are usually small for liquid drops) is propor- 
tional to its volume, since all molecules deep in the 
liquid are in the same environment. Hence, by analogy, 
the leading term in the mass equation representing the 
binding energy will be -aA, where a is some undefined 
constant. All constants introduced into the mass 
equation will be defined to be positive. 

A second term in the binding energy also arises 
from the analogy with a liquid drop. The forces which 
keep the liquid drop together (van der Waals forces) 


GENERAL PROPERTIES OF NUCLEI 15 


are short range, so that the volume binding energy 
term is dominated by the interactions between near 
neighbours. At the surface a molecule will have fewer 
near neighbours than a molecule in the interior and 

so will be less tightly bound, giving a surface-energy 
term and a surface tension, the former being propor- 
tional to the surface area. It results in a lower 
total binding of the drop and therefore appears as a 
positive term in the mass of the drop. Since the sur- 
face area of a sphere «= oe and, reverting to the 
nucleus, r = aif 
as +pa2/5 


the surface-energy term will appear 
in the nuclear mass equation. As a large 
fraction (% 3) of nucleons are near the surface this 
will be expected to be a considerable fraction of the 
volume term. 

The analogy with a liquid drop neglects the charge 
on the nucleus. On attempting to remedy this by char- 
ging the drop, the analogy breaks down, since the 
charge migrates to the surface of the drop whereas, 
from the previous section, it is known that the charge 
density of a nucleus is proportional to the matter 
density. This illustrates the difficulties in making 
classical models of quantum-mechanical systems; ob- 
viously the uncertainty principle operates to reduce 
the energy when both protons and neutrons make full 
use of the space available. Assuming that the nuclear 
charge is uniformly distributed over the volume of a 
sphere then its electrostatic potential energy can 
easily be shown to be . (Ze) */ane 5k. If the discrete 
nature of the charge on the proton is taken into ac- 

2 should be replaced by Z(Z-1). This 
refinement is probably an improvement but, except for 


count, perhaps Z 


small nuclei, it makes little difference in the mass 
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equation. In any case there are larger factors of ig- 
norance: the charge distribution is known to vary 
with radius within the nucleus, and there could be a 
correlated motion which allows protons to avoid each 
other yet retains the required charge distribution. 
The electrostatic term is therefore introduced into 
the mass equation as +eZ7 alé3 where the (unknown) 
constant € is used to account for these uncertainties. 
At this stage the mass equation is 


M(A,Z) = Z(m, + mg) + (A-Z)m, - a4 + BA°/> + e274” 


1/3 
Minimizing this function with respect to Z at constant 
A gives, for stable nuclei, the approximate result 

x. gifs ? 
Z=A (nm, my 
which is far too slow a variation with A 
(see Fig.1.1). A = 32 has Z = 16 stable, so the above 


- m)/2e. The important factor is. 


Fig.1.1. Plot of 4 against Z of most stable nucleus for given A. 
For convenience only every sixth A-value is plotted. 
The full line shows W = Z. 
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formula would give Z = 32 stable for A = 256, to be 
compared with the naturally stable 23eu having A = 238, 
Z= 92, 

It is interesting to note from Fig.1.1 that, up 
to. A = 40, the line of stability has W = 2. Above 
this value W increases more rapidly than Z, presumably 
because of the electrostatic term in the mass equation. 
It would appear that a term « (w-2)*, i.e. (A-22)? is 
also required, which if of positive sign would mitigate 
in favour of WN. = Z for stable nuclei. One need not 
look far for the origin of such a term; the Pauli 
principle will ensure that the nucleus (in the absence 
of the electrostatic term) should have W v Z in order 
to utilize the lowest levels available to both species. 
To proceed further it is necessary to invoke a model. 
A convenient model, which incorporates the Pauli prin- 
ciple, is to consider the neutrons and protons as two 
non-interacting degenerate Fermi gases inside the 
nucleus. Neglecting the electrostatic term and the 
n-p. mass difference, which are already in the mass 
equation, the model instantly leads to NW = Z for mini- 
mum energy. Expanding the energy function from the 
point W = Z gives a leading term « (w-2)7/A. Hence 
a symmetry term +y (A-22)*/A is added to the mass 
equation. 

An extra term is needed before attempting to use 
the mass equation to determine which nuclei are stable. 
In Appendix B it is shown, under approximation, that 
the balance between terms so far developed in pre- 
dicting the heavy stable nuclei is so finely poised 
that the nearest nucleus (of given A) to the line of 
stability will be most stable. Thus, starting from an 
even-Z even-N stable nucleus, such that the next 
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stable mass is obtained by adding a neutron, then the 
second nucleon added to form a stable nucleus will be 
a proton more often than a neutron since, if (A ,Z) is 
on the line of stability then (A + 2, Z + 1) will be 
nearer that line than (A + 2,2). Thus, for any given 
A, only one stable nucleus will be expected, and of 
those with even A, stable nuclei will occur as often 
with odd Z as with even Z. 

Nuclear charts show that there are 166 even-Z 
even-N stable nuclei but only 4 of odd Z, odd WN, and 
fH, ane ‘OB, ISN. Conditions 
for very light nuclei do not conform to those under 


these are very light: 


which the mass equation was developed, and in addition 
the systems have few widely-spaced energy levels which 
could affect the balance between Z and W. Of the 105 
odd-A stable nuclei, 50 are found to have odd Z. 

In order to account for this distribution of 
stable nuclei, a pairing energy is postulated and is 
accounted for in the mass equation by a small term 
6(A,Z), where 6 is a function which has the value zero 
for all odd A and takes a negative value when both A 
and Z are even but a positive value when A is even and 
Z is odd. Its magnitude drops off as A increases, and 
this is to be expected since, on average, the two similar 
pairing nucleons will be further apart for large A than 
for small A. 

The atomic mass equation therefore reads 


M(Z,A) = 2(m, +m.) + (A-Z)m, - ad + par/> 4 
P (1.1) 
aaa ot, Wye , 


This five-parameter formula has been fitted to the 
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available mass data (several hundred nuclei neglecting 
the very light ones, A < 20) to give values: 


Parameter Ol B Y € 6 


alu 


Value (MeV) 15-6 17-2- 23°3 0-7 331 A 


The fit is shown in Fig.1.2, where it was more conveni- 
ent to present the binding energy per nucleon (B/A) 
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Fig.1.2. Plot showing the binding energy per nucleon as a function 
of A. Locations of magic numbers are indicated. The 
inset continues the plot and illustrates o-particle 
structure in light nuclei. (Based on Wapstra (1958). 
Atomic masses of nuclides, Handb.Phys. XXXVIII/I. 


Trhe values of the constants vary from source to 
source. Usually no errors are given so the signifi- 
cance of the variations cannot be assessed. The fits 
appear to predict the masses of stable nuclei, and 
their neighbours on the mass diagram, to about 1 per 
cent in binding energy. 
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against A. To give an idea of the relative importance 
of the terms, for the stable nucleus (125,52) the con- 
tributions to B/A are + 15°6 MeV (volume), - 3°44 MeV 
(surface), - 0°65 MeV (symmetry), and -3-0 MeV 
(charge). The variation with A is shown in Fig.1.3. 
Notice that, on the scale shown, the pairing energy 
term is quite negligible; for the neighbouring even- 


Volume energy 


Surface energy 


Asymmetry 


Coulomb energy 


Contributions to B/A (MeV per nucleon) 


50 100 150 200 


Fig.1.3. Contributions to B/A of the individual terms in the 


mass formula. 


even nucleus to 4=125, it contributes an additional 
8 keV to B/A. 


Stability agatnst B-decay 
8 -decay is represented by the reaction 


(4,2) + (A,2+1) +e + vy, 
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where Vv, and ye represent electron neutrino and elec- 
tion antineutrino (see Chapter 4). This reaction is 


energetically possible if 
m(A,Z) > m(A,Z+1) + m, (1.2) 
since m., 0. But 


m(A,Z) = M(A,Z) -2Z me 


leading to 
M(A,Z) > M(A,Z+1) . (1.2a) 
B’-decay is represented by the reaction 


AG (491) ee Ve 


and is energetically possible if 

m(A,Z) > m(A,Z-1) + M (1.3) 
leading to 

M(A,Z) > M(A,Z-1) + 2m, (1. 3a) 


For these two types of decay it seems that the 
nuclear masses produce a more attractive statement of 
condition. However, there is another type of decay 
- an alternative, competing, mode to 8” -emission - 
namely, electron capture. Nuclei are not required to 
be stable in isolation, but in their natural habitat 


when surrounded by atomic electrons. Any one of these 
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electrons can be captured, though if energetically 
possible K-capture is likely to be dominant, leaving 
the atom strongly excited. However, in discussing the 
limit of stability the possibility of capture of the 
outermost electron, leaving the atom in its lowest 
state, must be considered, since energetically it is 
most favoured. 

This reaction is represented by 


(A,Z) +e + (A,Z-1) + v 


e 
and is energetically possible if 
m(A,Z) + m, > m(A,Z-1), (1.4) 
leading to 
M(A,Z) > M(A,Z-1). (1.4a) 


These equations indicate that electron capture 
can proceed when 8*-emission is impossible and also 
that the atomic mass scale is again favoured: the 
condition for stability against the 8-decay process is 
that the atom (A,Z) should have least mass. 

For the heavier nuclei this can be expressed, 
approximately, as (oM/0Z) , = 0, which is more con- 
venient to use than the two inequalities. The corres- 
ponding condition on the nuclear mass scale 
(am/3Z) , SOS is probably not so convenient to remem- 
ber nor indeed is the differential condition on the 
binding energy, (9B/92) , St FS CM 

From these equations, two nuclei (A,Z) and 
(A,Z-1) can both be stable only if they have exactly 
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the same atomic mass. This statement is subject to 
approximations concerning the neutrino mass and the 
electron binding energies. Taking these into account 
might produce a band of stability a few eV wide. The 
chance that such a pair of nuclei should have equal 
atomic masses to within a few eV is remote, so it is 
concluded that, in the few such examples naturally 
occurring, one of the pair is unstable but has a very 
long lifetime, e.g. 87 Rb and 87cr - (87,37) and (87,38) 
- of which the former has a lifetime of 5 x 10°" years, 


emitting 8 -particles with maximum energy % 274 keV. 


Effect of the patring term on stabtltty 

When A is constant and odd the mass values for 
different Z will define discrete points lying on a 
parabola. For A constant but even, it is necessary 
to define two parabolae displaced from each other by 
26(A), which is roughly 2-5 MeV at A = 81, falling to 
~ 1 MeV at A = 256. The odd-Z mass points will fall 
on the upper parabola and the even-Z mass points on the 
lower. Fig.1.4 gives examples of the parabolae for 
odd A and even 4A. 

The displacement of the two parabolae for even A 
is sufficient to ensure, even when the lowest odd-Z 
mass point is at the minimum of its parabola, that 
both its neighbours (which will lie on the lower para- 
bola) have lower mass. This can be checked for the 
two masses given above by evaluating the second deri- 
vative of the mass formula without the 6 term; the 
mass difference between this odd nucleus, assumed to 
be at the minimum of the parabola, and its two even 
neighbours will then by 3(9¢m/az°) ,. At A = 81 this 
is ~ 1°5 MeV and at A = 256, ~ 0-6 MeV, so that the 
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—90 


Fig.1.4. Mass parabolae for odd A and even A. The mass excess 


(in MeV) is plotted against Z for A = 101 and 102. 


pairing term is amply large enough for its purpose. 

These considerations lead to the conclusions that, 
for odd A, only one nuclide is expected to be stable 
against 8-decay, but for even A two apparently stable 
nuclides is likely to be a common occurrence (and is) 
and even three is a remote possibility (one or two 
examples are found). The term ‘apparently stable' has 
been used because one such stable even-even nucleus 
should be able to decay to its less massive even-even 
neighbour by a double §-decay process, which as yet 
has not been observed. 


Stabtltty agatnst nucleon or a-parttcle emtsston 


From Fig.1.2 it is seen that for A > 20 the bind- 
ing per particle for nuclei near the minima of the mass 
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parabolae (i.e. along the 'valley of stability' on 
a three-dimensional M,A,Z plot) lies between 6 MeV and 
8 MeV. This large value ensures that even the most 
massive naturally occurring nuclei (A ~ 240) are stable 
against emission of proton or neutron. 

The a-particle, on the other hand, is itself 
bound by ~ 28 MeV, so that instability to a-emission 
is much more probable. The condition for instability 
is 


M(A,Z) > M(A-4, 2-2) + M(4,2) 
or 

B(A,Z) < B(A-4, 2-2) + B(4,2) 
However, Fig.1.2 presents B/A as a function of A; if 
we denote this by $(A) - along the line of stability 


Z is a function of A, and so B(A,Z) is a function of 
A only - then 


A o(A) < (4-4) 9(4-4) + 4 9(4) 
But 
o(A-4) ~ (4) - 4 do/da 
giving, after dividing by 4, 
o(4) < $(4) - (4-4) do/da 


From the graph d¢/dA v 0-75 MeV/100 and is roughly 
constant well above the maximum of the curve; and 


26 GENERAL PROPERTIES OF NUCLEI 


o(4) ~ 7 MeV. 

Thus the condition for a-instability is that 9(A) is 
less than ~ 7-8 MeV at A = 100, ~ 8-1 MeV at A = 150, 
and ~ 8-5 MeV at A = 200. a-decay is therefore pos- 
sible for A > 150 (though shell effects have an in- 
fluence; see Chapter 3). Because of the mechanism of 
the process (see Chapter 4) it becomes a significant 
mode of decay only above A %v 210. 


ANGULAR MOMENTUM (SPIN) AND STATISTICS 

Quantum mechanically, as classically, an isolated 
system has a well-defined angular momentum. It is 
therefore a characteristic property of the nuclear 
ground state - and indeed any other well-defined state. 
As such it is loosely referred to as spin of the nu- 
cleus, even though it is a combination of the intrinsic 
spins and orbital motions of the individual nucleons. 
The reason for this is probably historical. Nuclear 
spins of stable and long-lived nuclei were determined 
from their magnetic interactions with the atomic elec- 
trons or with externally applied fields. For such 
interactions the spin vector is an intrinsic property 
of the nucleus just as it is, at present, for an ele- 
mentary particle. 

Since the intrinsic spin of a nucleon is 3}, whilst 
all orbital angular momentum vectors are characterized 
by integer quantum numbers, it follows that the sum of 
A such spins and orbital terms will lead to an angular- 
momentum (or spin) quantum number J which will be an 
integer when A is even but (integer + 3) when A is odd. 
Examination of a table of nuclear spins confirms this 
and also reveals a striking fact: all even-Z even-N 


nuclei have ground-state spin 0, without exception. 
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Since a pairing energy term, favouring the pairing of 
like nucleons, has already been introduced into the 
mass equation, it seems reasonable to suggest that 
such pairs individually couple to zero angular momen- 
tum. This point will be taken up in Chapter 3. 


The Pauli principle requires that quantum-mechani- 
cal systems have wavefunctions which are antisymmetric 
to the exchange of like fundamental particles. If 
the system under consideration is atomic (for definite- 
ness, consider a molecule made up of two like atoms) 
then the electronic wavefunction should have a definite 
symmetry on interchange of two nuclei as a whole, 
whilst the nuclear wavefunction should be antisymmetric 
on interchange of two protons (or neutrons), either 
within a nucleus or between two nuclei. If the like 
nuclei have mass number A then their interchange is 
equivalent to A interchanges of nucleons. Thus the 
electronic wavefunction will be antisymmetric for 
odd A and symmetric for even A. The symmetry reflects 
in the statistical behaviour of an aggregate of such 
atomic systems. Antisymmetric states conform to the 
Fermi-Dirac distribution, whilst symmetric states 
conform to the Bose-Einstein distribution (Martin 
1975; Pauling and Wilson 1935); for a diatomic 
molecule of similar atoms, the intensities of lines in 
a rotational band spectrum reflect the statistics 
obeyed and so determine wnether the number of nuclear 
particle components is odd or even (see Burcham 1963). 
This may sound trivial now, but before the discovery 
of the neutron it led to rejection of the hypothesis 
that nuclei were composed of protons and electrons. 


In particular, for 14, (14,7) the N, molecule obeys 
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Bose-Einstein statistics, although on that theory 14 


was thought to consist of 14 protons and 7 electrons - 
an odd total of fermions. 

The symmetry rules impinge on nuclear reactions 
in a somewhat different way. Consider the system 
(ata); its wavefunction must be symmetric under the 
interchange of the two a-particles. If these particles 
can come together to form a state in 8Be then this 
state must reflect that symmetry. But BR e is more 
complex than a (2a) structure (e.g. it can also be 
formed from ig + p) and each state can be looked upon 
as a conglomerate of such structures. However, of all 
the states in Sie which exist only those with a certain 
symmetry can decay into ata. These states must have 
J= 0, 2, 4, etc., as can be shown by a consideration 


of the wavefunctions for ata. 


PARITY 

Parity is another aspect of symmetry in quantum 
mechanics, but one which has no direct correspondence 
in the classical limit. If the interaction potential 
of a system has the symmetry property 
V(wzy,2) = V(-x,-y,-2) = PV(x,y,2), where the operator 
P has been introduced and defined, and the Schrodinger 
wave equation has a solution y(z,y,2), then Pyi(z,y,2) 
is also a solution. For a non-degenerate case, the 
second solution must be a numerical factor times the 
first, i.e. PW(xz,y,2) = ab(xz,y,2). But it is obvious 
that pe = 1, corresponding to two inversions, so that 
a = + 1. Thus a non-degenerate solution corresponding 
to ths type of potential has either even parity 
(P = + 1) or odd parity (P = - 1). 

For a single particle in a central potential there 
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1S a one-to-one correspondence between the orbital 
angular momentum and the parity, which follows from the 
angular solution of the problem: 


even ¢ = even parity , odd 2 = odd parity. 


For a multi-particle assembly like the nucleus, short 
of a complete solution to the problem, the parity can 
only be deduced in terms of a model in which the 
parities of the individual particles can be specified, 
as in the shell model (Chapter 3). From the fact that 
the complete wavefunction can be expanded in terms 
which are themselves products of the wavefunctions of 
the individual particles, it follows that the overall 
parity is a product of the individual parities if + is 
assigned to even parity and - to odd parity. 

As a constant of the motion, parity plays an im- 
portant part in nuclear reactions. If in the reaction 
X + a > W* + Y + b, where W* represents a state in the 
compound nucleus (Chapter 6), the spins and parities 
of all five particles are known, then conservation of 
angular momentum and parity severely restrict the 
orbital angular momentum taken in by a and taken out 
by b, thus affecting both the angular distribution 
and the magnitude of the cross-section. As an example 
assume in the above that X has (spin, parity) = (07), 
a is a proton (47), b is a neutron eae Y is e. 
the compound nucleus W* is a Then, conserving angu- 

= lor 2, Ly = 0,1,2,3, or 
4, but parity conservation will reduce these possibi- 


lar momentum only gives La 


lities to La = 1; Ly = 1 or 3. For spinless particle 


and projectile, the accessible states in the compound 


12 16 


nucleus are severely restricted, e.g. for C+ a O* 
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then the excited states must be 0°, 1, 2°, 3, etc. 


ELECTRIC AND MAGNETIC MOMENTS 


Electrtc moments 

If a charge is distributed over a finite region 
of space, surrounding a chosen origin, through which 
an electric field derivable from a potential U(xz,y,2) 
has been imposed from outside, then the interaction 
energy can be expressed in a Taylor series thus, 


E = [u(e.y,2)o(e,y,2) a = u(0) | pdt + 


+ higher-order terms. 


The first term gives the interaction energy with 
the field of a point charge at the origin equal to 
the total charge - the monopole term. The next three 
terms represent the scalar product of a vector dipole 
- the dipole moment of the charge distribution - with 
the vector field at the origin; whilst the next six 
terms represent the scalar product of two tensors, 
one the generalized derivative of the vector field and 
the other the quadrupole moment of the charge distribu- 
tion. These classical definitions hold for a quantun- 
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mechanical system if we interpret p as e D vE(r)b(), 
where the summation ¢ is taken over the charged par- 
ticles of the system (protons in the nucleus). 

The fact that a nucleus has a well-defined parity 
makes p(r) = p(-r). Thus every odd moment - dipole, 
octupole, etc. - will vanish identically. The simplest 
deformation from spherical results in an ellipsoid of 
revolution, which, if the z-axis is the symmetry axis, 
produces a quadrupole-moment tensor which can be speci- 
fied by a pe parameter @) given by 
eay = [ (32? - 7” 2) dt - a definition which conveniently 
has 29 = 0 for spherical symmetry. A positive value 
for Qo corresponds to a prolate shape for the nucleus, 
i.e. elongated along the symmetry axis; whilst negative 
Qo defines an oblate nucleus, i.e. flattened along the 
symmetry axis. As defined the dimensions of @ are 
those of an area for which the convenient nuclear unit 
is the barn (b)(1 barn = 10°*4em* = 107282), Carrying 
out the ee — p = Peers ae within the ellipsoid 

R 


R- 
major and minor semi-axes, R the average nuclear radius, 


and AR = a-b. 
Obviously a state J = 0 has spherical symmetry 


gives Q@) = Sia -b ~< ZR » where a, b are the 


and so has zero observable quadrupole moment; less 
obviously, a state J = 3 also has spherical spatial 
symmetry and zero observable quadrupole moment. Table 
1.1 illustrates these points and also indicates that 
some nuclei can be strongly deformed with AR/R ~ 30 


per cent. 


Magnette dtpole moment 


A similar, but more complicated, expansion of the 
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32 
Nucleus B/MeV 
+ 
. ; 
tH : 
y 2-225 
jHe 28-29 
SLi 39-24 
Mo 127°6 
F3Cl 2982 
S7Fe 499-9 
tSb 1026 
*70Lu 1418 
neta 52 
7A Bi 1640 
253i) 1785 


TABLE 1.1 
Atomic 

binding s 
keV - 
17 
0-014 1° 
0-014 1” 
0-079 O° 
o-190 3" 
20 o* 
12 > 
34 1 
180 3* 
370 . 
390 5 
550 > 
690 I 


GROUND STATE 


+0° 


+3° 


+0° 


+Q- 


+3> 


+3° 


+2° 


+4> 


u/nm 


2564 


0905 


342 


180 


35 


080 


°35 


Q/d 


+0°00282 


-0°04 


-0°08 


#0 #35 


AR/R 
+t 


+0°35 


-0°02 


+0°25 
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Table 1.1 continued 


Footnotes: + B is given only to four figures; it is 
however known to the nearest keV for the light 
elements and somewhat less precisely for the 
heavy. 

tt The AR/R have been computed using 
R= penate fm. A relationship 
Q/4¢ = J(27-1)/(J+1)(27+3) is used to connect 
the observed quadrupole moment @ with the in- 
trinsic quadrupole moment @p. 


charge distribution and its interaction with an ex- 
ternal magnetic field, taking into account a possible 
circulation of the charge within the distribution will 
give a leading term dipole in nature. If in transition 
to the particle system the intrinsic magnetic moments 
of the individual particles are taken into account, 
this dipole term will have two parts: (1) from the 
circulating currents of the proton distribution which 
will be proportional to some suitable summation of the 
orbital angular momenta Ls and (2) from a suitable 
summation of the intrinsic magnetic moments of both 
protons and neutrons which will be vectors parallel 

to the spins Sy. The phrase 'some suitable summation' 
is used because the resultant magnetic moment of a num- 
ber of particles depends on how the particles are 
coupled together; recourse to a model is necessary, 

as will be apparent in Chapter 3. At this stage it 
should be noted that the magnetic dipole moment is a 
vector, and it is a fundamental result in quantum 
mechanics that all vector constants of a system are 
scale factors of one vector, which, for an isolated 
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system, will be the total angular momentum J. Thus 

vectorially u = gJ, where g is the suitable scale 

factor. It is convenient to represent the dipole 
moment by a scalar, which in the classical limit would 
represent the magnitude |u|. Two possibilities exist: 
gig (J+1) # and g(™) Wax = gJ. The latter has been 
adopted. 

For the case of a spinless particle of mass m and 
charge e, J = L and g is the gyromagnetic ratio e/2m - 
a relationship which can readily be verified for a cir- 
cular orbit, using Ampere's theorem connecting a cir- 
culating current with a magnetic shell. This relation- 
ship leads to a natural unit eh/2m for the dipole mom- 
ent. In atoms the mass is that of the electron, de- 
fining the Bohr magneton (up); in nuclei it is the 
mass of the proton defining the nuclear magneton (uy) > 
which has the value 5-051 x 10°’ 3 T+. The intrinsic 
magnetic dipole moments of the proton and neutron are 
respectively +2-79275 Uy and -1-9135 Uy: A few nuclear 
magnetic moments are given in Table 1.1. 

PROBLEMS" 

1.1. Assuming that the nuclear size parameters are as 
on p. 12, calculate the quantum number 7 for which 
the circular Bohr orbit of K_ around *05Pb lies 
in nuclear matter at a density of 10 per cent of 
the central nuclear density. The mass of the K- 
is 494 Mev/c?. 


1.2. By assuming that the neutrons and protons in the 


Ts indicates that the problem is taken from, or 


based on, an Oxford Finals question. 


4% 


oe 
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nucleus behave as independent Fermi gases, and 
by neglecting the Coulomb interaction, derive 
the form of the symmetry term, « (w-2)*, which 
appears in the mass formula. 


. Verify the statement that the pairing term in the 


mass equation is sufficiently large to ensure’ 
no odd-odd nucleus (above A=40) is less massive 
than both its even-even neighbours of the same 


mass number. 


-Calculate values (in MeV/c") for y and € in the 


mass equation from the following facts: er emits 


positrons with a maximum energy of 4-95 MeV; 


ree Ba is the stable isobar of mass number 135. 


-The maximum energy of the positrons from the 


decay of TON is 1+24 MeV and there is no sub- 
sequent y-radiation. Deduce a value for the 
radius of nuclei of mass 13. (The n-p mass dif- 
ference is 1:29 MeV/e*.) 


Show that the angular wavefunction siné exp(id), 
corresponding to 7Z=1, m4=1, has odd parity, 
whilst sin’ @ exp(i2¢) and sin® cos@ exp(id), 
corresponding to 2=2, m,=2 and 1, have even 
parity. 


Deduce the gyromagnetic ratio for the case of a 
charged particle moving in a circular orbit. 
Hence derive the natural units for the magnetic 


moments of electron and proton orbits. 


2 : Nuclear forces 


INTRODUCTION 
From the last chapter a picture of the nucleus has 
emerged which is quite different from that of an atom. 
In the latter, the small massive central nucleus inter- 
acts more strongly (%Z times) with any electron than 
that electron does with any other electron. For the 
nucleus there is no well-defined centre; the interac- 
tion of one nucleon with its neighbour at any instant 
of time is probably stronger than its interaction with 
the nucleus as a whole - this latter being a time- 
averaged effect of all the individual interactions. 
Thus the problem is a many-body one, without the ob- 
vious approximations of the atomic case; in the next 
chapter models will be illustrated which simplify the 
many-body aspects. The shell model has been taken 
over from the atom, but initially with reluctance, 
since its basic concepts seem to clash with ideas of 
the nucleus. 

In this chapter, the basic nucleon-nucleon inter- 


action is examined in order to answer the questions: 


(1) Assuming that the strong interaction between nu- 
cleons can be represented by a potential function, 
what is its basic shape and does it depend on the 


type of nucleon? 


(2) Arguing backwards from a nucleus towards a basic 
nucleon-nucleon interaction, does this inter- 
action agree with that deduced from basic nucleon- 


nucleon scattering experiments? 
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(3) Knowing that the answer to (2) is indefinite, 
is there a need to introduce new forces in the 
nucleus which are not present in the basic two- 


body interaction? 


The answer to the third question can be given 
now. Theorists are not happy with the idea of many- 
body forces, although there appears to be no very 
good reason why they should be small compared with the 
two-body interaction. The subject has not yet reached 
such a stage of refinement as to need a small 'topping 
up' with three-body forces. If they are to be intro- 
duced at all, it may as well be done in bulk - in 
which case what of four-body forces, etc? The hope 
appears to be that averaging many-body forces in a 
large nucleus produces much the same effect as avera- 
ging two-body forces. Note that, in answering this 
question, addition of a third nucleon will always 
modify the interaction between two nucleons because 
of the Pauli principle - but this effect is always 
allowed for in any attempt to solve the problem of 
larger nuclei and is not to be looked upon as a new 
three-body interaction in the sense referred to here. 


INFERENCES FROM THE DATA PRESENTED 

In developing the nuclear properties of the pre- 
vious chapter, a number of facts have emerged which 
have a direct bearing on the forces which hold the 
nucleus together. Taken roughly in the order in 
which they were presented they are: 


(1) nuclear sizes are of the order of a few fm; 


(2) nuclear density, at least near the centre, is 


the same for all nuclei; 
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(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 
(12) 
(13) 


(14) 
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protons and neutrons have similar density dis- 
tributions inside the nucleus; 


scattering of nuclear projectiles gives the 
Rutherford cross-section up to a certain angle 
with a rapid fall-off above this angle; 


the basis for a successful analysis of electro- 
Static energy differences between mirror pairs 
is that nuclear forces are ‘charge symmetric’, 
1.€. pp = nn; 

the nuclear density falls off over a range of 
~ 1 fm at the periphery; 

the leading terms in the binding energy, 

aA - ga2/s | were arrived at by analogy with a 
liquid drop; 


there is a ‘symmetry term' in the mass equation; 


the existence of a pairing term, and the coupling 
to zero angular momentum of every ground state 
of even A indicate a pairing force; 


the a-particle is almost as tightly bound as a 
massive nucleus; 


some heavy nuclei have large quadrupole moments; 
even the deuteron has a quadrupole moment; 


the proton and neutron have ‘anomalous' magnetic 
moments, i.e. different from those for structure- 
less Dirac particles; 


the magnetic moment of the deuteron is close to, 
but not equal to, the sum of the proton and 
neutron magnetic moments. 
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RANGE AND SATURATION 

The terms 'long-range' and 'short-range' have 
both a general and a specific meaning. In the general 
sense long and short are relative to the size of the 
system; van der Waals forces of cohesion of a liquid 
drop are short-range in the sense that the binding. of 
a given molecule is dominated by the contributions of 
a few fairly close neighbours rather than by the over- 
whelmingly greater number of distant molecules - these 
forces « l/r’ (or so) in energy-dependence. Likewise, 
in nuclei, forces are short-range if the contribution 
of most nucleons to the binding of a given one is less 
than the contribution from its near neighbours. In 
the specific sense, forces are said to have a range 
if the potential function contains an exponential, e.g. 
Yukawa forces having a potential function 
U « (l1/r)exp(-kr) are said to have a range 1/k. Thus 
it will be seen that nuclear forces probably are short- 
range in both the general and specific senses, though 
arguments on saturation are concerned really with the 
general sense. 

The nuclear force field extends beyond the nu- 
clear surface by a distance of the order of magnitude 
of the range of nuclear forces, so (1) and (4) (of the 
previous section) taken in conjunction can be used to 
define this range. The estimate obtained is confirmed 
by (6) since we would, in a simple way, expect the 
fall-off of nuclear density to occur over the same 
range. 

(2) and (7) in conjunction with the liquid-drop 
model from which the mass terms were derived do much 
to define nuclear forces. The molecular force field 


in a liquid drop has two main components: the attrac- 
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tive short-range force, which can be oversimplified by 
describing it as an induced-dipole - induced-dipole 
interaction for which the potential energy « ps and 
the repulsive, even shorter-range force which occurs 
when the electron clouds overlap. By analogy, a simi- 
lar state of affairs may be expected to occur in nuclei; 
the short-range attraction would require a range of 

1 - 2 fm, whilst the shorter-range repulsion might have 
a range of less than 0:5 fm. 

Oddly enough, in the historical development of 
the subject, the concept of a repulsion at short dis- 
tance was not readily acceptable. Nucleon-nucleon 
scattering had not yet indicated the need for such a 
force, which appeared to require a complex structure 
for what was hoped to be a simple fundamental particle. 
Instead, recourse was made to the exchange phenomenon 
(see Appendix C) to provide forces with suitable pro- 
perties. Such forces were described as: Wigner (W) 
or ordinary (no exchange of quantum numbers); Majorana 
(M) or space-exchange (corresponding to exchange of 
charge and spin); Heisenberg (H) or charge-exchange 
(but not spin); and Bartlett (B) or spin-exchange (but 
not charge). All these forces can arise from pion-ex- 
change. From analogy with electron-exchange binding 
of molecules, a nucleon should be surrounded by a pion 
cloud. That this is so can be inferred from the mag- 
netic moments of proton and neutron, which for simple 
Dirac particles should be l Uy and 0 Uy 3 but the pion 
has a much smaller mass than proton and neutron so the 
cloud can, and does, make a large contribution to their 
moments. 

Exchange forces will depend upon whether the pro- 
perty being exchanged is symmetric or anti-symmetric 
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in the wavefunction representing the relative motion 
of the two particles. As an example the spin-exchange 
operator (see Appendix C), es 2(1+4s,.85), has the 
value +l for S = 1 and -1 for S = 0, where 
on tS 8, + 85, SO the spin-exchange force can be made 
attractive for S = 1 and repulsive for S = 0, or vice 
versa. 

How do exchange forces lead to saturation? In 
a heavy nucleus A, the number of nucleon-nucleon inter- 
actions is 3A(A-1). If the forces are ordinary then 
all these are increasing the binding of the nucleus, 
and a complete calculation leads to all nuclei having 
roughly the same size with a binding energy propor- 
tional to a? rather than A. However, the operation 
of the Pauli principle ensures that all nucleon pairs 
cannot have the same spin symmetry or charge sym- 
metry or spatial symmetry, so that exchange forces 
will as often be repulsive as attractive. Their over- 
all effect need not be zero, since the strength of 
the force will depend upon a spatial overlap integral 
which of course is largest when the spatial wavefunc- 
tions of the two particles are the same. The exclu- 
sion principle, in its restricted statement, limits 
the number of particles with the same spatial wave- 
function to 4 - two protons (spin up, spin down) and 
two neutrons. Thus Majorana (spatial) exchange forces 
might be expected to produce strong interactions with- 
in groups of four particles and weak effects on aver- 
aging over pairs not in the same group. The emer- 
gence of four particles is an attractive result be- 
cause (10) emphasizes the fact that the a-particle is 
about as strongly bound as any heavy nucleus. 

Much effort has gone into producing a recipe for 
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nuclear forces using these four types - but each re- 
Cipe, whilst successful for explaining the particular 
fact on which the ‘cook’ was concentrating, has had 

its shortcomings, and the fact of a repulsive core 

(see p. 54) has been finally accepted with relief. 

This does not mean that exchange forces are unnecessary; 
it is difficult to see how ordinary forces with re- 
pulsion could lead to a saturated a-particle. Indeed, 
from the range of the repulsive force as required from 
nucleon-nucleon scattering, it is possible to attach 

an effective size to each nucleon (7 .¢¢ ~ 3 xX range, 
since two nucleons cannot approach closer than this 
range). The effective volume of the nucleons in a 
heavy nucleus is then found to be ~ 1 per cent of the 
nuclear volume. At its face value this indicates that 
Saturation is largely accounted for by exchange forces, 
though no doubt the repulsive core also plays an im- 
portant role. 


ISOSPIN 


In (5) (see p. 38) the point is made that forces 
are 'charge symmetric', by which it means that the 
(p,p) and (n,n) interactions are the same, apart from 
the electromagnetic contribution, but the (p,n) could 
be different. In the analysis of ‘mirror paris' the 
numbers of (p,p) and (n,n) interactions are inter- 
changed but the number of the nuclear-force binding 
does not specify the last interaction. 


A further step to 'charge independence' is made 
by assuming that all three interactions are the same; 
this can be checked by analysing nucleon-nucleon 


scattering or by a more detailed examination of nu- 
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clear properties. Charge independence is expressed 
mathematically using the isospin formalism, which. was 
developed as a convenience before it acquired physical 
Significance. Rather than referring to two distinct 
types of particle within a nucleus, theorists pre- 
ferred to employ a nucleon to which an additional two- 
valued quantum number was assigned and given the value 
+ 3 for the neutron and - 3} for the proton. The rea- 
son for this choice is that all two-valued quantum 
numbers have a similar mathematical behaviour, so the 
values reflect the similarity with pee 3 for the 
Spin quantum number. 

The name 'isospin' is a shortening of ‘isobaric 
spin', indicating that the two components have the 
same mass and emphasizing the similarity with spin S. 
By analogy, there exists a vector t and a component 
t. of isospin, but the space in which this vector 
exists is not the same as the space in which ordinary 
spin exists. Thus, in this space, te +} represents 
a neutron and t, = -3 a proton.’ Pursuing the analogy, 
a vector sum 72? = ¢, + ty; corresponding to S$ = s)+s,; 
can be made for two (and more) nucleons. A nucleon 
pair can have f = 1 giving a function symmetric on 
exchange of nucleons, or fF = 0, antisymmetric. If 
the Pauli principle is extended to make the overall 
wavefunction antisymmetric with the isospin component 
included then it can be seen that the same states in 


Tht is unfortunate that particle physicists have 


later chosen precisely the opposite convention - in 

principle it does not matter, since the direction up 
or down in an unspecified space is hardly relevant, 

but it is confusing. 
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the two nucleon configuration are obtained whether or 

not this formalism is used; since (p,p) and (n,n) are 

symmetric in isospin their wavefunctions must be anti- 
symmetric in (space) (spin), but the (n,p) combination 
can be symmetric or antisymmetric in isospin (f = l, 

T = 0 represents (p,n) and is symmetric; 7 = 0, 


a 
T = 0 also represents (p,n) and is antisymmetric) so 


the (space) (spin) wavefunctions can also be either. 
Thus. far, the formalism is merely mathematical, 
but charge independence can readily be fed into the 
scheme by asserting that nuclear forces depend only on 
T and not on To. This achieves its aim because the 
selection of 7 determines the type of (space) (spin) 
wavefunction and so ensures that we compare like with 
like when testing the similarity of (p,p),(p,n), and 
(n,n) interactions. The equality of strong interac- 


tions in mirror pairs of nuclei follows from this 


statement, at least for |v-Z| = 1, because these mirror 
pairs have Tf, = + 3 and, moreover, have the same value 
for T, also 3. Extending to the case |NW-2Z| = 2; ae 

14 


and O, corresponding to interchange of neutrons and 
protons, have the same value for T(= 1), but = i 
But charge independence further requires that there 
should exist a state in lay having T= 1, 7, = 0 with 
similar strong-interaction energy as the ground states 
of =—¢ and 6), The same should apply to excited 
states in rhe and aor The state of affairs is shown 
in Fig.2.1, where, to emphasize the point, the nuclei 
have been moved relative to each other through ener- 
gies to allow for the electrostatic term and the p-n 
mass difference. An interesting feature is that the 


14 


ground state of N has f = 0; there is a general ten- 


dency for states of lower fT to lie lower which re- 
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sults in the Symmetry term in the mass equation (point 
(8) of p. 38). The method used to derive it 


implied charge independence since it treated 
neutrons and protons in a similar way, which again 
is implied in point (3). 

Thus another quantum number has been introduced 
into the description of a nuclear state but, as the 
shift of energy in Fig.2Z.1 has perhaps indicated, it 


Fig.2.1. Comparison of energy levels in aioe My and ae 


Allowance has been made for the Coulomb energy and 

p-n mass difference. Note that re and 14, have T = 1 
states only, whereas an has T = 1 states and a greater 
number of 7 = 0 states. The dotted lines guide the eye 
to 7 = 1 states in 14y. (Based on AJZENBERG - SELOVE 


(1970). Wucl.phys. A152, 1.) 


is only an approximate one, since the electromagnetic 
term, whilst splitting the isospin sub-states by a 
large energy (several MeV), is also of a form which 
mixes states with different Tf to some (small) extent. 
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This quantum number plays a significant role in nu- 


16, a 14, ia 


N in excited states which are characterized 


clear reactions; e.g. the reaction 
1 14 
eaves 
by f = 0, since 16, (ground state), d, and a all have 
T = 0. Since the state at 2.31 MeV belongs to the 


T = 1 triplet which includes the ground states of G 


and 14 


O, it can be populated in this reaction only by 
way of the small impurity in the wavefunction of the 

T = 0 component. This selectivity in (d,a) and simi- 
lar reactions had been noticed experimentally and led 
to the conclusion that fF was indeed quite a good quan- 
tum number. 

Finally, to return to the mass equation, the 
pairing term (point (9)) appears to refute charge in- 
dependence in favour of charge symmetry, since the 
term increases the binding for (p,p) and (n,n) pairs 
but decreases it for (p,n). However, we are failing 
to compare like with like; the proton pair, or neutron 
pair, have similar quantum numbers - this will be more 
obvious after reading about the shell model in the next 
chapter - whilst in a large nucleus the (p,n) pair are 
in dissimilar orbits. However, in the very light 
nuclei when the (p,p),(n,n), and (p,n) interactions 
all refer to nucleons in the same orbits (shells), it 
is interesting to note that the odd-odd nucleus is 
stable on four occasions. Thus the pairing term can be 
shown to be consistent with charge independence. 


THE DEUTERON AND THE VIRTUAL DEUTERON 

Since the deuteron presents a simple two-body 
problem, its study was hoped to reveal the properties 
of the nucleon-nucleon interaction. In hindsight it has 


proved to be somewhat of a disappointment; the solution 
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for any simple short-range potential function which 

might approximate to the nucleon-nucleon interaction, 
e.g. the square well having potential - Vy fromr = 0 
to b and zero for r > b (see Fig.2.2) indicates that, 


Fig.2.2. Typical square well for the deuteron, depth 40 MeV, size 
1-6 fm (the scale factor for the abscissa). Also shown 
is the centrifugal barrier for 2 = 1, and u(r) = rF(r), 
where F(r) is the radial wavefunction of the ground 


state of the deuteron. 


in the loosely bound deuteron, the nucleons are for 
most of the time outside the effective range of the 
force - farther apart than b in the square well quoted 
above. To pursue this example it turns out that the 
wavefunction can be very similar for a large range of 
variation of the square-well parameters, provided that 
such variation keeps V 9b? approximately constant. The 


value of the constant is %v 10° -4Mev cm? i.e. 1 MeV b 
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corresponding to something like Vy v 50 MeV, b ™ 1°4 £m 
This condition follows from that fact that Kb ~ }n, 
where K is the wavenumber inside the potential, for 
successful matching of che internal and external wave- 
functions to produce a stationary state. For a more 
realistic potential without discontinuities, a simi- 
lar condition can be deduced in the form of a phase 
integral in the region of positive kinetic energy. It 
need hardly be added that the analysis of the deuteron 
sheds no light on the existence of a repulsive core. 
Also, because of the fact that the deuteron size stems 
mainly from that part of the wavefunction corresponding 
to negative kinetic energy and having a (slow) exponen- 
tial fall-off, feeding into the problem the (mean) 

size from electron scattering does not serve to pin 
down the force range with any precision. 

An interesting result arises, however. The lowest 
State is undoubtedly an S-state. That this state is 
bound only by 2-2 MeV, in a well of depth ~ 50 MeV, 
makes no other bound state possible; the next S-state 
with a node in the wavefunction requires a wavenumber 
of approximately 3 times that of the present state, 
which is clearly impossible, and the lowest 2 = 1 
state is made unbound by the ‘centrifugal potential’. 
Classically, a circular rotor with angular momentum 
mur ® has kinetic energy = ime which becomes 
HCL 1) (one on quantizing the angular momentum to 
{1(141) }a. A term of this form automatically arises 
in the radial Schrodinger equation - it is shown in 
Fig.2.2 forZ= 1. Numerically, to bind a state having 
2 = 1 would require a well-depth of about 209 MeV for 
the same range. It is highly unlikely that a mixture 


of exchange forces could be produced capable of such 
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large variation with Z of the potential. 

However, it is known from low-energy neutron 
scattering by protons that there exists a virtual state 
in the deuteron unbound by »v 30 KeV. The detailed 
analysis of scattering is beyond the scope of this 
discussion. Suffice it to say that from measurements 
of the differential cross-section it is possible to 
deduce the presence of resonances (or virtual states) 
and to determine the orbital angular momentum 7 taken 
in by the neutron in forming this virtual state. The 
present experiments were performed at thermal energies 
and could only detect such a resonance for Z = 0 (see 
later in chapter). Thus the virtual state is also 
an S-state. 

For the simple potential V(r) there should be 
Spin-degeneracy giving two S-state wavefunctions, 
cee That they are 


not degenerate is proof of the spin-dependence of nu- 


namely, J = 1, "Si and J = 0, 


clear forces and therefore of exchange - but the analy- 
sis does not lead to a recipe for the mixture. The 
ground state is obviously spatially symmetric and 

spin symmetric; in isospin it must be the antisymmetric 
T = 0 state. The virtual state is spatially symmetric 
and spin antisymmetric; in isospin it must be the sym- 
metric T = 1, ae 0 state. From the operator forms 

of the different exchanges (see Appendix C) it is seen 
that the Majorana term should be the same for both 
states, but that both Heisenberg and Bartlett terms 
should change. In the absence of any other states 

of the deuteron system, no evidence is available on 

the Majorana term; the effect of the other two terms 
would appear to be small, but the overall picture is 
far more complex than presented here so no quantitative 
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conclusions are drawn. 
So far, a simple theory has assigned the charac- 


teristics 3g J = 1, parity even - denoted by i 2 


1’ 
to the ground state of the deuteron. An S wavefunc- 
tion is necessarily symmetric in space, so the deute- 


ron should have no quadrupole moment. But, as stated 


in point (12) (see p. 38), the deuteron does have a 
(small) quadrupole moment. There can be no doubt that 
1* is correct since it is soundly based experimentally. 
The only other possible wavefunction having 1+ is the 
5Di-state, i.e. the coupling of spin 1 to4te=2 (nel 
has odd parity) to give J = 1. This wavefunction 
cannot by itself form the ground state since it has al- 
ready been shown that no state possessing orbital 
angular momentum can be bound by any reasonable 
potential, but in any case it would give too large a 
quadrupole moment. The remedy is to mix a small amount 
(v 4 per cent probability) of this wavefunction with 
the basic 3s) wavefunction. But this mixture no longer 
has a well-defined value for 4; L is certainly a 
constant of motion for any central force, and so a non- 
central term must be included in the potential. Such 
an interaction is already present. The two nucleons 
have magnetic dipole moments giving a dipole-dipole 
interaction of the form 

Vaa ™ p?{3r* (uy x) (uy -2) - (Hy -Hy)}, i.e. 

c pr 3{ 3" *(s,.1) (89-2) - (s,-S7)} = oe 842° This 
magnetic interaction is, however, quite negligible - 
something of the order of a few hundred eV. It is 
therefore necessary to introduce a strong interaction 
of the form V(r) S39? 
with a range, replaces the electromagnetic function 
7 This term is obviously different for the parallel 


where V(r), some radial function 
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spin orientations of Fig.2.3. Since the deuteron is 
prolate, the left-hand configuration must lead to an 
attractive force, giving a negative sign to the poten- 
tial term above. Note that the large quadrupole mo- 
ments of heavy nuclei (point (11)) result from cor- 
related motions of a number of nucleons and can arise 


with central forces between nucleons. 


D 
o 


Fig.2.3. Schematic representation of the spin-spin interaction 
in the deuteron. The spins are shown directed along the 
symmetry axis. 

This admixture of °D, in the wavefunction is also 
required in accounting for the magnetic dipole moment 
of the deuteron. As stated in point (14), the mag- 
netic eae is almost equal to the sum Uy + U3 for 
a pure S, “state it should be exactly equal. Classi- 
cally this statement is obvious; quantum mechanically 
it is much less obvious. A calculation along the lines 
of the Landé approach to hyperfine structure in atoms 
(see e.g. Kuhn 1961) is needed to produce this re- 
sult and also to determine the magnetic moment of the 
°p,-state. Adding the few per cent admixture of °D, 
moves the magnetic moment in the right direction. 


At this stage a note of caution should be injected. 
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If the magnetic moments of nucleons derive in large 
part from the meson clouds, what happens when nucleons 
are in proximity? For heavy nuclei this could lead 

to difficulty, but for the deuteron the nucleons are 
never really close enough for it to matter. 

From the fact that the admixture of =D; is small, 
it would be hasty to conclude that the non-central 
part of the potential is also small compared to the 
central terms. The 'centrifugal potential' for the 
D-state for such a small system serves to reduce its 
effectiveness to such an extent that it becomes im- 
portant to keep the range of the tensor force at least 
as large as that of the rest of the potential, in 
order to keep its magnitude down to an acceptable 
value. Both the exchange central force and the tensor 
force can be derived from pion-exchange. The fact 
that the pion has spin 0 but negative parity means 
that it must exchange with Z = 1 in order to preserve 
nucleon parities. This transfer of orbital angular 
momentum can cause the spin-flip necessary in some 
exchange terms and can also produce the tensor term. 
The exchange of a pion pair (07) results in pre- 
dominantly S-wave transfer and an ordinary force (the 
combination of two pions gives f = 2 or T = 0 sym- 
metric and T = 1 antisymmetric; these must combine 
with a symmetric and an antisymmetric space function 
respectively. Thus 7 = 0, S$ = 0, L = 0 is the only 
possibility for exchange). Since it is also repul- 
Sive it is believed to be (partly) responsible for the 
repulsive core. 


NUCLEON-NUCLEON SCATTERING 
The (n,p) scattering system at thermal neutron 
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energies has been cited (p.49) as a source of in- 
formation on the 255 virtual state of the deuteron. 
Scattering at higher energies with this system and 
the (p,p) system gives much more information pro- 
vided that the energy is high enough. Wavemechani- 
cally, the projectile is represented by a plane-wave 
passing the target, placed at the centre of the co- 
ordinate system. Such a plane-wave can be analysed 
into an infinite series of spherical waves converg- 
ing on and diverging from the coordinate centre: and 
characterized by the value of the orbital angular 
momentum Z. This series is the quantum-mechanical 
equivalent of the classical sub-division of the wave- 
front into zones characterized by the impact para- 
meter of the collision; if this latter parameter is 

p then pmv = h{Z(2+1)F is the correspondence, where 

v is the velocity of the projectile of reduced mass 
m. It is interesting to put numbers into this equa- 
tion, using the nucleon-nucleon system to define m: 
for Z = 1, 20 MeV incident neutrons in the laboratory 
system will correspond to a classical impact para- 
meter of 3 fm. The Z = 1 partial wave of the plane 
wave corresponding to 20 MeV neutrons will therefore 
be unaffected by the collision if the interaction is 
weak at a distance of 3 fm (the corresponding distance 
at FE, = 1 MeV is 13 fm). Thus there is a need for 
higher and higher projectile energies in order to 
determine the interaction for higher 7-values. Al- 
though thermal neutrons were sufficient for the S-wave 
interaction, they could give no information on the 
interactions at higher Z2-values for which laboratory 
energies ~~ 30 MeV are necessary. At Bee 300 MeV, 
the impact parameter is down to 0-7 fm for 2 = 1, and 
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so information can be obtained on the potential at 
close approach. 

Note that this same equation with (1(1+1) # =] 
is also an expression of the uncertainty principle; 
if the wavelength in the plane wave is A, then this 
distance is also the measure of, the uncertainty with 
which a particle can be located and is therefore a 
limit to which the system can be investigated. Thus 
for the S-wave interaction an energy of v 300 MeV is 
also required to probe the repulsive core. It is from 
experiments performed at these energies that the need 
arose for a repulsive core in the nucleon-nucleon 
interaction. 

Another result can be discerned from a qualitative 
examination of neutron-proton scattering (see Fig. 
2.4). The scattering in the centre-of-mass system 


100 


— 
=) 


a(6){mb sr-? 


40 80 120 160 Q 30 60 90 
8(centre-of-mass) (deg) @ (centre-of-mass) (deg) 


Fig.2.4. (a) n-p scattering at high energies. (b) p-p scatter- 
ing at high energies. Note, because of similarity 
of striking and struck particle, the angular distri- 
bution must be symmetrical about 90° (centre-of-mass) . 
(Based on HESS (1958). Rev.mod.Phys. 30, 368.) 
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rises sharply toward 0° and 180° at high bombarding 
energies. Peaking at 0° is to be expected; scattering 
through large angles requires a head-on collision (or 
nearly so), whilst small-angle scattering can occur 
when the neutron passes the proton at a distance and 
so experiences the extreme tail of the potential. 

Many more neutrons will 'miss' the proton than ‘hit' 
it. This naive picture conforms with analysis pro- 
vided the potential is simply V(r), i.e. an ordinary, 
central potential. The rise at 180° is proof of 
charge-exchange, corresponding to a small-angle de- 
flection coupled with exchange of charge; the detected 
neutron, which was originally the proton, moves in 

the centre of mass in the opposite direction. The 
type of exchange has been described as charge (or 
Heisenberg); it could equally have been charge-spin 
(or Majorana), since the measurements shown refer to 
unpolarized target and projectile and therefore are 
averaged over spins. As has been stated previously 

a large Majorana term is already required. 

Another point arising from nucleon-nucleon scat- 
tering is that the charge-independence hypothesis 
appears to be wellfounded. This is not a result which 
can be arrived at qualitatively; indeed the quali- 
tative conclusion (see Fig.2.4) might well be the 
opposite. The (p,p) system is, of course, affected 
by Coulomb forces, and there is an added complication 
that the two particles are identical, giving rise to 
interference terms on correctly antisymmetrizing the 
wavefunction. Even when these differences are taken 
care of, the comparison of (p,p) and (m,p) is still 
not like with like: the (p,p) system must have 
T = 1 whilst the (n,p) has equal admixtures of 7 = 1 
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and 7 = 0. The analysis is lengthy and complicated, 
but, the final result appears to be that the 7 = 1 con- 
ponent of (n,p) behaves similarly to the (p,p) when 
allowance is made for Coulomb and symmetry effects. 

One final piece of information has arisen from 
nucleon-nucleon scattering with polarized beam and 
target. A spin-orbit coupling term (L.S) appears to 
be necessary. Again it has a magnetic analogue, but 
again it is far too weak. In introducing a strong 
interaction of this type, the radial function has been 
given a short range since the term seems to become 
more effective at higher energies. 


SUMMARY 
Nuclear forces have been shown to have: 


(1) a short range, of order 1-5 fm; 


(2) a considerable exchange component, largely of the 


Majorana type; 
(3) a repulsive core, of range v 0:5 fm; 
(4) a tensor component; 
(5) a velocity-dependent spin-orbit interaction. 


The greater part of these requirements can be met 
by a potential derived from exchange of a pion, giving 
a range of 1-4 fm and accounting for the exchange and 
tensor components. The repulsive core can be accoun- 
ted for by simultaneous exchange of two (or more) 
pions or by exchange of heavier mesons. Among these 
heavier mesons must be included a vector meson 
Cal = 1) in order to account for the spin-orbit term. 

A rough qualitative picture of the nuclear po- 


tential is given in Fig.2.5. Notice that in detail 
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MeV 


Fig.2.5. An indication of the nucleon-nucleon potential, with 
hard-core and pion-exchange forces. Note that, in 
detail, the potential depends upon the combination 
of 7,2,S,¢. 


it should be different for different combinations of 
the quantum numbers L, S, and T. 


PROBLEMS 
2.1. If the n-p potential is of the form 
V,@) + Vj(r)s)-So> show that ce parity, as 


and S“ are constants of the motion. 

2.2. Determine the magnetic moment of the deuteron 
in the 20: 
admixture of this state into the °g-state gives 


-state and show that a few per cent 


the measured magnetic moment. 


2.3. Show that the formal introduction of isospin to- 
gether with a generalization of the Pauli prin- 
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ciple does not affect the number of states 
available to two nucleons. 


2.4. The Serber potential is of the form V(r) (1+P*) 
(see Appendix C). Show that it gives no force 
in odd spatial states. 


2.5*.The ground states of vBe and rac lie respec- 


tively 0°56 MeV and 3-58 MeV above the 3° ground 


state of ee Estimate the energy and spin of 


one excited state of reB i Comment upon the pos- 
sibility of exciting this state by deuteron in- 


elastic scattering. 


2.6*.Discuss what facts you can deduce about the 
interaction between two nucleons from the fol- 
lowing observations: 

(a) in the elastic scattering of neutrons by protons 
the angular distribution of the scattered neu- 
trons is nearly isotropic in the centre-of-mass 
system for incident energies below 10 MeV; 

(b) when the incident energy is raised above 100 MeV, 
sharp peaks develop at 0° and 180°; 

(c) the binding energy per nucleon in medium-weight 
nuclei is approximately constant and independent 
of mass number; 

(d) the cross-sections for scattering of slow neu- 
trons by ortho- and para-hydrogen differ by a 
large factor. 


3 s Nuclear models 


In the previous chapters, simple nuclear models have 
been introduced, namely, the liquid-drop and Fermi-gas 
models, which are in effect classical and quantum as- 
pects of the same system. This chapter will discuss 
three models: the shell model, the collective model, 
and the optical model. The first is an extension of 
the Fermi-gas model, whilst the last is a dynamical 
model concerned with the nucleus at highish excita- 
tion, where the restriction of the Pauli principle 

is decreasing in importance and the system is under- 
going transition to something. like a liquid drop. In 
the collective model the restriction to spherical 
symmetry is removed and the nucleons correlate their 
motion in the well in such a way as to create the 


non-spherical well in which they are moving. 


THE SHELL MODEL 

The atomic shell model has had a great deal of 
success in accounting for the construction of atoms 
and the chemical properties of the elements. A brief 
description is therefore appropriate. 

Solution of the Schrédinger equation for. the 
simple system of one electron bound by a central nu- 
cleus of charge Z units results in energy levels charac- 
terized chiefly by the principal quantum number n. 
Neglecting the fine structure, each level is degenerate 
having orbital angular momenta Z = 0,1,...(n-1), which, 
coupled with spin degeneracy 2, gives a total degeneracy 
Dn: Also the scale factor of the radial wavefunction 
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is largely determined by n. In filling up such an 
atom with its Z electrons, account must be taken of 
the Pauli principle. Thus the first two electrons go 
into the lowest energy orbits (n=1) and couple to 
zero angular momentum; their combined wavefunction 

is therefore spherically symmetric. The next eight 
electrons go into orbits with n=2, which are less 
tightly bound. They are therefore on average much 
further from the nucleus than the other two, and in 
addition six of the eight have one unit each of angu- 
lar momentum which tends to keep them away from the 
centre. Thus these electrons rarely penetrate the 
spherical cloud which represents the other two and so 
move around a central nucleus of effective charge 
(Z-2). This electron group corresponds to a distri- 
bution of charge in the form of a broad spherical 
shell (in the classical sense), overlapping very 
little with the shell below. So we build up a complex 
atom shell upon shell, rather like an onion. This 
picture is oversimplified and leads to incorrect num- 
bers of electrons in shells Coat The reason is that 
the interpenetration of shells does take place. Be- 
cause of the centrifugal effect sub-shells of lower 

Z will penetrate the shell below to a greater extent 
than those of higher Z (and the same value of 1), so 
removing the degeneracy enough to move sub-shells 
effectively from one major energy shell to another. 

A complete sub-shell, however, is spherically sym- 
metric, so the overall picture is much the same but 
the numbers of electrons in major shells is now dif- 
ferent. This has a profound effect on the chemical 
properties of the elements. 


In trying to adapt this approach to the nucleus, 
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a number of difficulties arise. In both cases the 
potential well for the next particle must be built 

up from the interactions of particles already present. 
In the atom there is initially a strong attraction 

to the centre, which carries virtually all the mass 
and provides a basis for building up layers, in ad- 
dition ensuring that the layers are physically well- 
separated. A central force is therefore dominant at 
all stages, since a uniformly charged spherical shell 
behaves, in the region outside it, as though the 
charge is concentrated at its centre. From Chapter 1, 
however, it is known that there is nothing special 
about the centre of the nucleus except that it is in 

a region where the density is constant and where, from 
symmetry, the average force on a nucleon is zero. 

Thus the simplest form of potential which might approx- 
imate to the average nucleon-nucleus potential is 

a square well, whilst the next approximation is to 
make the potential well follow the shape of the den- 
Sity (see Fig.3.1) but extending further out by 1 fm 


Fig.3.1. Saxon-Woods (or Fermi) approximation to the nucleon- 


nucleus potential. 
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or so, to account for the range of the nucleon-nucleon 
force. With such a well the wavefunctions for the 
single-particle eigenstates overlap strongly. The 
centrifugal potential pushes the wavefunction away 
from the centre for Z # 0, but the effect is not mar- 
ked for 2 < KR, where K is the wavenumber of the 
nucleon inside the nucleus of radius #, and KR ~ 10 
for a medium-weight nucleus. 

In addition to greater mixing of wavefunctions 
enhancing the importance of collisions, there is also 
the fact that the nucleon-nucleon interaction is short- 
range so that its contribution to the average poten- 
tial is much less than its contribution to collision. 
It would seem improbable that the picture of single- 
particle orbits in an average central potential should 
approximate to reality. But so far no account has 
been taken of the Pauli principle. If one particle 
gains energy in a collision, then the other particle 
must lose it. In the ground state of a nucleus all 
the lowest-lying single-particle states will be oc- 
cupied, so no particle can lose energy. .The same 
two states must be occupied before and after the col- 
lision, so effectively there has been no collision. 
The mean free path of.a nucleon in the ground state 
of a nucleus is effectively infinite, and the concept 
of particle orbits is meaningful. 

Having established the picture of nucleons in 
well-defined states it is now appropriate to consider 
whether these states group into shells, characterized 
by having similar energies but not necessarily similar 
spatial configurations as in atoms. In the case of 
atoms, the periodic table preceded its explanation 
in terms of the shell model, so for nuclei stutistical 
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analyses of the properties of nuclear ground states 
revealed the 'magic numbers' which are the nuclear 
equivalent of the location of noble gases in the 
periodic table. We take a nuclear property such as 
binding energy B and plot it against A, or, better 
Still, make a three-dimensional plot against W and Z. 
Thus in Fig.1.2 (p.19) the measured values of B/A 
differ from the smooth curve predicted by the mass 
equation in the direction of increased binding in 
the region of certain values of W or Z. Statistical 


analysis reveals the following facts: 


1. A plot of binding energy of the last nucleon 
against W (or Z) reveals discontinuities at 
certain values of W and Z. (The neutron binding 
energy may conveniently be determined from the 


y-spectrum following capture of thermal neutrons.) 


2% Similar discontinuities in the binding of a- 
particles occur at W = 82 and 126. 


53 Certain nuclei appear to be naturally more abun- 
dant than their neighbours in the plot of stable 
species. Certain elements have more stable iso- 
topes than their neighbours, i.e. a greater varia- 
tion of W for a given Z, and correspondingly 
there are values of W for which a greater varia- 
tion of Z is permitted. 


4. The excitation energies of the first excited 
+ : ; 
states (2 ) of even-even nuclei are at maxima 
for certain values of W or Z, falling to broad 
minima between these values. In the broad mini- 
+ + 
ma the lowest states (0 , 2 


in which the energy spacings (0,2) and (2,4) are 


+ 
, 4 ) form a sequence 


in the ratio 6:14; near the maxima these states 
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tend to be equally spaced. 


Quadrupole moments of odd-A nuclei are large when 
neighbouring even-even nuclei have low-lying 
first excited states and are near zero when those 


states lie high. 


Nuclear reaction cross-sections reveal periodic 
properties, e.g. thermal-neutron capture cross- 
sections (as shown in Fig.3.2). 


I 
10 20 30 £40 50 60 70 80 90 100 110 120 130 
Number of neutrons 


Fig.3.2. Effective capture cross-sections for reactor neutrons 


Te 


plotted against neutron number of target nucleus. 
Detailed explanation of this curve is complicated, 
but the influence of shell closure is apparent. (Based 


on CODD et al. (1956). Prog.nucl.energy 1, 296.) 


The parities of nuclear ground states of odd-A 
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nuclei exhibit a systematic behaviour correlated 


with the above 'magic' effects. 


The nett result of such analyses is that nuclei 
containing 2, 8, 20, (28), 50, 82, 126 neutrons or 
protons (apart from 126) are particularly stable 
(cf. the stability of noble gases). By analogy, shell 
closure at these neutron or proton numbers is there- 
fore expected. The single-particle energy levels of 
nuclei should reflect such closure by a larger than 
usual energy spacing at appropriate places in the 
level diagram. Simply altering the shape of the cen- 
tral potential through permitted variations is unable 
to produce energy gaps at these numbers (see Fig.3.3), 
where the ordering of levels for four types of well 
is shown to be unsuitable but the addition to the 
potential of a spin-orbit term (= V(@)Z.s for a given 
nucleon) achieves success. This suggestion was made 
independently by Mayer and by Haxel, Jensen, and Suess 
in 1949. Such a term splits the otherwise degenerate 
sub-shells 7 = 2 + 3 and g = 2 - 3} obtained by combin- 
ing the spin and orbital angular momenta of a nucleon. 

Since 


then 


(s.l) = 1{7(g+1) - 1(4+1) - s(s+1)}, 
giving the value +327 for j = 2 + 3} and -3(2+1) for 
j = 24 - 4, except for Z = 0, when both terms are zero. 
Thus the splitting « (22+1); its increase with 7 is 
important for the success of the modification since 
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——————————————————— 
Harmonic Infinite Finite Square Plus 
oscillator _ square square well with spin-orbit 
well well rounded coupling 
edges 


Fig.3.3. Locations of energy levels in four types of well, fol- 
lowed by the effect of a spin-orbit term on the fourth 
type. Note that the spin-orbit term is necessary to 
locate the magic numbers, as shown. (Based on FELD 
(1953). A.Rev.nuel.Set. 2, 239.) 


it serves to throw down, from one shell to another, one 
of the components of the sub-shell with highest 27. 
By making the sign of the spin-orbit coupling opposite 
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to that arising from a magnetic term (as in atoms) the 
sub-shell corresponding to j = (l2+3) is correctly 
thrown down. The order of the filling of neutron 
shells is then given by the sequence 


Z Z 


(1s4)?1 (ap) *(ap,)?| (14g) ©c2s,)7 (145) 4] c1,) 8! cap.) 
2 2 2 2. 2 2 Z 2 
2 8 20 28 | 


ie Gp.) Gag) | dyes) (as, - Cd ns) ee) 
Zz 8 2 z z z > 2S Z 
50 82 


(1hy)*°(3p,) 4(2£,)9 (3p) 7 (1844) |(2e9)" (3de)° C14, 4) (289) 
z i i i a a a a ee: 
126 


This sequence can be varied slightly by small changes 
in the spin-orbit potential without affecting the 
magic numbers. 

The notation is slightly different from the atomic 
convention. The Z. convention is the same but the 
number preceding does not have the significance of n 
in the hydrogen-like spectra; it merely numbers the 
states of given Ls in the order in which they occur 
(it can be connected with the number of nodes in the 
radial wavefunction). The superscript is the number 
of nucleons which the sub-shell can accommodate (2j+1). 
For protons, because of the influence of the Coulomb 
potential, the ordering is slightly different; it is 
the same up to 50 when the next shell looks like 
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(les) 26)" (nia) de) Gee) 
2, 2. ra 2 bX 


and above 82 the sequence is 


Cina Eo) pa) eaartes 
Z 2 2 

These sequences, together with the pairing tern, 
which is taken to couple pairs of like nucleons to 
zero angular momentum whenever possible, account for 
all the parities of the ground states of stable and 
nearly stable nuclei and, with a few exceptions, their 
spins for all even-A nuclei and up to A v 100 for 


odd-A nuclei. Two notable exceptions are as: and 23Na 


3+ 


which are and whereas the sequence of protons 


predicts d,, l.e. 3* for both. However, the Sy> d.- 
a. oe 


states lie in the same shell and it was also discovered 


SS) Sa Se) fae 
+ 


that a 3" excited state lies low in each nucleus, so 


it was felt that 'residual forces' (i.e. forces left 
over after accounting for the well by averaging the 
nucleon-nucleon interaction) were moving these sub- 
shells relative to each other. Lately it is considered 
that nuclear deformation is responsible for these dis- 
crepancies. 

Above A ~ 100 the ground states of odd nuclei do 
not appear to reflect the order in which the shells 
should fill; e.g. four isotopes of Sn (Z = 50) have 
ground states > instead of only one as predicted. 

The reason lies in the pairing term; adding a second 
neutron into one of the isotopes a should produce 
the pair Gas but the pairing of neutrons in a 


a 
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higher lying state with higher Z, Probably the lh 1 


2 
state, could give an increased binding which more than 


off-sets the higher energy of the intrinsic state. In 
which case the pair of neutrons will be found in the 
1h,,-state in the next even-even isotope, with the 

a | 
3s,-State unoccupied. So the next odd isotope of Sn 


pole 


will again be a etc. 


Also, just below this region, ground states is 
should be found corresponding to neutron odd numbers 
57 - 63. None is found, but in each nucleus where 
expected, a low-lying is state occurs. A Similar ex- 
planation can probably be made, though it is possible 
that the ordering of the lg.- and 3s,-states should be 


2 Z 
inverted. For the more massive nuclei, nuclear de- 


formation is expected to occur away from the closed 
Shells, and this would produce a different ordering 
of the ground-state spins (see next section). 

A triumph of the shell model concerns the loca- 
tion of ‘islands of isomerism'. It had been noticed 
that long-lived isomers (excited states of nuclei 
which decay by y-emission or its equivalent, internal 
conversion, with lifetimes long enough to be classi- 
fied as radioactive) tended to cluster at certain 
regions of the periodic table. To understand this it 
is necessary to know something about y-ray selection 
rules (see Chapter 5). At this stage we merely need 
to know that the bigger the spin change the slower 
the transition. In certain regions of the periodic 
table, the first excited state is expected to differ 


in spin from the ground state of an odd-A nucleus by 
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a large amount. For example, at Z = 39 the last pro- 
ton is ina 2p,-state, whereas the first proton ex- 


2 
citation is 1g; it lies at 0-381 MeV in aan at 


4 
0-91 MeV in B9y and at 0°551 MeV in 
lifetimes of 14 hours, 16 s, and 50 min respectively. 


Sly, giving 


Thus it is seen that the spin-orbit term which pro- 
duces the right sequence of shells also puts sub-shells 
of high spin in juxtaposition with sub-shells of low 
spin and so is responsible for isomerism. 

Some mention must be made of the magnetic moments 
of the odd-A nuclei. From the assumption of maximum 
pairing, the properties of the nucleus should be those 
of the odd nucleon. The magnetic moment contributed 
by a single nucleon is 


m 
uo= 9, (Lj) + gf s.5) FIGEEDI with m4. = d; 


where 2(¢2.j) = g(j+1) + Z2(2+1) - s(s+1), and corres- 

pondingly for ‘(s.j). This can readily be derived 

using the vector model (see e.g. Kuhn 1961). By and 

large the agreement with experimental data is poor. 

It is best when the odd nucleon immediately precedes 

or follows shell closure; examples for the odd proton 
3,, 15, 39,41, 89 

are 18> aN; 19*> 39%» and for the odd neutron 


3 17 207 41 45 91 | 
jHe, g0> g2Pb, but even here 99f8; 215°» Ag2T show 


large deviations. ‘SE might have been quoted above 

were it not for the fact that theoretical computation 
has come up with an untidy, three-particle mixture of 
S4> d., and d, wavefunctions outside an 16 


a. ae Zz 
that happens to give the magnetic moment of an s 


O core, 


1 Pro- 
2 
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ton. The deviations of magnetic moments are as yet 
unexplained; they could arise from three sources: (1) 
strict pairing does not occur, (2) nuclei may not be 
spherical, and (3) the pion field around a nucleon in- 
side a nucleus should differ from that around a nucleon 
in isolation; this could give rise to a magnetic moment 
differing from that obtained by vector addition of 
the individual moments as ascribed to isolated nucleons. 
Finally, how does the shell model cope with ex- 
cited states? The general ordering of levels indicated 
in Fig.3.3 should represent filled levels up to a cer- 
tain value and single particle excitations beyond this. 
If strict pairing does occur, the excited states of 
an odd-A nucleus should follow this sequence of ex- 
Citation. Such single-particle excitations do occur, 
especially when they refer to large Z-values in regions 
of sub-shells of low Z. However, the number of ex- 
cited states up to, say, 3 MeV above ground is usually 
far greater than the number of single-particle levels 
available. We are forced to conclude that strict 
pairing has broken down (it may well take only 
~ 1 MeV to break a pair) and that extra states are ob- 
tained by the coupling of three or more particles. An 
extreme example of this is a Its first excited 
state, at 108 keV, is > whilst the lowest single- 
particle states of odd parity (1f£,), (2p) should be 


2 Zz 
several MeV higher. Merely decoupling the last two 


neutrons will not produce odd parity since all three 
particles are in even-parity orbits - the Sy> dc; and 


2 Z 


d., which lie close together. It has therefore been 


[ur 


sumed that the closed proton p-shell of the +26 core 


"M bo 


a 
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has been disrupted, and one of the protons has gone 
into the s-d shell to create something like an a- 
particle; the (p,) "hole' provides the odd parity for 


2 
this state and others at somewhat higher excitation. 


Thus, where excited states are concerned, even the 
closed shells must be considered vulnerable. 

The shell model forms a working basis for attempt- 
ing to determine the level scheme of a nucleus, but 
there are interactions between the nucleons, termed 
"residual forces' since they remain after accounting. 
for the well by averaging out the main body of the 
forces. These residual forces are often arrived at 
empirically in building up nuclei progressively from 
a double-closed shell core by comparison with experi- 


mentally determined level schemes. 


COLLECTIVE MOTION 

In this and the previous chapters, reference has 
been made to a collective nuclear motion. For this 
to occur, motion of closed shells is implied, and this 
in turn means that such shells cannot be spherically 
symmetric - quantum mechanically there can be no mo- 
tion about an axis if the system has symmetry about 
that axis. Thus the occurrence of collective motion 
is characterized by a departure from spherical sym- 
metry of the nuclear system; it may be the rotational 
motion of a permanently deformed nucleus or the motion 
of a nucleus vibrating (possibly about an undeformed 
mean structure) through different degrees of deforma- 
tion. 

The possibility of permanent deformation was in- 
vestigated by Nilsson, who calculated the single-par- 
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ticle states in a spheroidal potential well having a 
Symmetry axis. The distortion of the spherical well 

1s then characterized by the single parameter & = AR/R, 
as defined in Chapter 1. For such a well ftxed tin 
Space it must be appreciated that, because of the ab- 
sence of spherical symmetry, the total angular momen- 
tum J is no longer a constant of the motion, but be- 
cause the system does have symmetry about one of its 
principal axes the component of angular momentum along 
this axis, designated by K, is indeed a constant of 
the motion. The behaviour with 6 of the low single- 
particle levels is shown in Fig.3.4, where the value 

of |K| is shown on each curve. Since the deformation 
is quadrupole no distinction is made between up and 
down, so that each curve has degeneracy 2 corresponding 
to + K. In the limit 6 = 0, the spherical states 


appear; in particular the d,-state appears as the con- 


fluence of three states of different |K| and there- 
fore has degeneracy 6, which is in accord with its 
spatial degeneracy, 2j7 + 1. 

If the problem is made physically more realistic 
by isolating the well in space rather than fixing it, 
then the well itself will be capable of rotation, but 
the rotation must have no component along the symmetry 
axis. If the rotational motion is slow compared to 
the single-particle motion (the adiabatic limit), then 
the latter will have a well-defined angular-momentum 
component along the nuclear axis K, but g itself will 
not be well-defined. The total angular momentum J will 
also be well-defined (an isolated body) and will have 
a component K along the nuclear axis and a component 
M along the external axis of quantization. J will be 
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Single-particle states in a fixed ellipsoidal potential 
well, characterized by 6 = AR/R. The states are 
labelled with |X| and parity (they are doubly degenerate) 
when 6 # 0. At 6 = 0 they correspond to the sub-states 
of the usual shell-model states. The +s, -state has 


been omitted. (Based on NILSSON (1955) .° Doan.mat. 
-fys.medd. 29, No.16.) 


a vector sum of j and R, the nuclear rotation (see 
Fig.3.5). The rotational energy will be given by 
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= (fixed in space) 


(fixed in 
TS See nucleus) =’ 


Fig.3.5. Coupling scheme for a deformed nucleus. J = j +R, 


der = Jar = kh, J. = Mh. Note that j itself is not 
a constant of motion, nor is 7? well defined; e.g. 
T 7- : 

J = z «can be constructed from fo hg» hij orbitals 


coupled to # = 0 or 2, and, except foe extremes of 


coupling, all components will be present. 


a = (h2/21)¢R2> = (n°/27)¢ (J-j)*> 


try 
t 


(h7/21) (J (J41) + (9%) - I-49) 


(h2/21) (s(s41) + (72) - 2K%), (3. 


since it can be shown that 


(Jd + Sap = 0, 


1) 
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unless there are components. of K + 1 as well as XK in 
the wavefunction. Since there are components of + K 
due to up-down degeneracy, these terms are not zero 
for |K| = 4, but this case will be ignored. The term 
in (32) can be incorparated into the other single- 
particle contributions to the energy, and forgotten. 
This expression defines a set of rotational states 
based upon a common wavefunction relative to the nu- 
clear axis. The lowest-lying state will have J = K 
(the lowest value J can have for a component K along 
the nuclear axis) and a degeneracy 2K+1. So, for the 
isolated nucleus, the confluence at 6 = O for the 


7 = d, configuration will no longer have the correct 
Pa 
number of sub-states and would appear to have three 


values of ee > >) instead of the 3 state of the 


spherical nucleus. 

The reason for the discrepancy is that the rota- 
tional energy has been ignored in Fig.3.4. How to 
evaluate this depends upon the nature of the rotation. 
It could be a rigid-body rotation, or an effective 
rotation of a hydrodynamical fluid, an example of 
which is the tidal deformation of the earth by the moon 
On the surface of the earth this appears like a rota- 
tion, but the motion is in fact up and down. Experi- 
mental investigation indicates that the latter is 
nearer the truth and I % 18": where Ty is 2a a 
body moment of inertia of a sphere, namely, =p MR. 

The rotational energy therefore becomes very large at 
small 6 so the states J = j and ; produced from a d. 
particle move up to very high energies as 6 eeonee 
small. More correctly, the coupling scheme breaks 
down; it is no longer valid to couple the single-par- 
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ticle motions rigidly to the well and then to allow 
the whole system to rotate. The figure only approxi- 
mates to reality when 6 is sufficiently large for the 
rotational energy to be small compared to the spacing 
of the different particle configurations. 

At sufficiently large 6, the diagram correctly 
predicts the ordering of states. The particles couple 
strongly in pairs with equal and opposite K, so an 
even-even nucleus will always have K = 0 and a ground 
state J = 0. Putting J = 0, K = 0 in the rotational 
energy gives a rotational band of levels of energies, 
relative to the ground state, J(s+1)h7/21, where be- 
cause of symmetry requirements (up = down). only even 
values of J can occur. This symmetry requirement is 
met by expressing the wavefunction of the system as 
a combination of the wavefunctions (J,M,+K) and 
(J,M,-K) thus: (Yn ex + (ane = Uc) For K = 0 
it is seen that the two components are the same except 
for the term (21)7 which eliminates all odd J. For 
other values of K no such restriction arises, so J 
increases in unit steps up a band. For odd-A nuclei, 
the sequence of levels in Fig.3.4 determines a value 
of |K| for the ground state. A rotational band is 
now built up having energies relative to ground 
h*{o(J+1)-K(K+1)}/21. It should be noted that K=} pro- 
duces a more complicated ordering of states which will 
not be discussed, except to mention that it arises from 
the (J.j) term in the rotational energy. 

A typical level scheme is shown for an even and 
an odd nucleus in Fig.3.6. It will be seen that ro- 
tational states up to large values of J are obtained 
and that they tend to be lower lying than the next 
particle configuration. Their locations usually con- 
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Fig.3.6. Examples of rotational bands in even-A and odd-A 
nuclei. In the latter nucleus the level schemes 
representing the two non-interacting bands have been 
separated laterally. For clarity some levels un- 
connected with the bands have been omitted from both 


schemes. 


form well with the simple expressions above, but there 
is usually a small discrepancy which increases with J 
and which is ascribed to an increase of I due to a cen- 
trifugal stretching of the nucleus. Experimentally, 
levels in a rotational band are revealed by the very 
large E2 transitions (see Chapter 5) between members 

- the intrinsic strength can be as large as 200 times 
single-particle strength in heavy, deformed nuclei 

in the rare-earth region. Such large enhancements 
arise from the rotation of the nucleus as a whole, 
corresponding to a correlated motion of up to Z pro- 
tons and therefore capable in principle of producing 
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enhancements up to z* 


Finally, nuclei undeformed in the ground state 
can have excited states corresponding to vibrations 
through varying degrees of deformation. As in the 
case of the harmonic oscillator, the different vibra- 
tional states are equally spaced, but because the de- 
formation is quadrupole, the vibrational quanta are 
characterized by = giving for an even-even nucleus 
the sequence of states (0), (2°) (0. 2° hy 
(092 i) ed 6 etc. where the brackets correspond 


to 0, 1, 2, 3, etc. vibrational quanta coupled using 


106 
jr aoPt keV 
4+ 1228-9 


0+ 1133-3 
oF 1127-8 


+ 511-8 


or 


Fig.3.7. A vibrational system. The next levels are located 
around 1-6 MeV and could include members of the three- 


quantum vibrational group of states. 
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Bose statistics. The bracketed states are degenerate 
in the model, but departure from the model removes the 
degeneracy as Fig.3.7 shows. 


THE OPTICAL: MODEL 

The previous models refer to a static or quasi- 
Static system since the excited states so far intro- 
duced are bound states decaying only by the relative- 
ly weak radiation process and therefore having life- 
times much longer than a characteristic nuclear time, 
which is the time taken for a nucleon to traverse 
a nuclear diameter. In this situation the quantunm- 
mechanical problem consists of setting up stationary 
states (just as for the ground state) and then treating 
the radiative process as a perturbation. When the 
energy of the system is large enough to permit nucleon 
(or composite particle) emission then the decay is 
via the strong interaction so the state can have a 
lifetime of the order of the characteristic time. 
The decay process is no longer a small perturbation 
and a different approach is needed. The optical model 
meets this need. 

This model attempts to cope with the quantum- 
mechanical problem of the scattering and absorption 
of particles impinging on a nucleus. Since a simple 
potential has been used to represent the nucleus in 
the shell model, it is appealing to do the same for 
the scattering problem. An argument for the use of 
such a potential was based upon the long mean free 
path of a nucleon inside the nucleus as a consequence 
of the Pauli principle. However, a nucleon coming in 
from outside has (ELan + 7) MeV (the incoming kinetic 


energy + binding energy) more energy than any nucleon 
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within the nucleus. It can collide with such a nu- 
cleon thereby losing energy and exciting that nucleon. 
Because of the Pauli principle it cannot interact with 
nucleons lying more than (Flan + 7) MeV below the 
Fermi surface of the nucleus. Thus as the kinetic 
energy of the incident nucleon increases, the limita- 
tions of the Pauli principle are reduced and the mean 
free path of this nucleon gets shorter. 

How is this concept, fed into the problem? A beam 
of particles traversing a region in which the mean 
free path of a particle is X will decrease in beam 
intensity with distance as I = I exp(-x/A). If the 
beam of particles, without absorption, is represented 
by the plane wave amplitude VI) exp(ikz), then with 
absorption it should be 


VIq exp(-x/2\) exp(ikzx) = VI exp{i(k+i/2A)z}. 


The wavenumber has therefore. been moditied from k to 
(kK + i/2d). But hk = {2m(E-V)}?7, where F is the total 
energy of the particle and V the potential energy. 
Putting in the modified expression gives 


Le EPID) <a SES Oe 


k 
Since the total energy of the particle must be real 
(it can move outside the region of absorption to be 
measured), an imaginary component must be added to 
the potential. Expressing this as V+iW gives 
2m(E-V) /h@ n?-(1/2a)7 me x? for reasonably large A; 
and 2m /h? -k/X connects the imaginary component 
with the mean free path. 


In its early stages therefore the optical model 
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used a complex square-well potential, -(V+iW) for 
ry < R, zero forr > R&R. This was later changed to the 
Saxon-Woods potential 


-(V+ivW) 


l+exp{(r-R)/e} ’ 


where V, W, R, and ec were looked upon as arbitrary 
parameters in fitting the experimental scattering data 
at a given energy. The latter two parameters should 
obviously correspond in some way to the density para- 
meters; the simplest way being to assume that the 
potential follows the density when allowance is made 
for the range of nuclear forces. The real term JV 
might be expected to approach the shell-model potential 
at low enough energy. There is no good reason why V 
and W should behave in the same way; indeed, an energy- 
dependence of W is expected from the previous argument 
which is unlikely to be similar to that of VY. Also 

it can be argued that the bound nucleons near the nu- 
clear surface should be those with greatest energy, 
for which collisions with the incoming nucleon are 
least restricted. Thus a peaking of W near the sur- 
face may be more realistic, and this can be achieved 
by associating W with the derivative of the Fermi 
(Saxon-Woods) function, or by simply using a Gaussian 
function peaked at the surface. In addition, for 
treatment of polarization data, it has been found 
necessary to introduce a spin-orbit coupling tern, 

« 2.5, where £4 and s refer to the incident nucleon. 

It is instructive to consider s-wave scattering 
from a square-well potential of radius R in the limit 
of very small W. In order to avoid the complexities 
of determining elastic scattering cross-sections, the 
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variation with energy of the reaction cross-section 
is estimated; for W=0, there is no absorption and 
therefore only elastic scattering, but for W small 
enough the simplifying features of W=0 can be used. 
The wavefunction inside the well (for W=0) is of the 
form (A/r)sinkKr , satisfying the requirement of a node 
at the origin in r times the radial wavefunction, 
whilst outside the well the solution can be written 
(B/r )sin(kr+6), where K and k are the interior and 
exterior wavenumbers: K = 2m(E+V)? 2 th and k = (2mE)? 2h, 
The standing-wave solutions indicate that no absorption 
has taken place, i.e. as many neutrons are moving away 
from the nucleus as toward it (and the flux of neutrons 
toward the nucleus « zn |B) 4, though the direct con- 
nection is complicated since the neutrons proceeding 
toward the nucleus constitute a plane wave). 

The continuity of the wavefunction and its deri- 


vative at the nuclear surface lead to the equations 


Ktan KR 


Ktan(kR+6) 


and 


Bk/(K*cos“KR + k-sin?xry 1/2, 


mS 
I 


If for simplicity it is assumed that K > k then A/B 
is seen to vary between maxima of ~ 1 when KR WV (n+3)T 
and minima of ~ k/K when KR ~ nt. Also lal? will fall 
to half its maximum value when KR ~ (n+}3)T7 + kK/K. 

If now a very small imaginary term is introduced 
into the potential, then absorption will take place 
within the nucleus, but for small enough W this will 


hardly change the wavefunctions. From the connection 
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between W and mean free path it is obvious that the re- 
action yield will be « w\a|?. The reaction cross- 
section will therefore be « wia|*/k|B|* = wk/(K2cos?KR + 
+k*sin“KR), showing maxima and minima with variation 
of neutron energy. If, as in Chapter 4, the problem 
were approached in a different way, namely, by in- 
serting the neutron inside the nucleus and observing 
its transition out, then it will be seen that the maxi- 
ma correspond to the formation of ‘virtual states' of 
the system. The reaction is said to exhibit 'resonance' 
when the neutrons are at an energy appropriate to one 
of the virtual states within the nucleus, cf. forced 
oscillations of a pendulum near its natural frequency 
or of an electric circuit. The ratio of resonance 
half-width to resonance spacing can be estimated from 
the above equations as (2k/K) x 1/m or (4k/K) x 1/2mT. 
The ratio has been expressed in this latter way since 
the factor 4k/K is the mismatch factor at the poten- 
tial step (see Chapter 4 and Problem 4.1), whilst 
the factor 1/27 represents an important theorem con- 
cerning widths of virtual states, namely, that on 
average the decay width of a state into any one channel, 
after allowing for matching the internal and external 
wavefunctions, should be 1/21 times the local spacing 
of levels with the same spin and parity. For a medium- 
sized nucleus and a well-depth of 50 MeV, the level 
spacing of s-wave single-particle neutron states 
(7 = y* for J" = 0° targets) is about 15 MeV and 
the single-particle width about 2 MeV. Thus the single- 
particle resonances are broad but well separated. 

In going to a more realistic model of the nucleus, 
as well as the single-particle states there are also 
found vast numbers of levels of the same spins and pari- 
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ties but corresponding to excitations of more than 
one particle. In fact the single-particle states 
alone cannot produce reaction but only re-emission of 
the same particle. That reaction does take place in- 
dicates that the single-particle state is mixed in 
(by some unspecified perturbation) with all these 
more complicated states and W must in some way re- 
present the degree of mixing. In the mixed-up system 
the theorem concerning level width and spacing will 
on average still hold (see below), though the mixing 
may reduce the spacing from a few MeV to a few eV. 
Inclusion of a large imaginary term in the po- 
tential will broaden out the resonance even more (the 
absorption term is representing other channels into 
which the state can decay, so the width should rough- 
ly increase by the amount W (see Chapter 6)). Since 
absorption has taken place, followed by decay into 
another (reaction) channel, a broad resonance is to 
be expected in this reaction channel, e.g. the emis- 
Sion of inelastic neutrons. Data on total neutron 
cross-sections are presented in Fig.3.8 and also the 
results of analysis in terms of the optical model. 
A feature of the model is that the parameters which 
describe the potential should vary only slowly from 
nucleus to nucleus, and this appears to be the case. 
Typical values for the potential for neutron scat- 
tering from a wide range of nuclei at neutron ener- 
gies up to ~ 50 MeV are 


Voe= 52:5 - 0:6 EH (MeV), 
W = 2:5 + 0-3 Ey (MeV), 
R = ae, where ro = 1°20 = 1:25 fm, 


e = 0-5 - 0:6 fm 
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Fig.3.8. Plot of total neutron cross-section against energy for 
C, Al and Cu. The full lines show optical-model fits. 
(Based on BOWEN et al. (1961). Nucl.Phys. 22, 640.) 


The success of the model is illustrated better 
by measurements of angular distributions as shown in 
Fig.3.9. The situation is paralleled in optics by 
diffraction through a circular aperture or, better 
still, the scattering (diffractive) of a small trans- 
parent drop of liquid of refractive index different 
from unity. Also, perhaps paradoxically, the model 
has had a great deal of success in the analysis of 
slow neutron data, where compound-nucleus (Chapter 6) 
resonances of width typically measured in. fractions 
of 1 eV are observed. Representing the nucleon- 
nucleus interaction by a simple central potential 
cannot possibly account for these manifestations of 
the detailed structure of the nucleus, but it has 
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Fig.3.9. Differential cross-section (expressed as a ratio to 
Rutherford) for elastic scattering of 30:3 MeV protons, 


showing optical-model fits. (Based on SATCHLER (1967). 
Nuel.Phys. A92, 273.) 


been shown that the optical model predicts average 
cross-sections where the range of averaging covers 
many compound nucleus resonances. For W very large 
(~ single-particle spacing) then the ratio of level 
width (without barrier) to spacing, which determines 
the average cross-section, will be of the order of 
1/27; but for W small this ratio, known as the 
"strength function’, will follow the shape of the 
optical-model resonance. 


PROBLEMS 


3.1. The lowest states in ae gPb are (energy in MeV in 


brackets) 4 (0), 3° (0-57), 5 (0-90), 4>"(1°63), 
and x (2:34). Interpret these in terms of the 
shell model. 

9 13 
Pe» 6 
- 1-17 nm, +0°07 nm, and -1°89 nm respectively. 


17 


Sues C, “20 have ground state magnetic moments 


Deduce values for (a) the spins and parities, 
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(b) the magnetic moments of these nuclei, as pre- 
dicted by the shell model. Comment on any dis- 


crepancies. 


Under what circumstances are (a) the single-par- 
ticle shell model, (b) the rotational model ex- 
pected to give a good description of the proper- 
ties of nuclei. Describe two features of nuclear 
behaviour for which either model makes charac- 
teristic predictions and show how these predic- 
tions can be tested experimentally. 


.The energy levels of the K=0 ground-state rota- 


tional band of 238) are given in Fig.3.6 (p.78). 


Estimate the moment of inertia I by equating the 
excitation energies with rotational energies. 

Why does I increase with increasing angular 
momentum? Why does the sequence contain only 
rigid — é AMR” 
for the appropriate rigid-sphere rotation and 


even values of J? Compare I with I 


comment on the difference. 


. Use Bose statistics to show that three i vibra- 


tional quanta can couple to form the states 
T + + + + 


re a, Mo A 

The optical model potential for neutrons on a 
medium-weight nucleus (A ~ 125) is given by the 
combinations (0,50,2),(30,40,8),(100,25,8), and 
(300,10,5), where the first number in a bracket 
is the neutron energy, the second “Vo> and the 
third “Wo all in MeV). Determine the mean free 
path at each energy near the centre of the nu- 
cleus; show that the nucleus is reasonably trans- 


ds 
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parent at each end of the range, and even at maxi- 
mum absorption there is ~ 10 per cent chance of 

a neutron passing diametrically through the nu- 
cleus. Explain the general trend of V and the 
mean free path. If the nucleus is transparent 
for near-thermal neutrons, how is it that they 
can be captured for long periods of time (see 
Chapter 6)? 


Work through the problem outlined on p.83 and 
verify the two equations given there. 


4. Spontaneous decay of Nuclei I: 
a and & decay 


In Chapter 1 it has been shown that the stable nuclei 
are surrounded by neighbours which are unstable to 
B-decay, whilst further away from the stability line 
nucleon emission can take place. These latter nuclei 
can only be observed in nuclear reactions. For mass 
values greater than 150, even the nuclei on the stab- 
ility line become unstable, not to nucleon emission 
but to emission of a-particles because the binding 
energy per nucleon of the o-particle is similar to 


that of a nucleus. 


a- DECAY 

Historically, a-particle decay was observed in 
the heavier nuclei (A > 206), and the lifetime for 
decay was seen to be very strongly correlated with the 
energy of the a-particle emitted (because the process 
is a break-up into two particles, the energy of the 
a-particle is well defined for unique parent and 
daughter states). For these nuclei all observed a- 
emissions lie within the energy limits 
4 MeV < es 9 MeV, the lower value being associated 
with a lifetime ~ 107!%s and the upper value with 
hor es according to the Geiger-Nuttall rule which can 
be expressed as 1n A = a+ b In £, where 1 is the 
mean life. This rule explains why, although almost 
all nuclei in the range A = 150 - 200 are unstable 
to a-emission, only a few have been observed to decay 
in this way; in general the available energy is very 
low and the lifetime too long to permit detection of 
the process. An example of an observed decay in this 
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‘24 Gd, with a lifetime of 10! 


an a-particle energy 2-24 MeV. These figures. are not 


region is 


years and 
compatible with the Geiger-Nuttall relationship as 
applied to the natural radioactive series, but it will 
be seen later that the connection between lifetime 
and energy depends very strongly on the nuclear charge 
ie 

Another historic experiment was the bombardment 
of the daughter nucleus with o-particles of the same 
energy as (or even greater than) those emitted by 
the parent. The scattering was observed to be of the 
pure Rutherford type, which, on the face of it, indica- 
ted that the nucleus was still behaving effectively 
as a point in its interaction with an a-particle 
approaching it from outside, even though this particle 
was closer to it than the region from whence (in the 
classical sense) an a-particle must have originated 
in the decay of the parent. This paradox was explained 
by Gamow in his theory of barrier penetration. 

In Fig.4.1 an idealized form is presented of the 


Fig.4.1. Idealized potential for s-wave a-particle plus product 


nucleus. Also shown is the potential step approximation. 
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potential to which an incoming or outgoing a-particle 
is exposed. The shape of the potential is dominated 
by the Coulomb interaction at long distances and by 
nuclear forces at short distance. These latter have 
been simplified :-to a square-well form, though they 
should conform to something like the optical potential, 
which will smooth out the sharp discontinuity and also, 
by means of the imaginary term, produce absorption 
within the nucleus. Also shown in Fig.4.1 is a cruder 
approximation to the potential in terms of potential 
steps. 

It is instructive to solve the wave-mechanical 
problem, in one dimension, of.a plane wave of wave- 
number K impinging from the left onto this latter 
potential, being partially transmitted, and finally 
emerging as a wave of wavenumber k in the far right 
region; in the intermediate region it must be represen- 
ted by real exponentials (multiplied by complex ampli- 
tudes) if, in this region, its total energy is less 
than the height of the potential step. This exercise 
is left to the reader. 

If in this problem the incident wave is of the 
form Texp(ikXr) and the transmitted wave Aexp(ikxr) then 
the transmission coefficient 7 is defined by the ratio 


of transmitted to incident flux 
Z 2 2 2 
r= vglal?/oplzl? = xlal?selzl?, 
where Vi» % are ingoing and outgoing velocities. For 
this particular problem 


a. 
To = —OekY  exp(-2yA), (4.1) 
(Kv +y") (ko +7") 
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where y = {2m(v-E) }4/2 


barrier, (R,-R). 


/n and A is the width of the 


There is a direct correspondence between the 
solution of a one-dimensional problem and the radial 
dependence of the solution of the corresponding 
three-dimensional problem with spherical symmetry. 
The one-dimensional wavefunction U(r) is simply con- 
nected to the radial dependence of the three-dimen- 
Sional wavefunction F(r)S(8,¢) by U(r) = rF(r), but 
the potential for the one-dimensional problem is not 
just V@) but V(r) + L(1+1)n?/ 2m 2, 1.e. allowance 
must be made for the centrifugal potential. Thus the 
solution above is appropriate for transmission of an 
~Z = 0 a-particle through the spherically symmetric 
potential step. Also the transmission coefficient 
is unchanged since the flux at r is v(r) .r?|F(r) |? in 
the three-dimensional case and therefore « v|u\? as 
before. It has been assumed that S(0,%) is correctly 
normalized, i.e. £13(8,6) |7an =]. 

However, considered as a three-dimensional pro- 
blem, there must be a node (in U(r)) at the origin. 
In principle, no simple solution as obtained can give 
such a node since, to conserve flux, the reflected 
wave on the left (i.e. inside the spherical potential) 
must have smaller intensity than the incident wave. 
But the difference (~% exp(-2yA)) may be very slight, 
and to good approximation can be ignored. The re- 
quirement of a node at the origin in addition to the 
(amplitude and) phase relationship between the two 
components derived from solution of the previous pro- 
blem overdetermines the solution within the potential, 
leading to a discrete set of solutions corresponding 
to energy levels of the a-particle. These solutions 
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are not quite stationary states since they are decay- 

ing via transmission through the barrier.' However, 

if exp(-2yA) is small then the decay is slow when meas- 

Po on a nuclear time scale. The decay constant Ay (in 
will be the outgoing flux when the ee rer in- 

side the 'nucleus' of radius R is normalized to unity: 

this normalization is realized by putting 


R 
f |Texp(ikr) + Jexp(-ikr)|*dr = l, 
0 


where Jexp(-ikr) is the reflected wave. Because of 

the node at the origin the integral becomes 
4)r|2sin?Ke dr giving Fale ~ 1/2R, since the point 
0 

R is very close to being an antinode. From the un- 

certainty principle, 


Tay more accurate solution of the problem intro- 


duces difficulties of normalization. The outgoing wave 
must slowly increase in amplitude as the distance 

from the nucleus increases, representing the fact that, 
at an earlier time, the probability of finding an a- 
particle within the nucleus was greater, and therefore 
the leakage through the barrier was greater. Thus 

the normalization integral diverges outside the nu- 
cleus. The approximate solution given can be looked 
upon as a perturbation approach in which an infinite 
confining potential has been switched off at t = 0 and 
the rate of increase of the final state (corresponding 
to a free a-particle) is determined when the initial 
state is known to represent one a-particle confined 
within the nucleus. 
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Ty = RA, = hvy lal? when |z|* = 1/22, 


a o} 
2 Dei 9 2 
Ok oe ig Ee Der, 4,2) 
o m K 2 m 2R Z 
eal 2mR 


This expression separates the a-decay width into two 
parts: the part n?/2mR?, which has dimensions of 
energy (width) and is the nuclear scale factor for 
width, and the dimensionless bracketed term, which is 
the transmission coefficient suitably modified for an 
effectively bound a-particle rather than a flux. This 
transmission factor is the product of a major term, 
the barrier penetrability, and a minor term which 
here represents the product of two mismatch factors 
arising from the two discontinuities in the potential 
(cf. reflection of optical waves at a glass surface 
or vibrations at a joint between stretched strings). 
In a more realistic calculation, introducing the 
Coulomb potential, it is reasonable to expect (and 


more detailed calculations confirm this) the pene- 
RR 
trability term to be replaced by exp{-2 fe v(e)dr}, 


since y is now a varying function of r. The limits 

of integration F and Ry are the points where y(r) = 0. 
The slowly varying y effectively changes adiabatically 
to an equally slowly varying k without discontinuity 
and so eliminates the mismatch factor at Ry: At R 

the discontinuity of the nuclear square-well poten- 
tial will introduce a mismatch factor, but a (rapidly 
changing) Saxon-Woods-type potential will reduce the 
mismatch though not necessarily eliminate it. This 
accounts for the fact that analysis in terms of a 
square well appears to require a potential of approxi- 
mately zero inside the nucleus whereas analysis in 


terms of the optical model usually gives a negative 
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internal potential (-25 MeV or so), though with con- 
siderable latitude since there are other variables in 
the analysis. Note that the modified barrier pene- 
trability may also be applied to the case 12 # 0, since 
Y(r) will be computed with the centrifugal barrier 
included. 

The chief source of uncertainty in the calcula- 
tion of decay widths and lifetimes is in deciding 
what fraction of the single particle width HOP Dan 
must be used for real nuclei. The original calcula- 
tions employed the full single-particle value, as 
given here, which for nuclei with A > 200 is about 
0.1 MeV, though the precise value will depend on 
the nuclear potential adopted. It was later suggested 
that 1 eV may be more realistic; this estimate was 
based upon widths of resonances in slow neutron re- 
actions (see Chapter 7). Such a large reduction in 
nuclear factor was off-set in the barrier-penetrability 
term by employing a larger nuclear radius. Nowadays, 
these radii so deduced are considered to be too large, 
which is to be expected since the levels from which 
the width was inferred are at an excitation of order 
7 MeV and statistical theories of nuclear levels re- 
quire that as level densities increase then level 
widths (for a particular channel, though not necessari- 
ly total widths) should diminish. Somewhere near 1 
per cent of the single-particle width is now con- 
sidered a reasonable estimate. Obviously, in face 
of such uncertainties, variations in mismatch factors 
due to details of the nuclear potential are of little 
Significance. 

An examination of the a-spectrum of a given 
radioactive nucleus reveals in general more than one 
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transition energy. The ground state of the parent 
nucleus can decay to a number of states in the daughter 
nucleus (of which the ground state may be one). In 
general too the higher the excitation of the daughter 
state the smaller the intensity of the transition, 
since the penetrability is so strongly dependent upon 
the a-energy. For the naturally occurring a-emitters 
a rough estimate is a factor of 3 decrease in transi- 
tion probability for 100 keV decrease in a-energy. 

In addition, the angular momentum requirements will 
contribute to intensity variations. Table 4.1 gives 
a rough guide to this effect. 


TABLE 4.1 


1 (ot 2 = & & 6 7 8 


Reduction 1 0-7 0:37 0-14 0:04 7x10° 


factor 


1x10° 1x10* 7x10° 


As an example consider the decay 75 4Pu > <u at OLe, 
Table 4.2 gives the measured transitions to states of 


known J (since 7 ePu has J = 0 in the ground state, 


the final J value is also the Z-value for the emitted 
a-particle). 

The first five transitions are to members of a 
rotational band. Such states have similar wavefunc- 
tions since they represent the same nuclear configura- 
tion in different degrees of rotation; it is therefore 
gratifying to see that the decrease in transition 
strength is largely accounted for by the variation of 
penetrability with energy and angular momentum. The 
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TABLE 4.2 
J ae al 4* 6* g* 1” 
(keV) 0 43 143 297 499 786 
Relative strength: 
Estimated 10°25 9x10°2 3x10? 3x1078 1x107 
Observed 1 0-39 1-2x107> 7x107> 1x1078 2x107 


sixth transition (others are given in the literature 
for this decay, but this one is typical) has a much 
weaker strength than would be accounted for by these 
means. It therefore reflects the nuclear contribu- 
tion to the process. Thus the conclusion to be drawn 
is that for members of a rotational band (or for that 
matter a vibrational band) the simple estimates of 
relative strengths should be fairly good, but in 
general the dominant factor is the nuclear intrinsic 
width when low-lying states are being compared. 
Finally, mention must be made of 'long-range' a- 
particles (this historical term is synonymous with 
"high-energy'). An excited state can, in principle, 
emit a-particles; indeed its decay probability into 
the a-channel will probably be greater than that of 
the ground state. However, it will be competing un- 
favourably with y-emission (see Chapter 5) at least 
for states near ground; at high enough excitation 
the reverse will occur but such states are not 
excited in the natural a-decay series. If the 


energy available to the a-particle is very large 


3 
7 
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then the decay probability will be correspondingly 
large. If also the y-transition probability is 
reduced because of high spin change, for example, then 
it may be possible that the excited state can have 

an observable a-decay. The key word is observable; 
modern techniques have revealed quite large numbers 
of 'long-range' transitions giving low-intensity a- 
groups above an intense ground-state group. From what 
has been said they are most likely to be found in con- 
junction with short-lifetime ground-state decays. An 
example is PO: with a ground-state half-life of 

304 ns and Be 8-79 MeV, having three eran 
tions, of which the most intense, 1:8 x 10 of the 
ground state strength, is from a state at an excita- 
tion of 1°8 MeV. 


B - DECAY 

In Chapter 1, 8-decay has been introduced as a 
process involving the simultaneous emission of a B- 
particle and a neutrino. The neutrino was postulated 
to satisfy the requirements of conservation of energy, 
momentum and angular momentum in this process. The 
B-spectrum is continuous up to an upper limit which is 
the energy expected for the 8-particle in the absence 
of an undetected particle. Similarly experiments in- 
volving the direction of recoil of the daughter nu- 
cleus indicate that momentum is not conserved - or 
preferably that an undetected particle has also been 
emitted. Also, to conserve angular momentum a 
neutron cannot decay into two spin-j particles, 
whether isolated or inside a nucleus. The neutrino 
must be neutral since charges balance without it. It 
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must have intrinsic spin 3} and a very small rest 

mass, since the B-energy at maximum accounts very 
closely for all the energy available. Efforts were 
made in the past to detect the neutrino, but failed. 
Now that recent efforts have been successful, it 

is not surprising that former methods failed since it 
is now known that neutrinos can pass through the Sun, 
the centre of which is a prolific source - indeed 
neutrino emission accounts for a considerable fraction 
of the energy emitted by stars. 

A simple theory of 8-decay was provided by Fermi, 
who exploited the similarity with emission of radia- 
tion. A proton (or neutron) changes into a neutron 
(or proton) by simultaneously creating the 8-v pair 
just as an excited state of an atom changes to a 
lower state by creating a photon. Note that in iso- 
i mass, but inside a nucleus a bound proton 
may be in a state of higher energy than an equivalent 
bound neutron, so the process can take place; if the 
total masses are favourable then the reaction can go. 

Considering first the decay of an isolated 
neutron (n>pte +v,) then the wavefunctions of the two 
light particles can be represented by plane waves in 
the directions of emission, Ve = Ne exp (i k,-1g) and 
vy) = N exp(ik.r.)-. The origin of the coordinate 
system is assumed to be the location of the nucleon 
both before and after the event (i.e. recoil has been 
neglected for simplicity). The plane-wave approxi- 
mation should be adequate for the neutrino, but 
for the 8-particle the long-range Coulomb interaction 
with the proton (or nucleus when at a later stage the 
nucleons are considered to be part of a nucleus) will 


SPONTANEOUS DECAY OF NUCLEI I 101 


distort the wave - this point will be taken up later. 
In the absence of guidance, Fermi assumed the simplest 
possible interaction namely that it is proportional to 
the simultaneous overlap of all four particles, i.e. 


He, = a [¥Rv,ae = 98 (0) VE (0) 


in this simple case of point nucleons at the origin. 
From perturbation theory (Fermi's golden rule) the 


transition probability per unit time is 


where dn/d#F is the density of final states. In 
applying this formula it is convenient to enumerate 
the final states when the system is confined to a 
specified (large) volume 2. The normalization factors 
for the plane waves, Ne and W, are then both 1/27, and 
the density of plane-wave states of a particle having 
momentum between p and (p+dp), with the particle any- 
2 apa /20eh?. The number of states for 
the 8-momentum in the range Pe to (pg*dp.) and the 


where in Q is p 
V-momentum in p,, to (p,,+dp,,) is then 

_ Z 2.3 2 a 

dn = (p gdp ,/2n h ).(pjdp/2a h~)Q™. 


In multiplying these two functions together it has 
been assumed that there is no correlation between the 
directions of emission of electron and neutrino. This 
is partly because of the initial simple assumptions, 
but chiefly because the final nucleon is available to 
take up the requisite momentum required for conserva- 


tion; because it is so massive it can do this and still 
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have little effect on the energy balance. It does 
not appear directly in the calculations, but its 
presence has been felt nevertheless. 

To determine dn/dFk it is necessary to transform 
from variables Pp» p,, to Pps E, using the equation 


2,1/2 


_ Z 2 
E = 6 ae + mec ) + E 


B? 
where £ is the total available energy. The trans- 
formation is 
dp 
ae ae -i 

dp gdp, = TE dp dz = Dy, dp dE = pe 
since the rest mass of the neutrino is either zero or 
close to it. Very close to cut-off of the 8-spectrum 
this neglect of the rest mass of the neutrino could 
distort the spectrum, but such effects have not been 
observed. Thus the density of states is given by 


222 22 
dn _ © Pep, 2 2°pi(E-Ea)* : 
= = P ’ 
dE an4n®, 8 ann ® 69 : 


using Eo = ep = E-Ee. 
The transition probability per unit time to pro- 
duce electrons in the momentum range Pe to (patdp.,) 


2 
_(g @ av2 
P(pg)4g = [S73] (F-Bg) PadPg, (4.3) 


or in terms of energy 


is 


dpe 
PRE Bg EUG? ide | 
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Note that the result is independent of 2, the volume 
in which the system was confined. In the non-relativi- 
stic approximation the energy spectrum is proportional 
to (E - are and this accounts for the general 
shape of a 8-spectrum. 

To extend this treatment to complex nuclei it- 
is necessary to consider the decay of each nucleon 
as an independent event and to sum the amplitudes of 
each event. So far the interaction is at a point, so 
the 8-v combination cannot take away orbital angular 
momentum; since no mention has been made of spin the 
combination must be assumed to be the spin=0 state. 
In the nucleus the decaying nucleon merely changes 
its charge but remains in the same nuclear state - 
the instantaneous effect can be represented by 
2 lies where a is the isospin raising or lowering 
operator acting on the kth nucleon in the nuclear 
wavefunction y, (t for 6 emission, t for a a 
This instantaneous wavefunction must be projected 
onto the final wavefunction to determine the nuclear 


contribution to the process, giving the matrix element 


M 


+ 
fete 


The summation may be taken over all nucleons since T. 
acting on a proton gives zero - it changes neutrons to 
protons - and correspondingly for t. Obviously, from 
what has been said above, the selection rule AJ=0, with 
no parity change, holds for the nuclear states. 
Extending the theory to include the possibility 
of the 8-v combination taking away one unit of spin 
leads to a parallel decay obeying the selection rule 
AJ=0, +1, with no parity change (but J=070 excluded). 
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The nuclear matrix element must contain a spin opera- 
tor capable of changing the spin orientation of the 
nucleon kK - it must be in the form of a vector in order 
to take away (spin) angular momentum. There is no need 
to go into its form here; the simple picture of the 
B-v combination being emitted in a wave having L=0, 
S=] is sufficient to arrive at the above selection 
rule. The spinless emission is known as the Fermi pro- 
cess and the spin-1 emission is ascribed to Gamow- 
Teller. The processes so far discussed are termed 
"allowed'. 

The treatment of nucleons as point particles is an 
approximation which allowed the 8-v combination to 
take away no orbital angular momentum. Whilst there 
is little to be gained by giving the nucleon structure 
in the decay of an isolated neutron (to preserve 
parity a two-unit change would be necessary and this 
is unlikely for such a small structure), it should be 
noted that even for point nucleons it may be possible 
to change the orbital angular momentum of a nucleus. 
If the emitting nucleon is at the periphery of the 
nucleus then, classically, its recoil momentum will 
produce a change of angular momentum. The greatest 
change will occur when an electron of maximum energy 
is emitted tangentially from the surface, and will be 
of order (E2/e) Ry. Taking Ee ~ 1 MeV (for which 
Bele = Pe is only approximate) and Ry % 6 fm (a medium- 
Size nucleus) gives 0-03h for the classical angular 
momentum in this extreme situation. Thus the taking 
away of one unit of angular momentum must be looked 
upon as possible but improbable. Mathematically, the 
two plane waves of 8 and v are amalgamated to a single 


plane wave exp(ik.r), where k = ke + ky, and, since 
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the interaction is still assumed to be at a point, 
fg = r, = ©. This plane wave is expanded within the 
confines of the nucleus giving 

exp(ik.r) = 1+ ik.r - }(k.r)* ... = . The matrix 


element for the process will now become 
+ 
= * 2: | 


where for convenience the possible spin factor is 
omitted. If Ve and v5 have opposite parities then 

the leading term in the bracket will integrate to 
zero, but the next term, which is itself odd, may not 
do so. If polar axes are chosen such that k lies 
along z then k.r, becomes kr) cos 6. The angular wave- 
functions of nucleon & in Ve and v5 are in the form 
of spherical harmonics Y, toro) » and their properties 


are such that [YF (6,¢)cos 6 Y, o (6,¢)dQ will be 
| f 14 


om 
zero unless Me = mt and Le = Ls + 1. Thus the kth 
nucleon changes angular momentum by one unit (cf. 
emission of a photon by an electron in an atom) and, 
by recoupling to the rest of the nucleus, the selec- 
tion rule follows: AJ=0, +1, with change of parity 
(J=0 + 0 excluded) for the Fermi process. By coupling 
up to S=1, the corresponding Gamow-Teller selection 
rule is Ad=0, +1, +2, with change of parity. 
Consideration of the radial integral indicates 
that these transitions will be down by a factor (kRy)* 


2 ; 
term arises from 


on the allowed transition - the k 
the integral and the R- 


2 2 
fyery,r dr must be % Ry lel v, |r dr, and we assume 


from the fact that an integral 


that all integrals of the latter type have approximate- 
ly the same value. These types of transition are 


known as 'first forbidden' since they are down by the 


Ty different choice of axes results in Me = ms come Ie 
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factor (KRy)”, which, as previously estimated, is 

~v io or 107. The next term in the series would 
produce second forbidden transitions and will be down 
by the same factor again. 

Before examining whether this classification con- 
forms with experimental fact, it is necessary to con- 
sider the effect of nuclear charge on the process 
Since this will vary from nucleus to nucleus. The 
wavefunction of a particle in a Coulomb field is con- 
plicated, but here we merely need to know the behaviour 
at r = 0, where the process occurs, and r = ~, where 
the products are observed. If at infinity the wave- 
function for the 8 is the original plane wave then 
lv. (0) |? will be muleTPi ied by the factor 
mtr, E) = 2mn{l-exp(- 2mn)} 1 , where n = +e 7/4m€ hv, 

(+ for electrons, - for positrons), Ve being the speed 
at infinity and Z the atomic number of the final nu- 
cleus. The Coulomb correction enhances the probability 
of electron-emission and decreases that of positron 
emission, especially at low energies. In the latter 
case it is effectively the Gamow barrier-penetrability 
factor since the positron may well be created with 
negative kinetic energy and need to penetrate the 
barrier before emerging with positive energy. The 
enhancement of the 8 -process occurs because the 
attractive potential increases |¥g(0) | 4 relative to 

| vg (~) | 4. The normalization of the wave is largely 
determined by |vg(0) |? since 2 is very much larger than 
the nuclear volume, but the transition probability 
depends on |vg (0) |. The effect of the Coulomb 
potential on §*-spectra is shown in Fig.4.2. 

With the Coulomb correction the momentum spectrum 
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Fig.4.2. Effect of Coulomb potential on 8-spectra for cer 


decay, which shows g and B° (and electron-capture) 
decay. Note that the 8 energy spectrum has finite 
intensity at zero kinetic energy, but the momentum 
spectrum does not (the effect of the relationship 

dz = vdp). (Based on EVANS (1955). The atomte nucleus, 
McGraw-Hill, New York.) 


for allowed transitions is of the form 


= _ 2.2 
P(pa)dp,g = CF(Z,E,)(E-E,) padp,, (4.5) 
where the constant C contains terms independent of the 
B-energy. Thus a plot of {P(p,)/Epat against Ee 


is a straight line intersecting the abscissa at E£. 
This, or its equivalent for the energy spectrum, is 
known as a Kurie plot, an example of which is shown 
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Fig.4.3. 
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(a) Kurie plot for allowed decay of tritium. Inset 
shows fit near end point. Lack of deviation indicates 
that the rest mass of the neutrino must be small. 
(Based on LEWIS (1970). Nucl.Phys. A151, 120.) 

(b) Kurie plot for first forbidden decay of 895. The 
curved plot is the usual plot with Fo in the ordinate 
function; the linear plot contains a theoretical cor- 
rection for the first forbidden transition. (Based on 
WOHN and TALBERT (1970). Wucl.Phys. A146, 33.) 


in Fig.4.3. The linear plot results for allowed tran- 


Sitions only; in forbidden transitions terms in 


k 


Z 


aed 
= ke 


+ Ke + 2k Ky appear, which depend upon the 


B 


relative direction of the 8-v. On integrating over 


all directions, the term k,.k, averages to zero, and 


B 


the two remaining terms inject a factor {pe+(E-Eg)"/e°} 


into the momentum spectrum for a first forbidden 
transition. The Kurie plot will not be linear, but 
a modified Kurie plot, including this factor, can 
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restore the straight line. This technique can be used 
to determine the type of decay (see Fig.4.3). 

A convenient classification of 8-transitions is 
based on the lifetime of the decay 


1 


2 2 
= [P(pa)dp.p Sg |Me, | 


f; 

which defines f as an integration over all the energy 
dependent terms in the spectrum and contains all the 
fundamental constants excluding the coupling constant. 
Thus for an allowed transition 


3 Se. Lac 1 max 2 2 
= 1 h B-E BE 
fo (2 Cc ) * - ( 3) paF(Z, g) IP 2 - 


The product fobs: where ty is the half-life (=tT ln 2), 
is therefore a measure of Cae ke for an allowed 
transition and is found to vary between 10° and 10° 

in magnitude. Transitions having values near 10° are 
known as 'super-allowed' since they correspond to 
(almost) complete overlap of the nuclear wavefunctions. 
They are usually Fermi transitions between states of 
the same isospin multiplet (e.g. B-decay between mirror 
nuclei), but there exist a few Gamow-Teller transi- 
tions in very light nuclei where the overlap is also 
complete. The remainder of the range corresponds 

to a diminished overlap. The same product fot for 

a first forbidden transition will be of the form 

{97 (1M; ¢|?/Rq) (By) 2} >, where the bar denotes a suit- 
able average over the 8-spectrum. For complete over- 
yD in value on an 
allowed transition giving a range of order 10° - 10 


lap this is up by a factor (kR 


Similarly, a second forbidden transition will give 


fot values > 1019, Thus there is no strict division 
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between the different classes, but for certain ranges 
of fot values it is possible to deduce the type with 
a reasonable degree of confidence. 

Note that these last considerations should really 
apply only to 8 transitions since, as has been men- 
tioned in Chapter 1, an alternative decay mode to 8° 
is available. This will be discussed briefly in the 


next section. 


ELECTRON CAPTURE 

Since atomic nuclei are in general surrounded by 
atomic electrons an alternative process to B*-decay 
can take place, namely, 


p+p 7 n + Ve (within the nucleus). 


This process is of greatest importance at low avail- 
able energies - indeed there is a range of available 
energy (see Chapter 1) of 0-1-02 MeV over which ce 
decay is impossible and e -capture only can take place. 
Also it is expected to be largest for K-electrons (pro- 
vided the energy balance is not critical), since they 
have strongest overlap with the nucleus, and for nu- 
clei of large Z, since the electron orbits are more 
tightly bound, again leading to greater overlap. 
Applying Fermi's golden rule to this process, 
the density-of-states function will be that of the 
neutrino alone - the nucleus is constrained to recoil 
in the opposite direction and both particles have 
fixed energy. The probability of finding the electron 
at the origin of the coordinate system is now no longer 
inversely proportional to the volume of the enclosure 
but to the volume of the K-orbit (for K-capture), so 
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once again 2 cancels (it has also disappeared from the 
density-of-states function for the electron which is 
now unity - a defined state). If F is used, as be- 
fore, to denote the maximum available energy for the 
B*-decay, then the energy available to the neutrino 
will be (F + 1:02 - By) MeV, where By 


energy of the K-electron. There is now no integral 


is the binding 


to perform since no variation is allowed of the pro- 


ducts, and 


2 
Cease 
2nch-o> 


27 2 (4.6) 


ie 
tT HY 


1 2 3 (E+1°02-B 
K 


Mel” = (% 
where ais the (hydrogen) Bohr radius, so that a/Z 

is the mean radius of the K-electron wavefunction in 
the parent nucleus, and where the matrix element Me 
is the same as before, and 1/t, is the probability per 
second of K-capture. Thus the probability of the pro- 
cess increases as z>, whilst that for g” -decay de - 
creases with Z (due to the Coulomb effect); for light 
nuclei B° -decay will dominate except for an energy 
range in FE from the negative value (-1°02 + Bx) MeV 

to slightly above zero. But for high Z, K-capture is 
overwhelmingly predominant for B*-energies commonly 
encountered. 


RELATIVISTIC THEORY AND PARITY VIOLATION 

The account of 8-decay so far described has been 
essentially non-relativistic (NR) in its approach, 
apart from the use of a relativistic equation connect- 
ing momentum and energy - a trivial aspect. For re- 
lativistic invariance it is essential that time and 
space be treated on an equal footing, and this is not 
the case with the Schrodinger wave equation, which is 
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of second order in spatial derivatives and first 
order in the time derivative. Putting F = ih 3/dt, 
p. = (h/i) d/dx into the relativistic (R) equation 
pe = ety + me of gives the Klein-Gordon equation which 
is second order in all four derivatives. However, 
its solutions are found to conform to Bose statistics, 
e.g. photons, and the equation cannot be used for 
spin-} particles. Dirac looked into the conditions 
under which this second-order equation factors out 
into two first-order equations and further examined 
the solutions of the resulting equation. He found 
that the solutions were four-component vectors in a 
space which was closely allied to the two component 
space of a NR spin-3 particle. Indeed two of the 
components are down by a factor v/e on the other two 
such that in the limit v + 0, the R system reduces to 
the NR system. In addition the 'small' components 
have opposite parity to the 'large' components, crea- 
ting the possibility of parity violation at relativi- 
stic energies. 

In 8-decay, the approximation v/e < 1 is rarely 
valid for the electron and never so for the neutrino 
- even for tritium decay the 8-particle has maximum 
kinetic energy of 18-6 keV corresponding to v/e v 0°27. 
Thus the product of Vev,)> which in the NR limit gives 
four terms which can be regrouped into a scalar (S=0) 
and a three-vector (S=1) gives in the R limit six- 
teen terms which can be regrouped into two scalars, 
two four-vectors, and a six-component antisymmetric 
tensor. The two vectors are respectively polar and 
axial; the former changes sign under the parity opera- 
tion whilst the latter does not, e.g. r and p are 
polar three-vectors whilst r A p is an axial 
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three-vector. Of the two scalars one is a pseudo- 
scalar in that it changes sign under the parity opera- 
tion; an example of such a scalar is the triple three- 
vector product a.(bsc), where a, b, c are polar three- 
vectors. Such a product can be used classically 

(e.g. to represent a volume) but the sign property 

is usually ignored. 

Thus a full relativistic treatment is much more 
complicated than the picture presented here and has 
built into it terms which mix the parity of the re- 
lativistic components. It is possible to choose an 
interaction in such a way as to eliminate the parity- 
violating component, and this was believed to be the 
case until specific measurements were made to test 
parity conservation in §-decay. Such measurements 
are designed to determine whether pseudo-scalars are 
necessary in the description of the reaction. Such 
a scalar product is o.p, where o represents the NR 
spin of the electron and p its momentum or J.p, where 
J is the spin of the initial nucleus and the §-momen- 
tum p is now measured relative to the same axis as J. 
Obviously, if parity is conserved then(o.p) and 
(J.p) must be zero. The first measurements were made 
on the decay of 800, with the initial state polarized 
by cooling in a magnetic field. It was then found 
that the number of 8-particles directed into the hemi- 
sphere with a component along the field was different 
from that directed into the opposite hemisphere, in- 
dicating that (J.p) was not zero. The final conclusion 
to be drawn from such measurements was surprising, 
not so much because parity was not conserved but be- 
cause its violation appeared to be the maximum pos- 
sible when extrapolated to the limit v > ec. 
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When polarization measurements are ignored as in 
the simple determination of the spectrum, then the 
Simple predictions of the Fermi theory are largely 


correct, but, in more detail, separation into 'large' 
and 'small' terms of the nucleon wave functions is nec- 


essary. Part of a first forbidden transition can arise 
from the 'small' term of an allowed transition, and vice 
versa, though in the latter case the effect will not be 
discernible since it is also reduced by the factor (kRy)? 
relative to the main component; in the former case the 
factor serves to enhance this extra component. Thus, 
whilst under most circumstances the simple theory copes 
with allowed transitions, great caution is needed in 
extending its results to forbidden transitions. 


PROBLEMS 


4.1. Solve the one-dimensional problem of transmission 
through the rectangular potential as given in 
Fig.4.1 (p.91). Show that in the limit of the 
barrier width approaching zero, the transmission 
becomes that calculated for a potential step. 


r 
4.2. Solve the integral f 0 ydr (see p.95) for the 
r 


more realistic sGtentiall of Fig.4.1., where the 
limits of integration are-.the classical distance 
of closest approach and the nuclear radius (use 
the transformation r = a cos “8). Show that, 

in the limit coo barrier height, the Gamow 
factor exp (-20Z2e*/4ne hv) is the important term 


in s-wave barrier penetration. 


4.3. Show that the non-relativistic formula for the 


B-spectrum for an allowed transition gives a 


4% 


ory 
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mean kinetic energy of the 8 of z the maximum 
energy, whilst the relativistic formula, for 


a very high maximum energy, gives a factor }. 


Derive expressions for the maximum recoil energy 


of the nucleus of mass A when the 8 is emitted 
at the mean energy in both cases. (Neglect Cou- 
lomb effects.) 


.Derive an expression for the electron momentum 


spectrum in allowed 8-decay (ignore Coulomb cor- 
rections). Find the approximate dependence of 
the total decay rate on the maximum electron 
momentum Po when Po > mec. 

Compute the value of the partial decay rate 
for the pion decay T > ne +e + Mesh 4-5 MeV 
given that the Fermi decay of 145 to ay (excited) 
has a momentum end point Py = 2°26 MeV/e and a 


mean lifetime of 103 s. 


Explain the meanings of the following terms: (a) 
first forbidden, allowed and super-allowed 8- 
transitions, and (b) Fermi and Gamow-Teller matrix 
elements. Classify with regard to (a) and (b) the 
following: (1) n> p; (ii) He (0°) > le 


(Gi) "000" NCO Vg. Poe = Saoee 10°: 

nee Sees. eae eT 7 

(iv) P2s(5") > 730105"), fot = 1-0 « 10°; 
S6scetee a SOs 

(v) 1701 (2 ) > ight (0 aes 


;Be decays by electron capture chiefly to the 


ground state of ably and the recoil energy of 

the Li has been measured to be 57 eV to 1 per 
cent accuracy. If the mass difference between 
the two atoms is 0-862 MeV/c? show that the re- 
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sult iS consistent with a neutrino rest mass 


anywhere in the range of zero to % 10 keV/e*. 


5 . Spontaneous decay of nuclei Il: 
electromagnetic transitions 


Thus far the discussion has been restricted to the 
decay of nuclear ground states, apart from the oc- 
casional intrusion of an excited state decaying, for 
example, as an isometric state, or giving a 'long- 
range’ a-particle. In the next chapter states at 

high enough excitation to emit nucleons will be con- 
Sidered, but here the electromagnetic decays of states 
within 1 MeV or so above the ground state are treated. 
The electromagnetic process is, in general, a stronger 
one than the weak 8-process, so that such excited 
states only rarely 8-decay, though there are a few 
instances when the excited states lie close to the 
ground. 

The treatment of radiation presented here is 
largely based upon the. corresponding process of elec- 
tric dipole radiation from excited atomic states. 
Later in the chapter it will be shown that such radia- 
tion can take place only between states satisfying the 
criteria AJ = 0, + 1 and change of parity. But, ex- 
perimentally, using techniques to be described briefly 
in the next chapter, it is found that in nuclei y- 
transitions frequently take place between states not 
conforming to these criteria. It is therefore neces- 
sary to extend the treatment tO cover the case of 
emission via the higher electric multipoles and via 
magnetic multipoles. Such transitions can be ignored 
in atoms since they occur so seldomly and then only 
under very special conditions. If the rapid process 


of electric dipole (E1) radiation cannot take place, 
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then the state usually de-excites by collision. In 
nuclei, on the other hand, it transpires that El 
emission is not a great deal. faster than the other 
processes, and in any case most low-lying states 
have the same parity as each other since they pro- 
bably arise from different couplings of the same 
shell model particles or from collective motion. 

But nuclei are isolated from each other by vir- 

tue of their charge, and are coupled only weakly 


to the atomic electrons; a state once formed can only 
decay via the electromagnetic process, though it need 
not necessarily lead to y-emission since the Coulomb 
coupling to the atomic electrons can result in the 
emission of an electron termed, misleadingly, a 
"conversion electron'. This latter process is electro- 
magnetic and gives the same information concerning 
nuclear states as does y-emission. 


y -RADIATION 

The difficulty in calculating the spontaneous 
decay of an excited state is that one does not know 
the perturbing term in the wave equation which is 
responsible for the transition. However, the reverse 
process of absorption can be treated by using the 
electric vector of the radiation as a time-dependent 
perturbation. Einstein, using a statistical argument 
(known as ‘Einstein's As and Bs'), obtained a re- 
lationship between the two processes. A fuller ac- 
count is given in Pauling and Wilson (1935), but 
very briefly the argument is as follows. 

The coefficient of absorption is defined such 
that the probability of a transition from a lower 


state n to a higher state m is given by B p(v_)> 


+m 
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where P(v dv is the energy density of radiation in 
the range Vee CO dv) with ae the frequency of 
= (BE )/h. The probability 
of transition down is given by A + BoP, Jes 


where As is the coefficient of spontaneous emission 


n 
the transition, i.e. v 


and Bas, is the coefficient of induced emission. If 
we consider the equilibrium of this two-state system 
in a bath of radiation at infinite temperature then 
ete Te es since the two states will be present in 
equal proportions, and the induced processes will be 
much more probable than the spontaneous process. It 
is assumed that the two states are non-degenerate; 
this can be achieved by applying a magnetic field to 
separate out sub-states of » and m - the transitions 
between a single sub-state in each level is then 
under consideration. If now the equilibrium is con- 
sidered at a finite temperature, then the ratio of 
abundances will be given by the Boltzmann factor, 
N/M exp{-(£_-E)/kT}, whilst p(v) is given by 
Planck's radiation law, 


o(v) = (81h v/e") Lexp(hv/kT)-1} 1 
But at equilibrium 
7 -] 
Nt Na 7 ore ¥ B aon Og asm PI gt 


Ts 
if) 


3 3 
ae (8Thv_ /e 1 ee 


3 3 
(Zhu /te ee 
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It is not proposed to go through the complete 


calculation of B as given in the above reference 


x 
for El transitions. but to show how to modify this 
approach in order to produce the other multipoles. 
The perturbation resulting in absorption was taken 
to be due to the electric vector of the radiation 
field, Ep exp{i(k.r-wt)}. Since k = mv/h = E/he, for 
a photon of ~ 1 MeV and a nuclear radius of 5 fm, 
then kRy ~" i: Thus the simplification Ey exp (-iwt) 
is a fairly good one, though not so good as in the 
atomic case. Now a uniform electric field is given 
by a dipole at infinity, so this approximation, giving 
a perturbing potential - eEy-r exp(-twt), leads to El 
absorption of amplitude proportional to the matrix 
element LV ez,)p dt, where the summation 7 is taken 
over the protons of the nucleus and where Ey has been 
taken as directed along the z-axis. This gives 

‘nn a (a 
3he 9he 


p &®) 2 


mn 


A D 


mn 


m>n ’ (5.1) 
where in the latter expression averaging over all 
directions of Eo has been made in forming the dipole 
matrix element - and also averaging over all sub- 
states m and , which is correct for the initial 
state but not for the final state. A further factor 
(27, + 1) should therefore be present, but it is a 
degree of refinement which will be ignored. The suit- 
ably averaged dipole matrix element becomes e times 
the radial integral contribution to the moment for a 
particular choice of sub-states. 

The higher electric multipole contributions can 
be obtained by using the expansion 
exp(ik.r) = 1 + ik.r - (k.x)? ate in the expression 
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for the radiation field. The leading term has given 
a component of the El matrix element; the next term 
gives a component of the E2 matrix element 
funkeQ x.2,)W dt if the propagation vector k is taken 
along the x-axis (it must be perpendicular to Eo? 
since electromagnetic waves are transverse). And 
Similarly for higher orders. 

From the matrix elements the selection rules 
arise. Since z is an odd function, the product 
vey must also be odd to give a non-zero dipole tern, 
i.e. there must be a change of parity. Similarly, «xz 
is even on reflection through the origin, so the E2 
transition is characterized by no change of parity. 
In more detail, if the particle making the transition 
is in a state of orbital angular momentum (7,7, ) 
before and (Z,m,) after the transition, then the 
angular dependence of the above component of the El 
matrix element will be a eee ree eae 


and the properties of the spherical harmonics (Pauling 
and Wilson 1935) are such that this integral is zero 
unless Me = Ms and Lent; = + 1] (there are two other 
components which lead to Me - m, = + 1). Similarly 
the above component of the E2 matrix element, 


Jy% (8,¢)cos @sinOcos oY (8,¢)dQ, will be zero 
Em bams 
unless Me - m;, = 4 1 and if - Ls = 0, + 2, with 0 > 0 


excluded since it cannot satisfy the former condition. 

Again, there are other components, for which 

Me ~ m; = 0, +2, e.g. the xy component giving Am = #2. 
For nuclei consisting of a single proton out- 

side closed shells, the selection rules on (J,M) will 

be the same as those obtained on (Z,m) for transitions 

between states corresponding to excitations of this 
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proton, apart from the addition of a spin vector. In 
the general case it is necessary to decouple the rest 
of the nucleus from the particle making the transi- 
tion and then to recouple after the transition. The 
mathematics of these coupling operations is beyond the 
scope of this book, but it is not difficult to see 


that the selection rules will become 


El AJ = 0, + 1 (0 + 0 excluded), with parity 
change, 
E2 AJ = 0, t+ 1, + 2 (0 0,4 > 3} excluded), 


with no parity change. 


The exclusions arise because all components of a 
given moment must be allowed since they transform 
into each other on rotation of axes. Thus Am = 2 
(and therefore AM = 2) must be possible in E2, and it 
is patently not so for the excluded transitions. 

So far no mention has been made of magnetic 
multipole transitions for the very good reason that 
the interaction has been taken to be that of an elec- 
tric field. But electromagnetic radiation also pos- 
sesses a magnetic vector Ho» 
E 


0 and k. Again the leading term will arise by 
assuming the (magnetic) field constant over the nu- 


perpendicular to both 


cleus giving rise to the perturbation term 
“y-By) exp(-iwt), where the magnetic moment 

» t t : y 
u = By 2(g9,8; 1 972.) 5 where 8, is the nuclear mag 
neton and 95°91 are the appropriate spin, orbital 
g-factors for the nucleon z. To simplify matters, 
consider the case of a single neutron outside a 
closed-shell core then 

oe -1 

Uekg = 9 9 byl eg B yt ty Boyt 2 208)? MneTS 3? ey? “2 
are spin operators (see Appendix D). If the nucleus 
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conforms to the J/-J coupling scheme (and most do, at 
least approximately, and especially near closed shells) 
then the two states J = 2 + 4 are distinct and ortho- 
gonal. They can be expressed in terms of orbital and 
Spin product wavefunctions thus: 


|2+3,m) o|Z,m+3>]3,-4) ” B|Z2,m-4 >|4,+3), 


|2-3,m) = -B[Z,m+}) 3,-4) + af2,m-4)|4,4)), 
where a and 8 are normalized such that la] 7+] B17 = ] 
(they are known, but the values are not needed here). 
Since és. = LAG, an it follows that S| t+} ,m) contains 
both |Z+4,m) and |Z-41,m) as also does 8_|t-3,m)>. In 
treating oe and a it is a es to recombine them 

to form st = Sie and § = og with weighting 
factors which will depend upon B and B. . Now, 


- Ox Oy 
from Appendix D, s |Z+4,m) « |Z,m+} |4,+4) 


» which is 
a linear combination of |Z+3,m+1) and |2-4,m+1); and 
similarly for s . Thus the effect of s.B) is to mix 
the original state |j,m) into states belonging to 
both couplings of zZ and s and also to mix the orien- 
tations, Am = 0, + 1. Thus the matrix element will 
be non-zero if the final state contains one or other 
of these components. A similar argument applies for 
the term in 7.B) for the proton. The selection rule 


for magnetic dipole single-particle transitions is 


M1 AL = 0, Am = 0, + 1, [Ag] = 0, 1, no change 
of parity, 
where |Aj| = 1 allows the alternative coupling scheme. 


In the general case it is necessary to decouple and 
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recouple with the rest of the nucleus giving 
Ml : Ad = 0, + 1, no change of parity, (0-0 excluded). 


The next step would be to expand the exponential 
again to give higher magnetic multipoles, but we shall 
leave it at this point, except for observing that the 
parity will alternate with the higher magnetic multi- 
poles, as for the electric. It is interesting to 
note that, whilst El transitions are essentially from 
one shell to another, Ml and E2 transitions can be 
made between states with the same basic shell struc- 
ture. In fact most low-lying nuclear states are 
different configurations of the same basic states so 
one encounters Ml and E2 transitions much more often 
than El. 

From the foregoing it is possible to make crude 
estimates of the ratios of these three main types of 
transition. E2 is seen to be down by a factor 
(kB)? ~ 1/2000 on El for a medium-weight nucleus, 
whilst 


1/2 


(M1/E1) av hg. ,Bolvts Vdt/ek sSysay dt 


% (8\9,/2eRy) (Bp/EQ) 


(h /MeRy) (g,/4) 


From the uncertainty principle h/Ry ~ Mv leading to 

a ratio of intensities v v for which, for a typical 
kinetic energy inside the nucleus of order 50 MeV, 

takes the value %v ww . This is perhaps an overestimate, 
a factor %v — is closer. 


Blatt and Weisskopf (1952) make estimates of the 
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matrix elements under simple assumptions arriving at 
ary for El and eke for E2. For M1 they take 

M1/El %~ 10(h/MeR,)*, where the factor 10 allows for 
the large-spin g-factors of proton and neutron. In- 
serting these estimates and the values of the phy- 


sical constants gives 


3,2/3 
r = 0:07E 
(EL) 0°O7E A 
-8 5 ,4/3 
T = 4- E 
~(B2) = 4-9 x 10 yA 
oo e 3 
r (M1) = 0 0212, ; (5.2) 


where iy is in eV, Ey in MeV, and A is the mass number. 
The radiation width ry is connected with the decay 
half-life by 4720S) ~~ 4°6 x 1070s (see also Fig. 
Beds) 


105 


Half-life 11/s 
= 


a 
Oo 
l 


1 Q -15 
0-01 002 005 O1 02 OS 1 2 5 
y-ray energy /MeV 
Fig.5.1. Lifetime-energy relationships for y-radiation, based on 


Weisskopf's single-particle estimates. 
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Thus it is seen that the single-particle estima- 
tes lead one to expect El transitions to be stronger 
than Ml or E2 at ~ 1 MeV. But, as remarked previous- 
ly transitions between single-particle states will 
rarely be found at such low excitation. In addition, 
there appears to be a tendency for neutrons and 
protons to move together in nuclei, and, from the 
remarks below concerning centre-of-mass effects, this 
will very much reduce the effective charge of the 
radiating system for El transitions. Except for a few 
transitions in the light elements low-energy El transi- 
tions are many orders of magnitude weaker than single 
particle. Since Ml transitions are expected between 
similar components in complex wavefunctions, they are 
rarely found to approach single-particle strength, 
and are typically a few per cent of this value. On 
the other hand, E2 transitions are often found to be a 
factor of 100 stronger than the single-particle esti- 
mate. The matrix elements so far discussed have been 
treated implicitly in terms of the shell model in which 
only a small number (% 1) of particles can find suit- 
able initial and final configurations. If, however, 
all Z protons move coherently, then it is possible 
in principle to obtain transitions up to 2° stronger 
than single particle. In a rotational band the states 
have such a correlated motion in which at least 20 
nucleons are taking part. 

The switch from protons to nucleons is deliberate; 
it is found that neutron transitions are as strong 
as proton transitions for both El and E2, though for 
different reasons. In the former case it is because 
of the finite mass of the nucleus; when a proton moves 
outward, the rest of the nucleus must move in the 
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opposite direction to keep the centre of mass unchanged. 
By simply taking first moments (in the classical sense) 
about the centre of mass it is seen that the system 
behaves as though it were infinitely massive provided 
that the charge on the proton is taken as 

te(1 - Z/A) ~ + 3e. A similar argument gives the ef- 
fective charge on the neutron as -eZ/A ~ -je. Since 
first moments were taken this applies to El transitions. 
Performing the same operations with second moments 

gives only small departures from +e and zero for the 

two nucleons. Thus the centre-of-mass effect does not 
produce large E2 neutron transitions. That E2 neutron 
transitions are as large as proton transitions even 

179 ana Lf 
neutron or proton) must be ascribed to a collective 


in nuclei like F (closed shell plus one 
polarization of the core, and in both cases it must 
be a large effect. 

Since E2 are in general larger than single par- 
ticle and Ml in general smaller, they are found to be 
comparable in strength. The commonest transitions in 
nuclei are therefore found to be Ml or E2 or mixed 
M1/E2 transitions; in the latter case either type can 
be larger though with a tendency towards E2 in the 
heavier, distorted nuclei and to Ml in the lighter 


nuclei or nuclei near to closed shells. 


INTERNAL CONVERSION 

It is possible for an excited nucleus in an atom 
to de-excite by emitting an atomic electron by an 
electromagnetic process known as internal conversion. 
It has an atomic counterpart known as the Auger effect 
when, as an alternative to X-radiation, an outer-shell 


electron is emitted instead. This process could be 
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looked upon merely as the photoelectric effect in 
which the atom happens to be the same one as contains 
the radiating nucleus; indeed, such a process must 
occur and is more likely to occur for this particular 
atom than any other atom, since the radiation density 
is vastly greater than for any other atom. However, 
the photoelectric effect is not very probable for ‘any 
other atom', so even for the particular atom it is 
likely to account only for a very small fraction of 
decays. But internal conversion can account for the 
major fraction of low-energy transitions in heavy 
nuclei, so a single-stage process is indicated in 
which the energy is transferred directly from the 
nucleus to the ejected electron. The Coulomb field 
provides the perturbation for the process. 

The Coulomb potential at a given point pa in 
space due to the nucleus is of the form 


_ 1 
Vin eee “Gace 
ee eg ey AG 

where the summation is taken over the protons, of 
position vector r., inside the nucleus. If the field 
point r_ lies outside the nucleus, i.e. |r_|>|r.| 

—e —e —21 
for any r, inside the nucleus, then the term ye ee ey 
can be expanded in powers of r..r /r? . Whenr.-is 

—t*-e' e t 

averaged over the nucleus each term in the expansion 
produces a nuclear electric multipole moment just as in 
the corresponding expansion of the radiation field. There 
is also a similar expansion in terms of magnetic multi- 
poles arising from the magnetic interaction between the 
nuclear magnetic moment and that of the.electron. Thus 
Y-emission and internal conversion depend upon the 
same nuclear property. Without going into details, 
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the process is expected to be more probable for the 
heavier nuclei, since the electron being ejected is 
then on average, closer to the nucleus and subject to 
a larger electromagnetic field. Also from considera- 
tions of density of final states, the process is ex- 
pected to be larger, relative to y-emission, for lower 
transition energies, provided they exceed the binding 
energy of the electron. The heavier the particle 
the greater its momentum for a given energy (assuming 
that energy is not required to create the particle) 
and photons have zero rest mass. The ratio of densi- 
ties of states favours internal conversion most 
strongly just above threshold and least strongly at 
high energies. The details of the process also favour 
conversion just above threshold. 

The internal conversion coefficient a is defined 
as the ratio of probabilities of the conversion pro- 
cess to y-emission. It can be further divided into 


components qa = Ohi + a etc. corresponding to 


L * Op 
the shell from which the electron is emitted. Its 


behaviour is shown graphically in Fig.5.2. The co- 
efficients Air Op» Ans etc. have been extensively 
tabulated for all atoms over a wide range of ener- 
gies and multipolarities. Measurements of a, dys 
Gp > %, are used to determine the type of transition 
and even the amount of mixing in a mixed transi- 
tion. Since the process is one of two-body break- 
up, the energy spectrum consists of discrete 

lines separated by the binding-energy differences 
of the atomic shells. Such lines were first dis- 
covered when superimposed on continuous £-spec- 

tra and initially led to confusion in the inter- 


pretation of §8-spectra. 
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hy energy / mar” 

Fig.5.2. Internal conversion coefficients (a) as functions of 
energy and multipolarity for a heavy element. (Based 
on SEGRE (1959). Experimental nuclear phystes, VOL.III, 
Wiley.) 

J = 0 + 0 TRANSITIONS 

It is not strictly correct to say that the con- 
version process has exactly the same nuclear depen- 
dence as that of y-emission. The statement arose 
because of the similarity of two expansions and the 
expansion of 1/|r,-r,| is of the stated form only if 

z lies outside the nucleus. Whilst an atomic elec- 

tron is usually well outside the nucleus, it does 

overlap it for a small fraction of the time; when 
eel eis the expansion becomes quite different, and 
in particular the monopole component, which for an 
external point gives a time-independent field because 
centre of charge and centre of mass coincide, can now 
have a time variation. Looked at simply, a monopole 


oscillation of the nucleus is a pulsating mode, having 


no effect on an external point, but for an internal 
point the potential depends upon the instantaneous 
nuclear radius and this is oscillating in time. 
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Although the finite size effect is small for 
transitions capable of radiating, for 0 > 0 transi- 
tions it is crucial. There are examples of such 
transitions of less than 1 MeV which can de-excite 
by no other process - the simultaneous emission of 
two y-rays is possible but has never been detected: 


An example is the first excited state of 1266 


(07, 690 keV) which decays in this way with a half- 
life of 0°3 us. Above 1 MeV the process in the next 


section becomes possible. 


INTERNAL PAIR CREATION 

In this process a state of excitation greater 
than 2 m C4 (1-02) MeV) can de-excite by the crea- 
tion and emission of ane -e pair. Since the 
emission can take place from inside the nucleus, 0 + 0 
transitions can occur this way. Also the density of 
final states will increase more rapidly with energy 
for this three-body break-up than for internal con- 
version, so at high enough energies this process will 
dominate in 0 > 0 transitions. An example is the 
first excited state of 16) (0°, 6 MeV) which decays 


-l1 . : era 
S. Pair-emission in com- 


by pair-emission in 5 x 10 
petition with y-emission has also been observed. 

Again in principle it should dominate at high energies, 
from phase-space arguments, but long before it does 

so the levels have become virtual with many channels 


open for the emission of nucleons. 


PROBLEMS 


5.1. Estimate the lifetime of the first excited state 
of a (spinless) proton in a cubic box of side 
S fm. (Consult Pauling and Wilson (1935) for 
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the quantum-mechanical problem.) 


5.2*.What is meant by multipolarity of electromagnetic 
radiation? Explain why multipolarities other 
than El are commonly observed in nuclei but only 
rarely in atoms. 


The decay scheme of low-lying levels in 14, is 


+ + 


(A; 0; 1,0; - ), (B; 2°33; 0 , 1; A, 100 per 
cent), (E; 5°1; 2 , 0; A, 67 per cent; B, 33 per 
cent), (C; 3°9; 1%, 0; B, 96 per cent; A, 4 per 
cent) (D; 4°9; 0 , 0; A 100 per cent), (E; 5:1); 
2, 0; A, 67 per cent; B, 33 per cent), where the 
capital letter designates a level, of energy in 
MeV given next, followed by 3g", 7, and finally the 
y-decay of the level in percentages to the levels 
below. Sketch the level scheme and the 
transitions. Explain the significance of J", 7 
and give the possible and expected multipolarities 
of transitions. Comment on the absence of the 
transitions DC, DB, ED, and EC. What other decay 
processes can occur? 

160 160 


65 66>Y° 
Two of the y-rays resulting are associated with 


Tb decays by f-emission to states of 


an internal-conversion line spectrum as follows: 
32°9, 78-0, 84:7, 86°1, 143-0, 188-6, 195-1, 
196-5 keV. The K,L.M,N, edges in Dy are 53°8, 
8°6, 1°9, and 0°4 keV respectively. Find the 
yY-energies, identifying each conversion line with 
the nuclear transition and the atomic shell from 
which it arises. (Ignore the splitting of the 
L,M,N levels.) 


5.4. Determine the effective charges due to centre-of- 


mass effects for neutron and proton El and E2 
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transitions in W = Z nuclei, using the argument 
outlined on p.127. 


From the result of Problem 5.4, expressing the 
dipole operator in W = Z nuclei as -e y t,n,, 
show that no El transition is possible between 

T = 0 states. (Consider specifically the isospin 
component of the wavefunction of a proton and a 
neutron outside closed shells.) Return to 
Problem 5.2 to determine whether this new selec- 


tion rule has relevance. 


6 . Nuclear reactions 


INTRODUCTION 

In the earlier chapters the decay of unstable ground 
states was considered, and in the previous chapter the 
decay of low-lying excited states. The scope is ex- 
tended here to states of high enough excitation to 
decay by particle-emission and also to consideration 
in more detail of the setting up of the low-lying 
states which decay by y-emission. The experimental 
technique is to bombard a suitable material (the 
target) with a beam of nuclear projectiles and to ob- 
serve the reaction products, or to select them for 
observation using coincidence techniques. If the 
projectiles are charged particles (p, d, a, ne. OO: 
etc.) they can be produced in quantity in a well- 
focused beam of well-defined energy from an accelerator, 
e.g. a cyclotron or an electrostatic generator. If 
the projectiles are neutrons, then they must first be 
produced by a nuclear reaction, selected in energy as 
well as possible, and collimated (rather than focused) 
onto the target; the products from the second nuclear 
reaction are then observed - against the background 
produced by neutrons from the first reaction emitted 
in directions other than that of the second target. 

A copious source of neutrons is a reactor (see Chapter 
7), especially for neutrons of low energy (the keV and 
eV region); but for neutrons of energy a few hundred 
keV and greater it is preferable to produce them 

using a primary reaction induced by charged particles 


from an accelerator. 
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The study of nuclear reactions has three basic 
aims: 


(1) to locate levels in a nucleus; 

(2) to determine nuclear properties of levels; 

(3) to study reaction mechanisms, i.e. nuclear dy- 
namics as distinct from the (quasi-) static 


behaviour in (2). 


Experimentally there is much overlap of the three 
aims, especially (1) and (2), since the location of 
levels and decay mechanisms can often be determined 
in the same experiment. To locate levels we need to 
be able to handle the kinematics of nuclear reactions; 
to determine their properties we need to look at the 
decay modes and measure their lifetimes. For (3), 

it is intended to limit the discussion to a brief 
description of the direct and compound nucleus re- 
action mechanisms. 


THE KINEMATICS OF NUCLEAR REACTIONS 
A nuclear reaction can be represented as 


A+a->B+b + Q, 


where a is the incident projectile and A the target 
nucleus, at rest in the laboratory frame; b is a light 
particle emitted in the reaction and B the final 
nucleus (the distinction between light and heavy final 
product is one of convenience; it is not always pos- 
sible to make it - an example being fission (see 
Chapter 7)). The ‘heat of reaction’ @, as defined 
above, is positive if the kinetic energies of the 
final products exceed that of the initiating project- 


ile. The extra energy comes from the internal poten- 
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tial energies of the particles, and these are reflec- 
ted in their masses. Elastic scattering is a par- 
ticularly simple example of a nuclear reaction, having 
@ = 0. In general, the final nucleus need not be in 
its ground state; an increased state of excitation 
results in a reduced Q-value. From the kinematics 

of the reaction the @-value is determined and so the 
state of excitation of B, provided the ground-state 


transition can be recognized. 


MpEb @ 

N 

\ 
\ 
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Fig.6.1. Schematic representation of the kinetics of the 


reaction A+a-7.>B+tb. 


If the reaction is represented in the laboratory 
frame as in Fig.6.1, then using non-relativistic 
mechanics (accurate enough for most reactions using 


electrostatic generators and small cyclotrons): 


conservation of momentum: 
¥(2m ,£,) = ¥(2m. Ey )cos 8 + ¥(2mpE_) cos $ 


0 = ¥ (2m, F,) sin 8 - ¥(2mpE,)sin 6, 
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which on re-ordering, squaring, and adding give 
ml 


Bp = ™,F, +t m Ey, - 27 (mm EE) cos Qe, 


conservation of (kinetic) energy: 


giving 


(my +m, Ey - (mp-m JE, - 2v(m mE JE.) cos @ - my @ = 0, 


(6.1) 
conservation of (mass) energy: 
my, +m. = my + m rar (6.2) 
A a B b : ° 


Note that in the non-relativistic limit it is con- 
venient to express the relativistic conservation of 
mass-energy as two separate equations connected by 
the heat of reaction, for a reason which becomes ap- 
parent in the next paragraph. 

In the general experiment, we measure the energy 
of emission of b at a given angle of emission @ rela- 
tive to the beam. In equations (6.1) and (6.2) it 
is assumed that mis My, ™m, EL, and 8 are known (given 
data or by measurement) and mp and @ are unknown. 

Thus (6.1) and (6.2) can be used to determine Mp and 
@. In point of fact it is rarely necessary to solve 
these cumbersome simultaneous equations in the normal 
way. From examination of the mass values of the 
stable and nearly stable nuclei it will be seen that 
the error incurred by using the mass number A for 

the mass of an unknown nucleus (and A will be known by 
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conservation of nucleons) in (6.1) is probably of 

the order of 0°1 per cent, and certainly not greater 
than 0°3 per cent for all but the very lightest nuclei. 
It is therefore possible to solve (6.1) directly for 

@ to an accuracy ~ 0°2 per cent; since in (6.2) the 
term in @ is very small, the error incurred in deducing 
an accurate mass value Mp from (6.2) using this in- 
accurate @ is probably less than the errars arising 
from the quoted mass values of the other particles. 

In any case, if it is justified to work to greater 
accuracy relativistic effects should be considered. 

If we have measured a number of different ener- 
gies of particle b at a given 8, corresponding to 
different energy levels in the final nucleus B (as- 
suming that b itself has no excited state - true for 
n, p, d, «), then to each energy E, can be assigned 
a @-value @,, where Q9-@ is the energy of the level 
in nucleus B. Qo» the @-value of the transition to 
ground, is either known or recognizable from a pre- 
viously known pattern of excited states. Often the 
largest @-value is assumed to arise from the transi- 
tion to ground - occasionally an inaccurate assump- 
tion. 

Often equation (6.1) is used in reverse. Ina 
given spectrum we may suspect that there may be spur- 
ious lines due to impurities in the target. We there- 
fore wish to determine energies corresponding to 
levels in the final nucleus following reaction with 
the target impurity, for which Qo is known and the 
energy levels of the final nucleus. Thus in (6.1) 
all is known except Ey and this must be solved either 
as a quadratic in (Ey) 2 or by iteration. 


Finally, in measuring differential cross-sections, 
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it must be appreciated that the correct coordinate 
System to use is the centre-of-mass frame, not the 
laboratory frame. This affects the measurements in 
that the transformation from the latter to the former 
frame changes both the angle and the effective solid 
angle of the detector. The calculation has been left 


as an exercise. 


PROPERTIES OF LEVELS 

Having located (bound) levels in a nucleus using 
the above kinematic analysis, it is desired to deter- 
mine other characteristics of these levels such as 
spin, parity, y-decay spectrum to lower-lying levels, 
half-life for y-decay, magnetic moment, electric quad- 
rupole moment, and less obvious characteristics such 
as the extent to which a particular level looks like 
a single particle (proton or neutron) outside a core 
which could be the ground state (or an excited state) 
of the appropriate neighbour. 

A great deal of information can be obtained by 
detecting the y-radiation from the state, either 
directly or in coincidence with the particle which 
resulted in the formation of the state. The relative 
intensities of y-transitions to known levels can give 
useful information on application of the selection 
rules. By the method of 'delayed coincidences' it 
is possible to plot the decay curve of the state and 
so determine its half-life. Used in conjunction with 
Weisskopf estimates (see p.125), it may then be pos- 
Sible to determine the spin and parity of the state. 
A useful technique is to measure the angular correla- 
tion of the y-radiation either relative to the bean, 


with the outgoing particle detected in a certain direc- 
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tion (a most convenient direction is along, or against, 
the beam, since the problem then has axial symmetry). 
In order to see how information can be obtained 
from the direction of emission of particles from the 
nucleus, consider an idealized nuclear state which 
looks like an a-particle in an 2? = 1 orbit outside 
a closed shell core. Then, relative to a given axis, 
the angular components of the wavefunction will be of 
the form Ya (0s). If the reaction setting up the 
state was selective in m and ensured that it was in 
the m = 0 sub-state then the angular wavefunction is 
keds cos 8; if also this state is unstable to 
the emission of a-particles, then the wavefunction of 
the free a-particle must have the same angular depen- 
dence (a necessary condition for matching wavefunctions 
at a spherical boundary). So the probability of ob- 
serving an a-particle at an angle ® to the given axis 
will in this simple case be proportional to cos? 0. 
It is unfortunate that the example chosen is that of 
an unbound state when discussing the y-decay of par- 
ticle-stable energy levels, but the object was to 
avoid the additional complexities of spin when con- 
sidering effects of orbital angular momentum and this 
is not possible for y-transitions (nor nucleon emis- 
sion). Angular dependence is, however, implicit .in 
the discussion of y-emission in the previous chapter, 
where it was indicated that y-absorption by a 
(J=0,M=0) orbital state of a proton to form a (1,0) 
orbital state arose from the zg-component of the radia- 
tion field. Since the field is perpendicular to the 
propagation vector, the z-component will contain the 
factor sin®, where ® is the angle between the pro- 
pagation vector and the z-axis. Thus the probability 
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of absorption « sin’ 6 and, by time-reversal, the pro- 
bability of decay of a (1,0) orbital state to a (0,0) 
orbital state by emission of a photon in a direction 
8 to the axis of quantization « sin’ ®& Another way 
of presenting this is merely to state the classical 
result that a dipole along the z-axis emits a radiation 
field whose amplitude is proportional to siné. 

Thus, in general, if a level can be excited in 
a single sub-state relative to a given axis then the 
decay products (whether particle or y-ray) will be 
emitted with an angular variation relative to that 
axis. A corollary to this statement is that if the level 
has equal population of all sub-states then the radia- 
tion will be isotropic, because lack of knowledge con- 
cerning the formation of the state can only be ex- 
pressed in this way. It is rarely possible to set 
up a decaying level as a single sub-state, but neither 
is it necessary to do so to observe an angular depen- 
dence of the decay; a non-random population of the sub- 
states is sufficient, and this occurs whenever the 
particle initiating the reaction is absorbed out of 
a directional beam. Classically, the absorption of 
a particle moving in a certain direction cannot change 
the component of orbital angular momentum of the 
system along that direction and this holds good in 
making the transition to quantum mechanics. It is 
therefore convenient to make the beam direction the 
axis of quantization; the setting up of the state by 
absorption of a particle out of the beam is then 
characterized by Am = 0. Thus if an a-particle were 
to be absorbed by a 0” target to form al state then 
this state could only be formed as the (1,0) sub-state, 


and its subsequent decay by y-emission to a 0” ground 
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state will give an angular distribution «= sin’ 0. If 
the target ground state has angular momentum J and the 
a-absorption produces a state J’ = J + 1 then the sub- 
state (J,M) will result in the formation of (J + 1,M). 
Since nothing is known about the target condition the 
angular distribution arising from each initial (J,M) 
must be summed with equal weighting. But in the inter- 
mediate state, the sub-states (J + 1, J + 1) and 

(J + 1, -J - 1) cannot be formed at all, so the avera- 
ging in the initial state is not equivalent to random 
averaging in the intermediate state and is not likely 


to produce isotropy. To proceed further requires a 
detailed study of the mathematics of coupling of angu- 


lar momentum in quantized systems; here it suffices 
to state that an angular distribution can be cal- 
culated in terms of a few nuclear parameters, the 
J" of the states involved, the mixing of emitted 
radiations (e.g. Ml and E2), etc., and checked a- 
gainst measurement to give information concerning 


these parameters. 
Two important qualitative features can be appre- 


ciated. If the intermediate state has J = 0 then any 
radiation subsequently emitted must be isotropic - 
Since randomness is represented by equal populations 
of all sub-states, J = 0 must always be in a random 
condition. Less obviously, if the particle absorbed 
has 2 = 0, then all subsequent radiations will be 
isotropic - this follows from the realization that 
adding nothing (angular-momentum-wise) to randomness 
can hardly produce other than randomness. 

The above outlines the theme, but the variations 
are numerous and can be complicated. As an example of 
complication, if a reaction is observed in the presence 
of a magnetic field then, from Larmor's theorem, the 
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whole system can be assumed to rotate with the Larmor 
precession frequency; measuring the angular distri- 
bution as a function of time after formation of the 
intermediate state should reveal a rotation of the 
angular distribution by an amount which depends upon 
the magnetic moment of the intermediate state and the 
time delay and so gives the magnetic moment of this 
excited state. 

The methods for determining to what extent a 
given state looks like a shell-model state are best 
left until the next section. 


INTERACTION MECHANISMS 

Broadly speaking interactions can be divided into 
two categories, but the division is not sharp. The 
category of ‘direct reactions' refers to reactions 
which take place whilst the bombarding particle is 
traversing the target nucleus, i.e. in a time %v 2Ry/v, 
where v is the velocity of the particle (typically 
> e/10). This characteristic time is ioe - 10°74 
An example is the reaction A(d,p)B, in which the 
deuteron, whilst passing through the periphery of 
nucleus A loses a neutron to form B. In the second 
category, 'compound-nucleus reactions', the time of 
interaction is very long compared with the characteris- 
tic time above. Thus there is a long time interval 
between the absorption of the bombarding particle and 
the emission of the product particle. It is therefore 
reasonable to talk of the compound nucleus as an inter- 
mediate product in the reaction; if the reaction is 
A(X,Y)B then it effectively proceeds in two steps 
AtX > C* + B + Y. Since C* exists for a time long on 
the nuclear time scale it must correspond to a (quasi-) 


stationary state of the nucleus C and is in fact just 
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one of the many virtual states of excitation of C. 
As is usual, we define the processes in terms 
of extremes. The real situation is not always so 
clear. If an incoming projectile interacts with a 
nucleon near the surface of the nucleus, it or the 
nucleon could be scattered clear of the nucleus giving. 
a fast direct reaction; if it gives up some energy 
to the nucleon but insufficient to raise the latter 
out of the nucleus and thereby moves further into 
the nucleus, it can once again be scattered out by 
another nucleon or share energy with it and still 
remain inside the nucleus. At each collision it be- 
comes harder for the particles concerned to escape, 
and finally the energy is well shared among all the 
nucleons. At this stage the compound-nuclear state 
has indeed been formed. It can only decay when, by 
chance, one nucleon can acquire a sufficiently large 
Share of the excess energy to become unbound, and 
this can take a long time when the energy has been 
shared over a large number of nucleons. Thus virtual 
nuclear states tend to be longer lived (for a given 
excitation) in heavier nuclei. At low energies, and 
especially for charged particles, the compound-nucleus 
process will dominate in the reaction cross-section 
(the elastic Rutherford scattering is large and 
direct) since, having penetrated the barrier (an im- 
probable process), the particle even if it fails to 
collide has little chance of penetrating once again 
at the first pass. It is therefore reflected back 
and forth at the surface until finally it shares 
energy and forms the compound nucleus. At high ener- 
gies direct reactions become of increasing importance 
- at very high energies it is virtually certain that 
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both incoming and struck particle will escape from 
the nucleus. 


The compound-nucleus process 

The reaction can be described mathematically as 
arising from a first-order time-dependent perturbation 
setting up the compound nucleus C* which subsequently 
decays, or alternatively as arising from a second-order 
time-dependent perturbation between the initial and 
final states which is dominated by one channel pro- 
ceeding through C*. The two processes are in reality 
the same thing, but it is instructive to look at both. 

As always in time-dependent perturbation, we 
consider the amplitudes a,b,e of normalized wave- 
functions of states A,B, and C* to be functions of 
time. The wavefunction representing C* is taken to 
be of the form egexp(-t/2t)b.(r,8,6) in order to re- 
present the fact that, left alone, C* will decay in 
intensity with mean life t. First-order perturbation 
theory gives the rate of formation of C* as 


-(i/h A jexpti(E,-E,)t/hia, 


where Hon will be assumed to be independent of time, 
except for having been switched on at t = 0. The anm- 
plitude ec therefore satisfies the differential equa- 
tion 


SS = -(1/20)e - (i/R)EQgexpli(Ec-Fa)t/hta . 


If the probability of transition is small, then there 
is little absorption of the wave representing the 


initial state, so a = constant = 1 is a reasonable 
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approximation, and the equation integrates to 


7 - (i/n) 4, [exp{i(Z,-F,)t/h}-exp(-¢/2t)] Poe 
(1/8) (E-E 4) + 1/2T 


satisfying the condition ec = 0 at t = 0. Whilst T 
is large on a nuclear time scale it is small on a 


laboratory scale ( % 19°24 


s or less), and so the term 
exp(-t/2t) can be neglected under experimental con- 
ditions giving 
2 lHoal” JBcal” 
Ce re 
{(Eq-Eg)S#h’/4t°} — {(EQ-Ey) 40/4} 

where [t = h . Now the probability of decay of a 
state is the sum of its decay probabilities to the 


various decay channels: 
1/t = »(1/t ,) or [Tf = , r 
2 

2 lz.,|°T 
Rate of formation of state B = Hake = os 6 —»—. 
B h{(E,-E,)°+I"/4} 

CA 

Considered as a second-order perturbation the 


rate of production of B is 


| 2 : : 
(2m/B) | LH pApc/ Eg" Fo) | x (energy density of final 
states). 


The time-dependence of the (decaying) state C* is 
exp (-i£ .t/h)exp(-t/2t) = exp{-i(Z,-il/2)¢/n}. 


Thus the decaying state C* is effectively at an energy 
(Eo - if/2). If in the above summation one term 
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dominates and that one is C*, then the expression in- 
side the moduli becomes 
2 
Pa ary ee ree 
{(Ey-Ep) +0" /4} 
Comparison of the two expressions for rate of forma- 
tion then gives 


T = 2n|H 2 x (energy density of final states), 


B Bc! 


and, by applying a time-reversal argument, 


ry = 2d Ce x (energy density of final states of 


time-reversed system). 


From Chapter 4 (p.101) this is 


where the surrounding box has been taken to have unit 
volume (in any case its size drops out), and the 
Statistical weights of the states have been included. 
Thus 


(1/mh?) |x Pia m *y (27441) (27,41) , (6.5) 
and the rate of formation of B is then 


Z 
mh (2J ,+1) Pale 
ii Se Oe OC aes ° see 3 ie 
m0, (2d +1) (27,41) (E,-EG) +T2/4 
where the statistical weight of C* has been inserted 
since the calculation has previously referred to a 
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Single intermediate state. The statistical weights 
have been introduced rather too glibly. The factor 
1/(27,+1) (27,41) is perhaps obvious since the Is have 
been defined in terms of a summation over final states, 
whereas the initial state has unit normalization 
spread over the sub-states, i.e. an average over 
initial states is required. In any case the argument 
presented is sound; the density of final states in ~ 
the time-reversed system certainly contains the fac- 
tor, but the yield does not and therefore must con- 
tain PT y/ (2d ,t+1) (27,41). At first sight it may appear 
that, at least in the expression given for the second- 
order transition, a term like (27, + 1)? arises, 
since the summation over sub-states is done before 
taking the square of the modulus. However, the ex- 
pressions derived are for rates integrated over all 
directions of emission of the final particle and all 
cross-products integrate to zero, leaving only the 
terms which arise if the summation over sub-states 
were taken after forming the square of the modulus. 
The term (27,41) is the number of 'final' states pro- 
duced by the first-order process; that these states 
decay is not relevant to the first-order treatment. 
Hence the right side of equantion (6.4) should be 
increased by this factor. 


From Appendix A, the cross-section for a process 
is defined such that Trp Ya gives the rate of forma- 
tion of state B. Therefore 


: _ Th (27,.+1) rT, 
AB 2.2 ° (27 +1)(27,+1) ° Zia dgae 
mova a A (E,-EO) +7 "/4 


(6.6) 


NUCLEAR REACTIONS 149 


and the factor in the: first brackets can be written 
more familiarly as nX?, In this expression Ey is the 
energy of the incoming beam of particles (in the 
centre-of-mass frame) and is therefore the energy 
variable. The shape of the resonance curve is shown 
in Fig.6.2, where the variation of mK? over the re- 
sonance has been neglected. 


-—2T —-T 0 r 2r 


Fig.6.2. Shape of Breit-Wigner resonance denominator plotted 
against (E-Z,. ), with energy units of I. If the re- 
sonance is sharp so that neither x? nor the widths 
vary appreciably over an energy range vu 4f, then this 
curve will represent the hehaviour of the reaction 


cross- -section. 


In the derivation of the Breit- Wigner formula 
above it has been assumed that the decay of C* is 
characteristic of Gs alone and does not depend upon 
the way C* was produced. This independence of forma- 
tion and decay of the compound nucleus was first 
stressed by Bohr. It is an important point, but has 
been somewhat overworked in the textbooks. It holds 
only for an isolated resonance, but in tests quoted 
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different formation processes can populate the over- 
lapping resonances to different degrees, and the 
decay processes will then differ. If the number of 
overlapping resonances is very large then a statis- 
tical averaging process will reduce the difference. 
Thus independence of formation. and decay is good for 
a single resonance or a large number of overlapping 
resonance, but probably does not hold for a few 


overlapping resonances. 
Two outstanding examples of resonant reactions 


are shown in Figs. 6.3 and 6.4. In the (p,y) reaction 
the resonances are sharp because the proton width is 
small due to the effect of the potential barrier. 

From the above, IT is proportional to the square of 


ii t 
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Fig.6.3. Resonances in 2781 (p.) EP es The resonance widths 
indicated represent experimental limitations; the true 
widths are considerably less. (Based on BROSTROM 
et al. (1947). Phys.Rev. 71, 661.) 
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i020 30 40 50 
En|eV 

Fig.6.4. Resonances in 238 yn, y) OU for neutrons up to 50 eV. 
Note the large cross-sections at resonance, the high 
density of levels in the compound nucleus (% 1 every 
15 eV) and the narrow widths. Note also that the 
measurements were made in absorption. Off-peak the 
absorption is due to 'hard-sphere' scattering inter- 
fering with the scattering from tails of the resonances; 
on-peak the absorption is dominated by the capture pro- 


cess. 


the modulus of a matrix element which connects a 
state representing a particle within the nucleus with 
a state representing the same particle outside the 
nucleus (a plane wave at infinity) and the overlap of 
this latter wavefunction with the nucleus introduces 
the barrier-penetrability term. It is then surprising 
that the neutron resonances shown have even smaller 
neutron widths. This is partly accounted for by the 
fact that the separation between resonances has become 
small because in the example shown the target nucleus 


in (n,y) is a heavy one and in (p,y) a light one. But 
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why should crowding of levels result in smaller partial 
widths? This question can be conveniently answered 

by considering what would happen to cross-sections 

were it not the case. Assume that the cross-section 

is being averaged over an energy range AE large com- 
pared with resonance widths, then the contribution 

to the total reaction cross-section of a single re- 
sonance is obtained by integrating the Breit-Wigner 
expression (with [ for Tp) giving a result propor- 
tional to ry only. The number of levels in AF is 

wu AE/D, where Dis the average spacing between levels. 
The averaged cross-section is therefore « T,/D. As 

the energy is increased D will get smaller corres- 
ponding to an increase in level density with ex- 
citation of nucleus C. If Ty remains approximately 
constant with energy then the cross-section increases 
without limit. But cross-sections at high energies 
should approach something like the classical TR» 
and we are forced to conclude that (T ,/D) is roughly 
independent of energy. Thus the partial width into 

a given channel depends upon the local spacing of 
levels (see also Chapter 3). This does not mean that 
the total width of the resonance is similarly limited; 
it is a sum of all partial widths and each such par- 
tial width is proportional to the local spacing, but, 
as the energy increases, the number of open channels 
contributing also increases. Thus the total width of 
levels tends to increase with energy, but only slowly 
Since the two effects almost compensate. 

To return to (n,y), it was stated that the above 
considerations only partly account for the small 
neutron widths. A closer look at the expression for 
the partial width in terms of the matrix element 
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(p.147) reveals that the width is proportional to the 
velocity in the channel and this is very small for 
neutrons of a few eV - so small that in this region 
Te < rv (This surprising circumstance is also a con- 
sequence of the fact that y-radiation is taking place 
from a level at an excitation of ~ 6 MeV or more in 
a large nucleus, so there is an extremely large num- 
ber of available states of lower energy to which y- 
transitions can take place.) In addition, the width 
includes the penetrability through the angular-momentum 
barrier (and for charged particles it is augmented 
by the Coulomb barrier) with the result that Ee in 
the eV region is far too small, when Z # 0, to give 
observable resonances. 

In the case of the (p,y) resonances, their sharp- 
ness is determined largely by the barrier factor. 
The ‘widths without barrier' in light nuclei are 
quite large fractions (~% a few per cent) of the single- 
particle proton width, but the barrier penetrability 
reduces them to a few eV giving a total width for 
the resonance less than the experimental resolution 
cf the beam. The peak yield is not therefore the 
yield at resonance, but a measure of the integrated 
yield over the resonance. 

Finally, the Breit-Wigner formula can be used 
to show that neutron cross-sections vary as 1/v at 
very low energies. So far it has been assumed that 
the [Ts are constants in the formula but this is an 
oversimplification; the matrix elements can have an 
inbuilt energy-dependence and, as has been mentioned, 
there is a factor Vo in the expression for tos For 
a sharp resonance at high energy in the incident 
channel the variation is not expected to be large over 
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the region of large cross-section, but at very low 
energies the nk? term becomes large and enhances the 
low energy tail of the resonance. It is therefore 
necessary to consider the contribution made by a re- 
sonance far from the peak of the resonance denominator. 
The cross-section is then of the form o « veT P/Q)", 
assuming that £ < E a and also lr < Eas Now T_ is un- 
likely to vary whilst the energy changes by a few eV 
since the level in C is 6 MeV or higher in excitation, 
so the energy variation is contained in ae and 


since X « 1/v, and ea 20) then o « 1/v,- 


Direct reacttons 

The characteristic of a compound-nucleus reson- 
ance is a cross-section which changes rapidly with 
energy over a small range of energy. The range over 
which the cross-section changes rapidly is of order Tr 
and is correlated with the time scale of the process 
by Tt ~h. A direct reaction has been defined as one 
taking place whilst the bombarding particle passes 
through the nucleus and therefore in a characteristic 
time » 1072s, Using the above relationship, it is 
associated with a width ~ 10 MeV. Thus a character- 
istic of direct reactions is that the cross-section 
should not change rapidly with energy, or rather should 
not exhibit sharp maxima, since it is possible for a 
direct reaction involving charged particles below the 
barrier to increase exponentially in cross-section. 

Another characteristic concerns the angular dis- 


tribution of the process. If a reaction proceeds via 
a single (and this must be emphasized), well-defined 


resonance then, at the intermediate stage of the process, 


it has a well-defined parity (the parity of the state 
C*). In setting up the state a definite direction has 
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been singled out (the direction of the beam), but 

this does not correspond to a well-defined parity, so 
in forming C* the information concerning direction 

of motion is lost but not the axis of the beam. (To 
form a state of well-defined parity it is necessary 

to superimpose upon the state defined by the bean, .a 
corresponding state with the beam moving in the op- 
posite direction.) Thus the subsequent decay of C* will 
give an angular distribution having symmetry about 90° 
(fore-aft symmetry). The direct reaction has no such re- 
striction, the symmetry axis contains the information con- 
cerning direction and angular distributions do not 

in general have fore-aft symmetry. In fact there is 

a strong tendency to forward-peaking, an-extreme ex- 
ample being the case of absorption of a fairly high- 
energy neutron beam which results in the diffractive 
scattering of the beam. As in the classical case of 
diffraction of light, the forward-directed scattering 
peak has angular width h/Ryrad. Many other types 

of direct reaction conform to this simple diffractive 
picture. 

The stripping reaction is an interesting example 
of a direct reaction. It is the reaction A(d,p)B (or 
equivalently (d,n)), in which one nucleon of the in- 
coming deuteron interacts with the nucleus whilst the 
other passes it by. The probability of the process 
occurring depends upon the square of the modulus of 
spevy,dt, where the initial wavefunction p, is that 
of the deuteron in the form of a plane wave combined 
with its internal function, multiplied by the target 
wavefunction Vane whilst the final wavefunction Ve 
will be that of the proton in the form of a plane wave 
multiplied by the wavefunction of the final nuclear 


156 NUCLEAR REACTIONS 


state bp» which for simplicity will be taken as Pavan? 
where ve is a single-particle wavefunction of the 
neutron about the core A. Notice that the Coulomb 
distortion of the plane waves has been neglected here, 
though it is included in the analysis of research 
data. For the perturbation VY we can put in the inter- 
action between the neutron and the nucleus A, which is 
taken to be of the form V(r) with the centre of A 
defined as the centre of the coordinate system. It 
is also assumed that any change of shape of the deut- 
eron during its journey past the nucleus can be ne- 
glected and that reaction only occurs when the proton 
and neutron in the deuteron are at their common centre. 
Another simplifying assumption made here is that the 
reaction takes place only at the surface of nucleus 
A - this is not assumed in a more serious calculation. 
With these approximations, in the integral, the 
internal wave-function of the deuteron becomes merely 
a constant (the probability that the proton and neutron 
should be found at its centre), and its plane wave be- 
comes exp(ik,.r), whilst that of the proton becomes 
exp(ik, .r,)> since r, = ty 
the same token V = V(Ry) = constant, as also is 


viv,- In the integral r, = (Ry,0',6') and the integral 


= Eas by assumption. By 


becomes 
constant x Sun (Cr, Jexp(iky.r,)Jexp(-ik, -1,) da". 


If k, - k = q, and the direction of q is taken to 
define the zg'-axis the integral 

« FYF (8's o' Jexp (igkcosé') da’, where the single- 
particle angular wavefunction has been inserted for 
the neutron, but its spin function has been neglected. 


NUCLEAR REACTIONS 157 


The integral over $' requires that m = 0, giving 


cos = oe 


1 7 f U 
a Y 9 68 Jexp(iqR,cos &)d(cos Qe’). 


cos 0 
To see how this behaves take the cases 2 = 0 and Zz = 1. 
The first integral, for Z = 0, can be performed at 
sight, giving constant x (sin Tey! ay ), which, when 
the substitution g? = (kg-k ~ “44k gk sin(0/2) is made, 
gives a maximum at 6 = 0° followed by subsidiary maxima 
of reduced intensity. (The function (sin qPy/qRy)” 
is recognizable as the diffraction function, so, pro- 
vided Ik : kl Fy is less than 1, the statement made 
is true. This is the case in most experimental appli- 
cations.) i in oe 9 ~ cos 8, for 2 = 1, gives 


an integral f porary e )ede. pa eines by parts 


gives (sin ay - qRy cos qRy) / (qRy)” which, on squaring, 
is seen to give zero at q = 0, rising to a large first 
maximum, and then having subsequent lesser maxima. 

With the same proviso as for 2 = 0, the angular dis- 
tribution will be non-zero at 8 = 0 (unless ka = kK 
which is possible, for low Eg and a high Q-value, 


Pp 


but not likely) rising to a large first maximum (at 

a fairly small angle in practical cases) followed by 
subsidiary maxima. The important point is that the 
location of the primary maximum can be used to deter- 
mine the Z-value of the transfer, i.e. the amount of 
orbital angular momentum taken in by the stripped 
nucleon. Thus the parity of the final state is known 
immediately on analysis of the angular distribution of 
the particle not captured, and the range of J/-value 
of the final state is limited. The degree of limita- 
tion depends upon J for the initial nucleus - if this 
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were zero then the final state must have d¢ = L ts Be 


where Lo is the transferred angular momentum. But in 
the general case J, can take the range (J; + La + 1) 


down to the smaller of |75 - 2. + 4|. Examples of 


| n 
angular distributions in (d,p) are given in Fig.6.5. 


dao(0) 
dQ 


20 40 60 80 100 
Ocm/ deg 


Fig.6.5. Typical angular distributions for the (d,p) direct re- 
action. If each process represented the full single- 
particle transfer, then relative to Z = 0, the curves 
for 2 = 1,2,3 should be multiplied respectively by 
1 1 1 


4° 12 ’ 30° 


A useful feature of the stripping reaction is 
that the magnitude of the cross-section is a measure 
of the degree to which the final state looks like the 
original state plus a single particle of the type in- 
dicated by the angular distribution. Of especial im- 
portance is the fact that the final state may be 
bound to the particle transferred (indeed the calcula- 
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tion is simpler for such bound states, and the simple 
assumptions made here would not be good for an un- 
bound state). Because the transfer is virtual rather 
than real, the real positive momentum q can be associa- 
ted with a negative total energy of transfer. 

The stripping reaction is a typical direct re- 
action, of which there are also many other examples. 
The inverse (p,d) or pick-up reaction, (*He, d), (a,p), 
etc., and their inverses are further examples of 
stripping. Other direct reactions are elastic and in- 
elastic scattering of charged particles (usually a- 
particles) off distorted nuclei; Coulomb excitation - 
the electromagnetic excitation of a nucleus by the 
rapidly changing field of a passing charged particle; 
transfer reactions, e.g. (a6: PPO). in which a 
neutron is transferred from the projectile to the 
target - this is really an example of stripping -with 


heavy-ion projectiles. 


PROBLEMS 


6.1. For elastic scattering of a projectile of mass 
m by a target of mass M show that the angle of 
scattering in the laboratory system (6) and 
the centre-of-mass system (¢) are connected by 
> = 6 + sin)! {(m/M) siné@} and that the ratio 
of solid angles between the two systems is given 
by aap /an(e) = sin2¢/sin“8 cos (¢-0). Ex- 
plain the significance of these transformations 


in the measurement of differential cross-sections. 


6.2*.A foil of ila is bombarded with protons of 
5°0 MeV, and it is found that at 90° neutrons of 


2°3.MeV are emitted. If the atomic masses of p, 


160 NUCLEAR REACTIONS 


n, and 7 


Li are respectively 1:0078U, 1:0087U, 
and 7°0160U find the atomic number and mass of 


the daughter nucleus. 


6.3*.A thin target of os is bombarded with 8 MeV 


protons, and the energy spectrum of the protons 
scattered at 45° to the incident beam contains 
strong peaks at 7:84 MeV and 6'99 MeV. What can 
be deduced about the levels of this nucleus? 


6.4. By using an argument involving parity, show that 
J = } for the intermediate state in a y-y cascade 
will, like J = 0, result in an isotropic angular 
correlation. Show also that the y-decay of an 
isomeric J = 3 state will be isotropic even when 
the state has been completely polarized into one 
of its sub-states (e.g. by application of a mag- 
netic field at low temperature), but will not be 
isotropic in this latter condition if the y- 
detector is sensitive to the state of circular 
polarization. 


6.5*.In a given nucleus, a state C (0°) decays by y- 
cascade through a state B to the ground state 
A (07). If the two y-transitions are El, what 
is the spin and parity of B? Will there be a 
non-isotropic angular correlation of the two 
y-rays? Alternatively, if the y-radiations are 
again El but B is known to be 0”, what are the 
spins and parities of A and C? Will there be 
an angular correlation? 

For the more general case, discuss qualitative- 

ly how angular-correlation measurements can help 
to give information about spins and parities of 


6 


.6* 


ps 
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nuclear states (a) if the multipole order of 
each Y-ray is known (e.g. from internal conver- 


Sion measurements), (b) if it is not. 


16 
-Suppose that the reaction tee ca,y) 0 
(2 = 7-15 MeV) shows an isolated resonance 
at By = 10°10 MeV with Dee = 0°20 MeV, 


Ty, = 0:15 MeV, and ry = 10 eV. The spin of the 
compound state is J = 4. The only other decay 
mode of 1% at this excitation is by proton 
emission to the ground state (J = 3) of 15y, 
When the same level is excited by 15u(p,y) 1%o 

(Q = 12:11 MeV) what will the peak cross-section 
be? (Remember to consider centre-of-mass 


effects.) 


(a) Write an equation embodying the principle of 
detailed balance as it applies to two-body re- 
actions with unpolarized beams and targets. (This 
should be deduced from Fermi's golden rule and the 
definition of cross-section - since these have 
been used to deduce the Breit-Wigner formula it 
can also be conveniently extracted from the 
latter, though it has greater generality.) Ex- 
plain how it may be used to deduce the spin of 

the pion from the measured cross-sections for 
ptp##d + T, (Note there is in addition a 
factor arising from the production of two iden- 
tical particles when the reaction proceeds from 
right to left.) (b) Explain why the Breit-Wigner 
formula is important in the theory of slow neu- 
tron capture by nuclei. Show that, at a reson- 
ance peak (neglecting any background), 

na’ = ga0, where o and o. are the total and 
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elastic scattering cross-sections for neutrons. 


7. Fission and Fusion 


In these days of energy crises no textbook on nuclear 
physics would be complete without some reference to 
these two sources of energy. Fusion has very little 
overlap with nuclear physics, except for the fact that 
the basic reactions are between light nuclei, but the 
nuclear reactor based on fission can truly be said 

to be the brain-child of the nuclear physicist. In- 
deed virtually the whole of nuclear physics is involved 


in describing how and why a reactor works. 


THE FISSION PROCESS 

In Chapter 1, the mass formula indicated that all 
nuclei above A ~ 150 are in principle unstable to a- 
emission. At the same time, the question of stability 
against splitting up in any way whatsoever might have 
been considered, but was not. It is so now. Splitting 
into two approximately equal parts is expected to 
be possible at large enough mass values, since the 
mass formula indicates that medium-weight nuclei are 
most stable. Expressing M(Z,A) - M(Z, Aq) - M(Z4,45); 
where Z = zy + Zo, A= 4, +A), in terms of the mass 
formula (p. 18), and maximizing this function, gives 
the result Z4 = 25, A, = Ay, and 


M(Z,A) - 2m(4z,34) = -642/3 (22/34) 4 e227 1/3(¢4-27 2/5), 


neglecting the trivial term in 6. The condition for 
stability is therefore 
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2 1/3 
a < 8 (2s ~ 0-78) ~ 17 

Thus, in principle, all nuclei above A = 90 
(Z = 40) are unstable to fission into equal fragments. 
This low value is reminiscent of the low value ob- 
tained for a-decay, which in turn leads to the reason 
why such light elements do not fragment, namely, 
the Coulomb barrier. It might be more realistic to 
expect fission to occur if the mass of the nucleus 
under consideration exceeds the masses of the two 
parts by more than the Coulomb energy when the two 
parts have just separated. Fission will then occur 
over the top of its barrier. This adds a term 


(Ze/2)? 
4ney.2Ry(4/2)° 


to the energy of the two fragments. Inserting a 
suitable value for Fy (1-3 fm), or alternatively as- 
suming that e¢ derives from a uniform charge distri- 
bution and using the value of e€ to determine the 
radius, modifies the condition to 27 /A > 50 for in- 
stability, and this would correspond to an extra- 
polation to Z ~ 130 along the valley of stability. 
These estimates are very crude. It has been ~ 
assumed that, in inverse, as they move toward each 
other the two fragments remain spherical up to con- 
tact. But nuclei in isolation can be non-spherical, 
and even if spherical it is not realistic to assume 
they are rigid. In the presence of each other it 
is possible that the Coulomb potential be less than 
the monopole term predicts because of induced dis- 
tortion or because the proton distribution inside 
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each nucleus may no longer be uniform. These and 
other mechanisms can diminish the Coulomb potential 
near the point of contact whilst not greatly affecting 
the potential of the single large nucleus. It is to 
be expected that spontaneous fission will occur for 
less massive nuclei than predicted by the above for- 
mula. 

To consider the mechanism in more detail, a po- 
tential function is needed to describe the inter- 
action between the two fragments. When the two com- 
ponents are distinct, the potential can be expressed 
as a function of the distance between their centres, 
but this parameter has little meaning before separa- 
tion. In this region a convenient parameter is AR 
for an ellipsoidal nucleus as defined in Chapter l 
(p.31), but this ceases to be useful near the point 
of break-up since AR + © produces a needle-like 
nucleus rather than two globules. If we ignore this 
difficulty the potential is shown in Fig.7.1, where 
the two dotted curves correspond to the limiting cases 
of a nucleus in principle unstable to fission but 
not observed to decay so and of a nucleus which would 
decay by the fission process in a time of order of 
the nuclear characteristic time. The interesting case 
is shown by the full line which represents a typical 
distorted nucleus in the region of uranium. 

If the well représenting the large nucleus is 
extremely shallow then it is possible for states in- 
side it to decay by penetrating the fission barrier 
(since the particle penetrating has large effective 
mass, the penetrability varies extremely rapidly with 
energy below the top of the barrier so that shallow 
means of the order of a few hundred keV). This leads 
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Coulomb 
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Fig.7.1. Schematic fission potentials. Note that the abscissa 
is the distortion parameter AR up to somewhere below 
separation and the distance between fragment centres 
above separation. The dashed curves are explained 


in the text. 


to a class of radioactivity not previously discussed, 
namely, decay by spontaneous fission. If the well 

is deeper than this, then the levels near ground may 
be stable to this mode of decay, but levels of higher 
excitation will not be, so fission will take place 

in preference to Y-emission. Such levels could be 
excited by inelastic scattering, e.g. (n,n') or (p,p').- 
At a greater well-depth, the location of the top of 
the fission barrier relative to neutron instability 
is of great importance. If the fission barrier 

lies lower then the nucleus one neutron less than 

the one here considered, on absorbing a thermal neu- 
tron, will decay by fission; if the fission barrier 


FISSION AND FUSION 167 


lies say 1 MeV higher then neutrons of energy greater 
than 1 MeV will be necessary to induce fission. 
Notice that a neutron threshold is well defined, but 
the fission threshold is effective in that just below 
it the process can take place but with very much re- 
duced probability. 

As an example of the above, spontaneous fission 
occurs in the following: 


235 17 238 1 
9 2U (Tr v~ 2 x 10° ‘years), 92U (Tr ~~ 10 Oears, 
240 11 250 


| 4 
g4Pu (Tr ~ 10°“years), gg lt (Te ~ 10 ‘years), 


254 


ggtf (Tr ~ 55days). 


Nuclei which undergo fission on capture of thermal 
neutrons are es ih ed and ao Di. whilst nuclei 
which require capture of fast neutrons are, with 


effective neutron thresholds in brackets, 


tecTh (1-3 Mev), °°" ane 


238 


U (0°4 MeV), 


U (1°1 MeV), and *OaNP (0-4 MeV). The two last 


sequences reflect the importance of even or odd 


U (0:8 MeV), 


neutron numbers, emphasizing the effect of pairing 
of nucleons. 

For nuclei in this region the fission reaction 
is highly exothermic with a Q-value of order 200 MeV. 
This is highly desirable for production of energy, 
but of itself does not lead to power production. The 
reaction must be self-sustaining as in chemical com- 
bustion when part of the energy liberated in the re- 
action is used to heat the fuel and so encourage 
further reaction. From information now to be dis- 
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closed, the fission reaction can be made self-sus- 
‘taining. Fission is not merely the splitting of a 
large nucleus into two smaller ones but, as is to be 
expected, a lot of debris is produced at the same 
time. Neutrons form part of .this debris, either pro- 
duced simultaneously with the large fragments or more 
probably by secondary emission from the fragments 
which will in general be formed in states of high 
excitation. In addition, the fragments lie on the 
neutron-rich side of the valley of stability (since 
W/Z increases with Z along the valley) and therefore 
have low binding energies for neutrons. These neu- 
trons will be emitted almost instantaneously 

(tr ~ 19°16 
This latter point can be understood in terms of the 


s or less) with a spectrum peaked at ~ 1 MeV. 


discussion on widths in Chapter 6; ro from a°given 
state to any final nuclear state is proportional 
to k or EN , but the density of final nuclear states 
has very roughly an exponential dependence upon 
energy of excitation. The most probable energy of 
emission of a neutron is therefore the maximum of 
(E-E*) 2exp(aF*), where £ is the energy of the neutron 
transition to ground and £* the excitation of the 
final nucleus. At maximum E-F* = 1/2a. Peaking at 
ED = EF-E*® ~ 0-7 MeV gives a = 0:7 in units of (Mev) 1, 
indicating how the level density behaves for a heavy 
nucleus. Note that exp (a&*) is an oversimplifica- 
tion, but it does represent roughly the crowding to- 
gether of levels at higher energies; a varies slowly 
from nucleus to nucleus. 

In addition to these 'prompt' neutrons there are 
also emitted some delayed neutrons amounting to 
~v 0:5 per cent of all neutrons produced. These were 


FISSION AND FUSION 169 


initially a source of mystery, but are now known to 
occur following B-emission; an example is the decay 


of the fission fragment 87Ry which B-decays with 55s 


35 
half-life to eek, with an available energy ~~ 6 MeV. 
As well as decaying to nucleon stable levels in 87 Ky 


a fraction of the 8-transitions are to states of ex- 
citation greater than 5-4 MeV, and these can decay 
by neutron-emission. Neutron-emission will therefore 
87 

Br. 


Although. relatively few, these delayed neutrons play 


be delayed by the time taken for §-decay of 


a large part in the control of a reactor. If. each 
neutron-induced fission event produces on average 
more than one neutron, the possibility of a self-sus- 
taining chain reaction arises. Typically the number 
v, per fission (average) is something like 2-2:5. 
Before going on to discuss how to make use of 
the fission reaction, comment should be made on the 
reference to the fragment A = 87. This is a long way 
from the symmetric splitting (A ~ 115) which has been 
suggested as giving the greatest energy release. In 
fact the most probable mode of fission is not the 
symmetric mode, but corresponds roughly to an A = 90 
and 140 split. The explanation of this is still re- 
quired; magic numbers probably play a part since 
WV = 50 for A % 90, W = 82 for A ~ 140. Note that the 
Coulomb term is less for the asymmetric split rela- 
tive to the symmetric, offsetting to some extent the 
mass balance which favours the symmetric split. Pre- 
diction of the 'mass spectrum' requires detailed 
knowledge of the fission process - any mechanism 
which can change the fission barrier height, as a 
function of the fragment masses, by © 0°1 MeV in 
200 MeV could account for the spectrum. It is inter- 
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esting to note that, close to threshold, the (90, 140) 
splitting predominates and there is virtually no 
Symmetric splitting. As the energy of the neutron in- 
ducing fission rises, the splitting becomes more 
nearly symmetric in accordance with the rough guide- 


lines above. 


THE CHAIN REACTION AND REACTORS 

If a material fissile to fast neutrons produces 
v neutrons per fission (on average) and the probability 
that, on absorbing a neutron (and the definition can 
be extended to include absorption in other materials 
present), fission occurs is p (again on average) then, 
if py > 1, a chain reaction will occur in a large mass 
of this material, culminating in an explosion. The 
intensity of the explosion depends on the rate of libe- 
ration of energy compared with the rate at which ener- 
gy can be transmitted away (possibly taking the 
material with it). If the interaction length for 
fission is X then the interaction time is A/v, and 
in that (mean) time the neutron multiplication is 
(pv-1), giving an exponential growth in neutron flux 


with time constant A/v(pv-1). For aun Or vu. 1b, 


9 1 


giving an interaction length ~ 20 cm, v v 10°cm s_ 
(Z, ~ 1 MeV), and for pv-1 % 1 this gives a time con- 
stant % 108s, For such a rapid rate of neutron 

growth and therefore energy production, heat dissipa- 
tion is negligible, and because of inertia the material 
has no time to get away. The result is the so-called 
‘atom bomb'. From what has been said the critical 

size for explosion is something greater than 20 cn, 

and this can be achieved by imploding (chemically) 


some smaller blocks of 2357), 
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On the face of it, it seems unlikely that such 
a system could be controlled. The fissile material 
could be diluted and pv could be made very close to 
unity, but it would be difficult to produce just the 
correct conditions for a slow enough growth rate in 
order to be able to introduce some form of mechanical 
control, and if by error conditions are changed to 
make pv slightly greater then the system runs amok. 
At this point the fact that ~ 0°5 per cent of the 
neutrons are delayed becomes crucial; if pv is just 
less than unity for prompt neutrons and just greater 
than unity for prompt-plus-delayed neutrons. then the 
time constant of the delayed neutrons takes control 
and the growth-rate time constant will be considerably 
greater than this (since pv-1 ~ 0-005). Thus there is 
a range of variation of pv of ~ 0-5 per cent over which 
the neutron flux is increasing at a controllable rate. 
For the fast reactor, which is effectively being dis--= 
cussed here, control would be achieved by moving 
parts of the fissile material away from the rest or 
by moving some of the surrounding material relative 
to the fissile core, thereby altering the number of 
neutrons being reflected back into the core. 


Historically the thermal reactor preceded the 
fast reactor. Why was this the case? The answer 
lies in understanding what makes a suitable fuel. 
Although Zo and some other nuclei listed on p. 
167 are fissionable above a neutron threshold they 
are not suitable fuels for a fast reactor, whereas 
nuclei like 2351. which are fissionable to thermal 
neutrons, are. For these latter nuclei the fac- 
tor p depends only upon the ratio o(n,y)/o(n,f), 


since o(n,n'), which is larger than o(n,f) at MeV 
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energies, does not remove neutrons but merely changes 


238), were used as a fuel 


their energies. However, if 
then (n,n') would rapidly remove neutrons below the 
fission threshold and so be akin to absorption, with 
a reduction in p. This process is enhanced by the 
fact that the reactor has to have a structure and @ 
coolant so the neutron spectrum inside a fast reactor 
is degraded to lower energies very quickly. It is 
therefore necessary that the fuel be fissile down to 
low energies. Thus natural uranium, containing only 


235 


0-7 per cent U, cannot be made to go critical; 


the fuel must be highly enriched in sca te a costly 
undertaking. 

It would therefore seem surprising that any 
structure containing natural uranium can be made to 
go critical. If, however, we look at the various 
cross-sections at thermal energies (~ 0:025 eV) we 
find that, for 2354 o(n,y) ~ 101 b and o(n,f) » 577 b, 
eo2 o(n,y) © 2-76 and o(n,f) ~ 0. Even 
taking into account that there is 140 times more 2995 
than foi in natural uranium, a value p %v 0°54 eet 


tained for thermal neutrons. Since v v 2°44 for U, 


whereas for 


it is possible to obtain criticality provided that not 
more than a quarter of the neutrons are lost in ther- 
malizing them. Before considering the thermalization 
process it is pertinent to ask why o(n,y) for ooo 
turns out to be so much smaller at thermal energies 
than the Bo cross-sections. This is partly a for- 
tunate accident, but not completely so. 

On comparing the excitation functions of U +n 
for the different isotopes, an odd-even effect is 
immediately apparent. The even isotopes produce 


narrow, well-separated resonances having very high 
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peak cross-sections; the odd isotopes produce broader 
resonances, poorly separated (the density of reson- 
ances is %10 times greater for the odd than the 

even isotopes) and with lower peak cross-sections. 
Previously it has been stated that level density in- 
creases exponentially with excitation, and the odd- 
even effect does ensure that the excitations produced 
by absorption in an odd isotope are about 1°5 MeV 
higher than in an even isotope reflecting the fact 

that in the former case the resulting ground state 

has one more pair of neutrons than the initial ground 
state (ei + Tip gives 6°48 MeV excitation, aij + Tip 
gives 4-78 MeV). But, because pairs must be split 

in order to create excited states there is an odd- 

even effect needed in applying the exponential func- 
tion. It is probably more correct to translate 
energies such that neutron thresholds coincide in 
comparing level densities of neighbouring isotopes. 
There is, however, a (2/+1) factor appearing in the 
level-density function, and this can account for the 
observed ratio, since 2351) + n produces states of 
J" = 3° and 47 (27354 ground state is 5 and Ln = 0 
for neutrons in the eV range) whilst ej + n pro- 
duces J" = 3" only. Merely increasing the level-_ 
density ratio should not (on average) enhance the 
capture cross-section ratio at thermal energies, 
since ony o rT /E (assuming that E > T - see the 
Breit-Wigner expression, p.148), and this can be 


written 
(T,/D) (1, /D) 
(E,/D)* 


Provided the top two factors remain much the same 
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from one isotope to the other it would appear that 

the thermal cross-section depends upon how close is 
the nearest level compared to the level spacing - and 
this is how chance creeps in. However, though 

(T /2) is not expected to vary much between neigh- 
bouring nuclei since it refers to a single channel; 
the same is not true for (T./D), because it represents 
a summation over a large number of y-transitions. The 
same (2J+1) factor, which produced a greater density 


~ 235 
of resonances in 


U +n, also provides a greater 
density of states to which y-transitions can be made 
since strong y-transitions change J only by one unit 
at most. It therefore turns out that total radiation 
widths ris following neutron absorption, tend to be 
Similar for neighbouring nuclei, rather than Tain: 
There is therefore a tendency for even-A isotopes to 
have lower thermal-capture cross-sections than odd-A 
isotopes, apart from any odd-even effects which may 
also occur because of pairing. This argument is 
little affected by the presence or absence of fission 
width, though it may need slight modification if the 
increased total width no longer satisfies Bs eT 
Returning to thermalization, neutrons always lose 
energy in elastic collisions, and the lighter the 
scattering nucleus the greater the loss. Hydrogen 
would therefore seem to be the most suitable material, 
conveniently perhaps in the form of water; but, un- 
fortunately it has a relatively high thermal-neutron 
capture cross-section (v 330 mb), and so absorbs many 
of the neutrons on slowing them down. It cannot be 
used with natural uranium, but uranium enriched in 
222i can go critical in a light (ordinary) water en- 
vironment, and modern power reactors have been con- 
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structed to this basic design. The next suitable 
material is heavy water since ay is the next best 
thermalizer, and, moreover, o(n,y) is very low (0°5 mb). 
Again, modern power reactors have been built using 
natural uranium fuel and heavy water. The next most 
convenient material is carbon, which is much poorer 

at moderating but still has a small enough o(n,y) to 
be of value (3-4 mb). The early experimental reactors 
were built with carbon (graphite) as moderator and 
natural fuel, and modern power reactors have been so 
based. 

The geometrical arrangement of fuel and moderator 
is important. The obvious homogeneous mixture is un- 
suitable for natural fuel since the neutron capture 
process in £29 j is important right down to its lowest 
resonance at EF, % 6-6 eV (see Fig.6.4, p.151). Sepa- 
ration of the two processes of fission and moderation 
is needed, but the next obvious step of having all 
the fuel at the centre surrounded by moderator is also 
unsuitable. At thermal energies the effective cross- 
section for removal of neutrons is 
2°7 + (101 + 577)/140 ~ 7-5 b, weighting the cross- 
sections according to the isotopic abundance, giving 
an absorption length of order 2 cm. Thus the thermal- 
neutron flux into the fissile centre is strongly 
attenuated, and most of the fuel is inefficiently 
used. The solution is therefore a matrix of natural 
fuel rods a few millimetres in diameter embedded in 
the moderator, through which must also be circulated 
a coolant with low neutron-absorbing properties to 
transfer the heat produced to the outside world. At 
least this is the case for the graphite-moderated 


reactors, the suitable coolant being carbon dioxide 
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gas; the liquid moderators can, in addition, be used 
for heat transfer. Note that, because the fission 
fragments have such a short range, virtually all the 
heat is produced within the fuel itself, with a small 
fraction arising from slowing down the neutrons and 
absorption of y-radiation. 

The control of a thermal reactor derives once 
again from the existence of delayed neutrons, though 
the time constant when critical to prompt neutrons 
is rather slower because the neutrons, when thermal- 
ized diffuse rather slowly through the moderator 


: e so diffusion speeds 


(thermal speeds are v 10° cms 
are lower and distances of traverse to get back to 
fuel are ~ 10 cm giving a (prompt critical) rise 
time ~ ms). The mechanical control is in the form 
of cadmium rods that can be inserted into the reactor. 
Cadmium is a very strong absorber of thermal neutrons 
by virtue of the large value of o(n,y) (~ 20,000 b) 
for the isotope ree of abundance 12:3 per cent. 
This high cross-section arises from a resonance very 
close to ee ce energies, as has been discussed pre- 
“U. 


So far only the rare isotope 


viously for 
2351 has been 
238) could 


also be used resources of power are considerably in- 


utilized as a nuclear fuel. Obviously if 


creased. This can be done, making use of the very 
process that has proved troublesome in the design of 
reactors discussed, namely, the capture process to 


239. 239 


form U undergoes two quick 8 -decays to 


239Du, which is long-lived (% 24,000 years) and has 
already been listed among nuclei which are thermally 
fissile. This breeding of nuclear fuel can be best 


achieved with a fast reactor since v tends to rise 
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with neutron energy and the greater its value the 
greater the surplus of neutrons for breeding over 
and above the requirements for sustaining the reaction. 
This surplus is made use of by making the core small 
enough to allow a large leakage of neutrons into a 

surrounding uranium blanket. Since the values of pv 
are 2:45 and 2:70 for pure ey and se 


reasonable loss incurred from other materials and 


Pu, even after 


coolant, within the core, it is possible to breed 
more new fuel than is consumed. It is necessary, of 
course, to remove the uranium blanket from time to 


239 


time to separate out chemically the Pu, which is 


behaving effectively as a catalyst in the fission of 
Zoey Thus, in principle, all of natural uranium 
can be usefully consumed and. at ~ 200 MeV per nucleus 


this works out at about lg per day per MW! 


FUSION IN THE SUN 

Before discussing the nuclear physics behind 
attempts to produce fusion in the laboratory, it is 
of interest to describe briefly the nuclear reactions 
going on inside the nearest working fusion reactor, 
the Sun. Prior to the discovery of nuclear reactions 
the Sun's source of energy was a mystery. The rate 
of emission of radiant energy from the Sun is well 
known and so also is its mass; gravitational collapse 
of this mass provides a large source of energy, but 
nowhere near large enough to keep the Sun radiating 
for a few thousand million years, and there is geolo- 
gical evidence on Earth to indicate that.it has done 
so for a considerable fraction of that time. On the 
other hand, gravitational collapse to the present size 


from a tenuous gas will provide enough energy to heat 
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the Sun up to a temperature % 10” K even allowing for 
radiation losses in the slow process. At that temp- 
erature, corresponding to a mean thermal energy 

~ 1 keV, protons are beginning to collide sufficiently 
violently with each other to penetrate the Coulomb 
barrier. The stumbling block to the onset of nuclear 
reactions culminating in the synthesis of helium from 
hydrogen would appear to be the initial reaction of two 
colliding protons. If we assume that two. protons 
momentarily make a compound nucleus, the diproton, 
then this nucleus will decay back into the entrance 
channel thereby producing no reaction. However, a 
parallel though improbable decay is by g*-emission to 
the deuteron. Improbable though it may be, the Sun 
contains a lot of colliding protons which cannot dis- 
appear by any other process, so this one must occur 
(though it has never been observed in the laboratory) 


initiating the sequence: 


ptprdte #v 


3 


p + dw-w+ “He + ¥y 


SHe + sHe > 4 


He + 2p 
Taking the first two reactions twice, the nett result 
is 


‘He + 2e* + 2v + Q, (7.1) 


4p > 
and from the mass values this liberates  ~ 7 MeV per 
proton absorbed (by comparison fission gives less than 
1 MeV per nucleon involved). The choice of reactions 
after the first is dictated by the preponderance of 


| 


I 
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protons. Thus the possible third reaction 


d+ He + 4 


He + p will not occur since all deuterons 
formed will be mopped up by the second reaction before 
they can get near a *He nucleus. On the other hand, 
the SHe + p interaction can only decay back into the 
same channel, so SHe is stable in the presence of 
protons and will build up in concentration until the 
third ‘reaction takes place at the appropriate rate 
for equilibrium. 

The rate at which a reaction proceeds depends 
upon 1 (V4) M4 (%Q) 90%, 9) 405 where ny and nN, are 
velocity distributions of densities of particles 1 


and 2 per unit energy range and v is their relative 


velocity. The two Maxwellian eee can be 
combined into a single such function in terms of the 
reduced mass and relative velocity giving (dropping 
the suffixes) n1(v4)n,(%)%1> o vy? exp(-mv*/2k?). 

o(v) can be written as Oy -P(E). (ER/E), where P(E) 

is the Gamow factor, exp(-272Z,2,e4/4ne) v), in the 
penetrability for the limiting case of an energy EF 
well below the barrier height Ex, and therefore for a 
nuclear radius which is very small compared to the 
classical distance of closest approach. (E,/E) 
includes the nae term in the cross-section. Expressed 
in this way 54 1s the cross-section at —F % Ep and is 
of the order of 0-1 b for a decay by particle emission 
and perhaps ~ 10 ub for decay by y-emission. 

The reaction rate is obtained by integrating over 
all v. The two key factors will be the rapidly falling 
exp(-2/E,), where Ens the thermal energy, is ~ 1 keV 
at 107 K, and the rapidly rising exp(-8£ 7) where, for 
Z, = Z, = 1 and EF in keV, 8 ~ 30. Maximizing the 


1 Z 
product of these two functions (the rate of variation 
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of other factors can be neglected) gives pol? ~ BE./2, 
giving E ~ 6 keV. Thus in the region around 6 keV 
there will be a sharp peak in the reaction rate (see 
Fig.7.2). At this energy the penetrability factor 


is v exp(-12) or 10°> and the cross-section for 
10 


ptd- sHe + y will be v 10° b; the reaction 


4 5 
EjkT 


Fig.7.2. An indication of the peaking of the yield curve of a 


thermonuclear reaction. 


4He + *He ag tHe + 2p, having Z, = 45% 2, will peak 
“10 4.) For 


the p + p reaction the cross-section will in addition 


at ~ 15 keV and give a cross-section %v 10 


contain a factor representing the probability of ae 
decay during a nuclear characteristic time. A g*- 
lifetime ~ 10° 
the 8° -decay of the neutron), and a nuclear time 

~ 1041s gives the factor Tn gia and a cross-section 
at peak of about 10°?" b. Notice that none of these 


cross-sections is measurable in the laboratory; 


s is a rough guess (about the same as 


measurements could be pushed down to v 30 keV in 


FISSION AND FUSION 181 


the second and ~ 100 keV in the third reaction and 
the energy-dependence used to extrapolate further 
down. In the first reaction the factor 10 24 pre- 
cludes measurement at any energy so the cross-section 
must be derived theoretically. 

An alternative set of reactions to the above 


is the carbon-nitrogen cycle: 


I ne 9 os PON , Diva oe cay 
ee a ae pence Mee 
eo ee eee Uy 4 ps og He 


#26 is regenerated and acts as a catalyst. 


Thus the 
The importance of the two cycles depends upon the 
temperature and upon the amount of reG present. 
(Presumably the amount present in the primordial gas 
from which the Sun condensed, since formation of 
heavier elements requires a considerable concentration 
of tHe in order to start the process with 

3"He > 120 + y. This reaction is necessary since 

p + que merely returns to the same channel, as does 
tHe + He?) At high enough temperatures (more massive 
stars) this latter cycle will take over, since none 

of the processes is particularly inhibited except 

for the Gamow factor, which will be larger because 


of the increased charges on C and N. 


FUSION IN THE LABORATORY 

Neither of the cycles described are suitable 
bases for terrestrially controlled fusion. The hy- 
drogen cycle begins with the improbable (weak) re- 
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action of two protons whilst the carbon-nitrogen cycle 
requires higher temperatures which must be maintained 
over periods long compared with the 10 min life-time 
for B* -decay of T5y, The chief problems are the 
heating up of the reacting materials, and, more 
difficult, the containment of the hot material whilst 
sufficient reaction occurs to make a nett liberation 
of energy. A gas at 107 K (or strictly a plasma, 
since it will be highly ionized) cannot simply be 
retained by walls since, no matter what the material 
of the container, it will not be able to withstand 
bombardment of particles at such energies either 
physically or chemically. Much work has been carried 
out using ‘magnetic bottles'; since the particles in 
the plasma are charged they are deflected by magnetic 
fields, and a great deal of cunning has gone into 
designing the field so that the particles are kept 
away from the walls of the vessel. In fact the mag- 
netic fields can be used to feed energy into the 
plasma and to compress it into a smaller and smaller 
volume thereby increasing the reaction rate both by 
temperature and density increase. Unfortunately, the 
system becomes unstable, and so containment cannot 
be maintained indefinitely. Containment times are 
at present measured in fractions of:.a second, which 
of course can be large enough to liberate enormous 
amounts of energy as in the atom bomb or its later 
development the hydrogen bomb, which is merely an 
atom bomb surrounded by fusible material; the atom 
bomb raises the temperature high enough (% 10° K) 
and fusion then occurs more rapidly than the energy 
and the material can be dissipated. 

Obviously the reaction, or cycle of reactions, 
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must be intrinsically fast. A suitable one for the 
laboratory would be 


d+de-> He +n 
(7.3) 
> T Dp, 


Since it uses deuterium only which occurs as 1/7000 of 
natural hydrogen and so is available in unlimited 
supply. The reaction involves unit charges and will 
therefore go at the lowest necessary temperatures 
(10/7 K in the Sun, but probably 108 K for short con- 
tainment times). These two reactions have rather 
low @-values (+ 3:27 MeV and + 4-03 MeV), releasing 
only ~ 1 MeV per nucleon. In addition, they have 
relatively low peak cross-sections, ~ 100 mb peaking 
at several hundred keV. On the other hand, the 
secondary reactions which could occur, namely 

d + "He cere tue and d + T >n + tHe, have higher 
Q@-values (18:4 MeV and 17-58 MeV), and higher peak 
cross-sections, especially the latter which peaks at 
~ 80 keV with o \ 7 b (see Fig.7.3). This latter 
reaction therefore, would, go very quickly at 108 K, 
corresponding to a mean thermal energy of » 10 keV, 
at which energy the cross-section is still as high 
as 10 mb. This mopping up, if complete, would boost 
the energy output to ~ 4 MeV per nucleon. 

From the above, the most convenient reaction to 
solve the problem of obtaining a useful source by 
fusion is undoubtedly d+ T but it may be difficult 
to use it as a basis for large-scale energy pro- 
duction. In principle, the tritium could be regene- 
rated by the reaction 
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n+d-+ T+ Y (6°26 MeV) 


or more conveniently by 
/ 


n=: Oe Te THe + 4°78 MeV. 
This latter reaction has a high thermal cross- 
section, so that Oli in the form of a solid 


compound could form a suitable regenerating 


blanket, as in some fusion bombs. There is no 


T(d,n)*He 


D (d.p) T+D (d.n)He 


3He (d.p)*He 


20 ‘40 60 80. 100 
E,/ keV 


° / Ld 
Fig.7.3. Thermonuclear reaction cross-sections. (Based on 


SEGRE (1964). Nuclei and parttcles, Benjamin, New York). 
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possibility, with these reactions, of breeding 
more tritium than is consumed but, by addition 
of Beryllium into the blanket, neutron 
multiplication can result from the dominant 


reaction 


Re +n>a+ai+ 2n - 1:66 MeV, 


leading to the possibility of tritium breeding. 

A great deal depends upon the supply of materials 
(assuming that the technical problems are solved). 
The abundances in the earth's crust in parts per 
million are, for some of the relevant elements: 


H, 1400(7H, 0-2); Li, 65(°Li, 4-8); Be, 6; 


235 


and compare with U, 4 ( U, 0-03). 


These figures do not give a complete picture, since 
hydrogen is readily available from water (the sea) 

and deuterium separation is probably the simplest of 

all isotope separations. Unless the breeding of 
tritium is relatively simple, it would seem that 
tritium-induced fusion will only be a stepping stone 
towards the final solution based on the d + d reaction. 


PROBLEMS 

7.1. If neutrons (mass m) of kinetic energy Eq are 
isotropically scattered by nuclei of mass M show 
that the scattered neutron energy-spectrum is 
independent of energy within the limits Ey and 


aE: where a = (M-m)/(M+m). Show further that 
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the mean value of lin (E,/E,), where Ey is the 
energy after a single scatter, 1s equal to 
1 + (a ln a)/(1 - a), denoted by &. 


7.2. Extend the above result to the case where the 
initial neutrons are not well defined in energy 
but have a spectral distribution. Hence show 
that, after n collisions, In (£)/E,) =nt. If 
the mean energy of the neutron fission spectrum 
is 2 MeV, show that 18, 25, and 117 collisions 
are required to thermalize neutrons in H, D, 


12 


and C respectively. 


ge and graphite is re- 


7.3. An intimate mixture of 
quired for certain experiments. The graphite is 
known to be contaminated with 1 p.p.m. by weight 
of eR: What is the maximum fraction by weight 
of £92) in the mixture if the multiplication 
factor at infinite size is not to exceed unity? 


= 0°04 b, 3800 b, and 700 b respectively 
10, 235 


(Cabs 

for eG and U at thermal energies; of 
the last cross-section, Oo, accounts for 580 b. 
Assume 2°5 neutrons per fission and that all 


reactions take place at thermal energies only.) 


7.4. For a gas consisting of two types of particles, 
show that the distribution of relative velocities 
of pairs. consisting of one from each type is the 
same as the absolute velocity distribution for a 
Single type of particle but with a mass equal 
to the reduced mass, m m/f (m4 + mo); where m4» 

m, are the masses of the two types. 
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(a) If the fusion reaction 
d+d-+ ° 


deuterons at rest, what is the kinetic energy of. 


He + n + 3-2 MeV takes place with 


the neutron? (b) If the deuterons must come 
within io 72 cm of each other, what energy 

must be supplied to overcome the electrostatic. 
repulsion? (c) If this energy is supplied by 
heating the deuterium to high temperature, what 
order of magnitude of temperature is required? 

(d) Discuss the effect of this heating on the 
energies of the neutrons produced. 

. A nucleus consists of a large number W of particles 
of spin 3}. Determine the number of ways, R(M,); 

a given M, can be produced assuming that M, < 3M 
and that the Pauli Principle can be ignored. By 
identifying the difference R(X) - R(X+1) with the 
number of states of total spin S = X, show that 
this number « (2X+1). Hence the plausibility of 
the statement made on p.173. 


Appendix A: Cross-sections 


To introduce the concept of cross-section, it is con- 
venient to consider classically a nuclear reaction 

in which a beam of particles falls on nuclei which 
behave like completely absorbing spheres of cross- 
sectional area o. Assume, that one such nucleus is con- 
fined somewhere within an aperture of area A through 
which passes a uniform beam of » particles per second. 
The number absorbed per second will obviously be no/A. 
If other nuclei are successively added to build up a 
practical target of thickness t and nuclear number 
density eo (i.e. the number of nuclei per unit volume 
of target material), then pAt nuclei will be exposed 
to the beam, and the number of incident particles 
absorbed per second will be pAtno/A = onpt, provided 
that opt <1. (If this condition does not hold then 
the beam flux is decreasing through the target and all 
nuclei are not exposed to the same beam - under these 
circumstances the beam decreases exponentially through 
the target, in the form n(t) = n(0) exp(-opt)). Notice 
that the result does not depend upon the aperture area 
A since, no matter what this area is, provided that 
the target area is larger, each particle in the beam 
passes through the same target thickness. 

Thus, for the absorption process the number per 
second absorbed is proportional to the cross-sectional 
area of the nucleus. Replacing this completely ab- 
sorbing nucleus by a hard sphere of the same size which 
scatters all projectiles hitting it, the total number 
per second scattered out of the beam will be the number 
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per second previously absorbed, onpt. The number per 
second scattered into solid angle dQ at (86,6), in the 
usual nomenclature, relative to the beam will be 
onotf(68,>)d2, where Sf(8,>)dQ = 1, to ensure that 

the total scattering has the correct value. f(8,6¢) 

is known as the angular distribution of the nuclear 
process. It is convenient to combine the product 
of(6,¢) into a single function do(6,¢)/d&2, which is 
known as the differential cross-section for scattering 
along the direction (6,¢). Most nuclear reactions 

are observed under conditions which are characterized 
by an axis of symmetry in the direction of the beam 

- so, in general, o(6,¢) becomes o(8). Experiments 
may be carried out under conditions in which two axes 
are defined, e.g. the beam axis and an axis of polari- 
zation, and for these the generalized direction (6,9) 
is required. 

These Simple pictures can be extended to cover 
more complex situations: the nucleus may not be com- 
pletely absorbing or scattering but partially the one, 
partially the other, and partially transparent; and, 
following absorption, something is usually emitted - 
the same, or another, type of particle or y-radiation. 
These processes can all be included in a general defi- 
nition of cross-section provided the equivalence of 
cross-section and nuclear cross-sectional area is 
revoked. This equivalence has arisen merely from the 
simple picture used to introduce the concept, and will 
in any case require modification when the wave nature 
of the projectile is taken into account. The fact 
remains, however, that the intrinsic nuclear property 
which determines the yield of a reaction has the 


dimensions of an area. 


APPENDIX A 191 


The general definition of cross-section for the 
reaction A + a-+> B + b is contained in the equation 


dy (8) do (6) 
aaj dQ = not iene 2 
dQ dQ 


age 2 (A.1) 


where d¥, ples/de is the number of particles of type 
b per second emitted into unit solid angle at an 
angle 8 to the beam, and oP a a is the differen- 
tial cross-section for the process (a,b). 

Also of use is the partial cross-section for 


the process (a,b) defined as 


do, pl) 


O = | dQ 
a,b do 


The total cross-section is defined as o =Zo j 
a,tot b a,b 


where the summation is carried out over all possible 
reaction products B + b including the incident channel 
(known as elastic scattering). This latter cross- 
section is very useful when the projectiles are 
neutrons, but not so for charged particles, when it 

is dominated by the comparatively uninteresting but 
very large Coulomb cross-section. 

Finally, notice that cross-section has been de- 
fined in terms of the flux of incoming particles that 
is their density times their velocity. In calculating 
reaction rates using perturbation theory the incoming 
particle wavefunction is normalized to represent one 
particle per unit volume, so in going to the cross- 
section a factor v is needed; ov is a measure of the 
reaction rate as calculated by perturbation theory. 

If the incoming particle is just above threshold, it 


may be possible that the appropriate matrix element 
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varies only slowly with energy, as will also the 
density of final states if the reaction has positive 
@-value. Under these circumstances ov % constant or 
o « l/v. This has been deduced in the text from the 
compound-nucleus formula, but as is seen above it has 


greater generality. 


Appendix B: The need for a pairing term in the 
mass formula 


In order to illustrate the need for a pairing term 

in the mass formula, it is convenient to neglect all 
terms in the mass equation but the Coulomb and symmetry 
terms and to simplify these to the form 


@ (A222) ° * b2 (2-1) (a, b both positive) 


The conditions that (A,Z) is stable relative to its 
neighbours (A, Z-1) and (A, Z+1) is 


a(A-22)* + b2Z(Z-1) < a(A-22-2)7 + b(Z+1)2 
and also 
< a(A-22+2)* + b(Z-1) (2-2) 
These can be re-expressed as 
b.2Z > 4a{A-2Z-1} 
and 


b.2Z < 4a{(A-22Z-1) + (1 + 


Adding a nucleon to (A,Z) gives: 


adding p, the energy change Eo ~a(2A-4Z-1) + bD.22, 


adding n, the energy change £, +a(2A-4Z+1) 


The condition for (A+1,Z) to be stable relative to 
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(A+1,2+1) is that E > Es which simplifies to 


b.2Z > 4a{(A-2Z-1) + 1} 


For the second neutron to produce a stable nucleus, 


the condition is that 
M(A+2,Z) < M(A+2,2Z+1), 


which is equivalent to increasing A to Atl in the last 


inequality, i.e. 
b.2Z > 4a{(A-2Z-1)}+ 2 


These conditions may be expressed in terms of 
an allowable variation of the function b.22/4a as: 


N 


Zn from the 


for (A,Z) stable, the range is 0 > 1 + 
value (A-22Z-1), 


for first neutron, the range is 1 > 1 + ee from the 


Z-1 
value (A-22Z-1), 


N 
“ed 


for second neutron, the range is 27> 1 + from the 


value (A-2Z-1) 


Since W/Z varies from 1 for light nuclei to % 1} 
for the heaviest nuclei, these ranges are, for the 
heavy nuclei 23 units, 1} units, and 3} unit respective- 
ly. Thus the probability of a neutron adding is rough- 
ly = and a further neutron roughly - A similar 
simple argument with protons would lead to é for the 
first proton and zero for the second. Using all the 
terms of the mass equation and the more exact functions 


for the two discussed does not alter the general con- 
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clusion that if (4,2) is stable then the addition 

of two more nucleons is more likely to be of the form 
(p,n) for stability than (2p) or (2n). Thus odd-odd 
nuclei should be as prolific as even-even. 

Obviously a term is missing from the mass equa- 
tion. It must be such that if A and Z are both even 
then, whilst not altering the probability that the 
first particle be a neutron for stability (2 is just 
about correct for a heavy nucleus), it must increase 
the probability for a second neutron giving stability 
up to unity. This extra term is not needed in con- 
Sidering the stability of odd-A nuclei, but must be 
chosen such that in going from even A to A+2, where 
both nuclei are stable, pairs of like particles will 
be added. Hence the form adopted. 


Appendix C: Exchange forces 


To introduce this concept it is interesting to consider 
the single ionized H, molecule - a system of two pro- 


tons and one electron. When the separation 


evat ry 


pj atr, 


Pjatr, 


Fig.C.1. Schematic representation of the hydrogen-molecule ion. 


[R] = |r, - r,| is very large the electron wavefunction 
will look like the hydrogen-atom bound state situated 
at r, or Yr, and designated by o(rz-r,) or o(rz-r,). 
However, the correct wavefunction must be symmetric 

(or antisymmetric) with respect to the two protons. 
Therefore, put w = ay? {6(r5-77) - o(rz-r,)}, where 

the normalization integral a takes the value 2 as 

R +o, whilst a, + 4, a. + 0 as R + 0. The energy 


integral is 
Ate an oe 2 - [vs + var = 
ge Irs-My] 9 [rg-raly 3 
by symmetry. Substituting for ¢ gives four integrals, 


and for ¢ real (corresponding to a bound state) two 
are the same giving (where r = Y,-Ty) 
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aney  - ge = -20yltse2cr)d at + so7(r) pode at # 
0 wy R a —Pr — [r-R] 7 

+ 2f(r-R)o(x)= dt}. (C.1) 


The first term represents the binding energy of 
a hydrogen atom, the second term represents the elec- 
tronic contribution to the Coulomb energy of a proton 
distant R from a hydrogen atom; for large RF this term 
exactly cancels the proton repulsion (1/F on the left- 
hand-side), but for small R it becomes the finite 
binding energy of the hydrogen atom. The last tern, 
unlike the other two, has no classical analogue; it 
effectively arises from not knowing to which proton 
the electron belongs. Its sign depends upon the over- 
lap of the two atomic wavefunctions when separated by 
R and, if the wavefunctions have nodes, can be posi- 
tive or negative. For the lowest s-state there are 
no nodes, so the integral is always positive, giving 
increased binding for the symmetric wavefunction and 
reduced for the antisymmetric. At large distances 
bound wavefunctions fall off exponentially. If the 
s-state looks like exp(-«r) at large r, the main con- 
tribution to the integral will occur nearly halfway 
between the two protons and will obviously contain 
exp(-kR) in the product of the wavefunctions (this 
qualitative argument neglects the effects of terms 
like 1/r in the integral). Thus the integral will 
also behave like exp(-«F) at large distances. 

Now assume that the existence of the electron is 
not known, but that experiments can be performed to 
determine the potential between a proton (charged) and 
a hydrogen atom (neutral) as a function of separation 
R. The classical term will be constant at large R (but 
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unknown since the structure of hydrogen is unknown) ; 
however, at small R, the proton gets inside the elec- 
tron cloud and experiences a net repulsion by the 
other proton. For the symmetric case, the non- 
classical term increases exponentially (approximately) 
to a finite binding at R = 0. The state of affairs is 
presented graphically in Fig.C.2. 


Ordinary’ ( = classical) 


Exchange’ (= non-classical) 
(for symmetric state) 


Fig.C.2. Potential of hydrogen-molecule ion. V=0 corresponds 
to infinite separation between ionized atom (p) and 


neutral atom (H). 


Thus there eixsts a potential between the proton 
and the hydrogen atom which depends upon the spatial 
symmetry of the (p,H) system. During the interaction 
the proton and the hydrogen atom cannot be individually: 
identified, and this is interpreted as exchange of 
charge and spin (in the form of an unknown electron). 
From Fig.C.2 this can result in attraction at large 
R, whilst the repulsion at small R arises because the 
'fundamental-particle' hydrogen atom in fact has a 
structure. 
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The analogue with nucleons is obvious, at least 
for the (p,n) system. A similar argument to the above 
can be applied if n is looked upon as p + m7 and gives 
the correct range for the nuclear force. Since the 
proton and neutron have almost the same mass, the total 
energy (rest + binding) of the bound m must be nearly 
zero. But « = |p|/h in the unphysical region, where 
|p| is given by oe = Be -m2o* with —F = 0. Therefore 
K = m_e/h = (1:4 fm) ~. In order to produce charge- 
independence, a Wo must be introduced (to give the 
(p,p) and (n,n) interactions) and also a 1 (though 
its need is not immediately obvious; it arises because 
the above picture has made the proton look 'more fun- 
damental' than the neutron - the 7 redresses the 
balance since its exchange appears to make the proton 
look like n + 1). 

Having gone through the above reasoning, it is 
convenient to forget the mechanism and express the 
different exchanges mathematically. The operator 
$1°82 = H{(s,+s,)7-s{-s4} = j{s*-s?-s3} takes the value 


2{5(5+1)~-s, (s,+1)-s5(s,+1)} for the spin state S, re- 
sulting from the combination of two spin states s, and 
85+ For 8, = 8 = 4, it takes the value + § for S$ = 1 
and -% for S = 0. 


po zi + 25,-S, = 1 for S = 1, triplet state, and 


Thus the operator 

= -] for S = 0, singlet state. This P° is the spin- 
exchange operator. If tY is used to denote the state 
in which particle 1 has spin up, particle 2 spin down, 
then it can be represented thus: 


Hy = U(MHHt) + ECM 4H) = ZL (S=1) + (5=0)}, 


therefore 
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P°sy = yi (8=1) - (S=0)} = E(t¥4e4) - U(te-¥t) = 4. 


Thus P° has exchanged spins. pe . which exchanges 
charge, can be defined similarly. The general poten- 
tial is then (see Chapter 2, p. 40, for nomenclature) 
Va Vy + Vy pt 4 V, Po + Vy P (where P = -p™ p? 

is a requirement of the Pauli principle). Since the 
exchange operator has a purely numerical part, the 


potential can be rewritten 
PV ep eek Vee ease Cape ea) 


and the operators are (loosely) referred to as exchange 
operators. 

The effect of the operators, p®, Pe: Pp”. is shown 
in Table C.1 for S- and P- states. Notice that it is 
the pattern of the changes of sign which is important. 
The sign of a whole column can be changed merely by 
changing the sign of the potential function accompany- 


ing the operator. 


TABLE C:1. 
ah : t 
WV Symmetry Sign of operator 
Spin Isospin Space W H B M 
3 
Sy S a Ss + = + + 
cn a S Ss + + 2 + 
3 . 
Po12 Ss S a + + + Bs 
Ry a a a + - - 7 
7 


see p.40 for definition of W,H,B,M. 


a denotes 'symmetric'; a denotes ‘antisymmetric'. 


Appendix D: Spin matrices 


Angular momentum can be introduced into quantum 
mechanics by inserting the operator form p = -ihV into 
the classical expression J = rAp, to give J = -ih(raV). 
It is left to the reader to show that JAJ = ihJ (hint: 
use Cartesian coordinates), i.e. M.. My - M, MS ihm, 
and two others, cyclically. 

Although these equations have been derived from 
considerations of orbital angular momentum which is 
quantized to integer values, Dirac has used them to 
show that odd half-integer values are also permitted 
indicating the existence of another type of intrinsic 
angular momentum. Here it is proposed to consider 
J = 3, and to designate it s, the intrinsic spin 
angular momentum of a fundamental particle. For this 
particular case it can be shown that 8) $85 = 0, 
etc. and therefore oy. = (ih/2)s_, etc. Since s, is 
two-valued, having values +3h, the wavefunction re- 
presenting an arbitrary state of spin can be written 
as a two-component vector [5] in an appropriate two- 
dimensional space. Ths components of s, namely, 

8 


s s, will be represented by 2 x 2 matrices, of 


x? y? 
which s. will be diagonal and of the form sh) 91 to 
give the required eigenvalues when acting on the 
eigenfunctions i, ! and (71 representing 'spin up' 

: : = pa a pe 
and 'spin down'. Putting s.. = C qi ey = lo Ar and 


using Oar oe zihs and $8. = zane results in 


a@az=ertdtzh=0,g ic, f = -ib. 
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There is some latitude in the choice of 8 ae the 
forms commonly chosen are 
_ 1,70 1 —~ 1470 -i 
oe zh Ly 0! and ey 2h[; 0! 
Notice that 3? + 3 + s? = s¢ = sn, which conforms 
to the usual expression J™ = n? J(J+1). 
The raising operator 
ol Oe ene ab ‘ _ 01 
h “s =h (6 ,,+18,,) = [o 0! 


obviously operates on ial to give [0] and on a to 
give re hence its name. Similarly the lowering 
operator n 4s” is defined as h”*(s,-is,). 

In a similar way the components of isospin may 
be defined. As for spin, the two-component isospin 
space has no connection with any classical space, 
but, and in contrast to spin, neither does the three- 
component isospin space in which the components 


T oe are defined. There is a one-to-one corres- 


, 
pondciee between the z-component of isospin and the 
charge of the nucleon, but otherwise there is no 

further definition of this space. The isospin raising 
and lowering operators have the property of changing 

one nucleon type to its counterpart, without changing 
any other part of its wavefunction. It should be 

stated that this is a rather restricted view of iso- 
Spin since only nucleons are under consideration as dis- 
tinct from its greater ramifications when dealing with 


other hadrons. 
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Answers to problems 


Des Be. <a eB 
14. € % 0°63, y %¥ 19°1. 
1.5. 3°7 fm. 


: “3 + ieid : 

3.6. A ~~ 16, 4°8, 6°3, 16 fm respectively. 
4 Re Br 1S i:g 

5.1. ~6x 10°28 5, 

6.2. 7-0191 a.m.u. 

6.6. 8-5 ub. 

7.3. 0°00115. 


7.5. (a) 2°4 MeV; (b) 14 keV in Centre-of-mass system; 
(c) ¥ 1-6 x 108 K; (d) thermal spread ~ + 100 keV. 


index 


O-radio activity, ; 
angular momentum barrier, 


barrier penetrability, 
Geiger-Nuttall rule, 


9, 90fFf. 


single-particle reduced width, 


Stability against, 


B-radioactivity, 
forbidden classification, 
Kurie plot, 
parity violation, 
selection rules, 
stability against, 
breeding of nuclear fuel, 
Breit-Wigner formula, 


collective motion, 
rotational, 
vibrational, 

compound nucleus, 

cross-sections, 
1/v law, 


deuteron, 
direct reactions, 
energy dependence, 
Einstein coefficients for 
emission and absorption of 
radiation, 
electron capture, 
stability against, 
electron scattering, 
exchange operators, 


fission, 
asymmetric, 
materials, 
neutron spectrun, 
stability against, 
fusion, 
in Sun, 
in laboratory, 
hydrogen cycle, 
carbon cycle, 
reactions, 


y-radiation, | 
effective charge, 
electric multipoles, 
magnetic multipoles, 
selection rules, 
Weisskopf units, 


143, 


143, 


122, 


internal conversion, 

o-o transitions, 
internal pair creation, 
isomerism and skills, 
isospin, 


kinematics of nuclear 
reactions, 


level density, 
and level width, 
(27+1) factor, 


84, 
173, 


magic numbers, 

mean free path in nuclear 

matter, 

mirror nuclei, 

multipole moments, 
electric monopole, 
electric dipole, 
electric quadrupole, 
magnetic dipole, 


nuclear forces, 
exchange, 
hard core, 
range, 
saturation, 
tensor, 


40, 


nuclear mass, 
semi-empirical formula, 
pairing term, 


nuclear reactions, 
angular correlations, 


compound nucleus, 143, 
direct, 143, 
stripping, 
stripping angular dis- 
tribution, 
nuclear state, angular 
momentum of, 
» parity of 


nucleon-nucleon scattering 


optical model, 
potential parameters, 
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reactor, fast, 171 
,» thermal, 171 
» thermal, control 
of, 176 
shell model, 59ff. 
shell sequences, 67, 68 
spin matrices, 199 
strength function, 87, 152 


thermalization of neutrons,185 


Physical constants and conversion factors 


Avogadro constant 
Bohr magneton 
Bohr radius 
Boltzmann constant 


charge of an electron 

Compton wavelength of electron 
Faraday constant 

fine structure constant 


gas constant 
gravitational constant 
nuclear magneton 
permeability of a vacuum 


permittivity of a vacuum 


Planck constant 
(Planck constant)/2x 


rest mass of electron 

rest mass of proton 
Rydberg constant 

speed of light in a vacuum 


Stefan—Boltzmann constant 
unified atomic mass unit (12C) 
wavelength of a 1 eV photon 


1A = 107!°m; 
O°C = 273-15 K: 


Lor Ny 6022 x 10? mol7! 


Lp 
ao 
k 


e 


9:274x 10°24 JT"! 
5-292 x107~!! m 
1:381x107~77 JK7 


—1:602x 107'°C 


Ac = h/m.c = 2-426 107'? m 


F 


9-649 x 10* C mol™! 


= Uye*c/2h = 7-297 x 1077 (a7? = 137-0) 


€ 


h 
h 
m 
m 


uo] 


8 


R 
c 


8:314JK~! mol7! 
6-673 x 107!1 N m? kg~? 
5-051 x 107?’ J T7! 
4xx 10-7 H m“ exactly 


8-854 x 107'*? Fm™! (1/4ze€9 = 
8-988 x 10° m F~?) 
6-626 x 10734 Js 
1:055 x 10734 Js = 6582x107 *% eVs 


9-110 x 1073! ky = 0-511 MeV/c? 
1:673 x 10°?’ kg = 938-3 MeV/c? 


= pam,e*c?/8h> = 1:097 x 10’ m7! 


2:998 x 108 ms™! 


o = 2nrk*/15h3c? = 5670x107? Wm? K~4 


u 


1 dyne = 107° N; 
1 curie (Ci) = 37x 10!%s7!; 


1J = 10’ erg = 6241 x10'8% eV; 1LeV = 1:602x107'°J; 


1-661 x 107-77 kg = 931°5 MeV/c? 
1:243 x10°° m 


1 gauss (G) = 107“ tesla (T); 


lcaly, = 4-184 J; 


In 10 = 2-303; In x = 2°303 log x; e=2:718; loge = 0:-4343; nx = 3-142 


